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Abstract 

The current status of electroweak precision observables in the Minimal Supersymmetric 
Standard Model (MSSM) is reviewed. We focus in particular on the W boson mass, M^, 
the effective leptonic weak mixing angle, sin^ 9cs, the anomalous magnetic moment of the 
muon, {g — 2)^, and the lightest CP-even MSSM Higgs boson mass, mh- We summarize the 
current experimental situation and the status of the theoretical evaluations. An estimate 
of the current theoretical uncertainties from unknown higher-order corrections and from the 
experimental errors of the input parameters is given. We discuss future prospects for both 
the experimental accuracies and the precision of the theoretical predictions. Confronting 
the precision data with the theory predictions within the unconstrained MSSM and within 
specific SUSY-breaking scenarios, we analyse how well the data are described by the theory. 
The mSUGRA scenario with cosmological constraints yields a very good fit to the data, 
showing a clear preference for a relatively light mass scale of the SUSY particles. The 
constraints on the parameter space from the precision data is discussed, and it is shown that 
the prospective accuracy at the next generation of colliders will enhance the sensitivity of 
the precision tests very significantly. 
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Chapter 1 
Introduction 



1.1 Motivation 

Theories based on Supersymmetry (SUSY) [1] are widely considered as the theoretically 
most appeahng extension of the Standard Model (SM) [2]. They are consistent with the ap- 
proximate unification of the gauge coupling constants at the GUT scale and provide a way 
to cancel the quadratic divergences in the Higgs sector hence stabilizing the huge hierarchy 
between the GUT and the Fermi scales. Furthermore, in SUSY theories the breaking of the 
electroweak symmetry is naturally induced at the Fermi scale, and the lightest super sym- 
metric particle can be neutral, weakly interacting and absolutely stable, providing therefore 
a natural solution for the dark matter problem. 

SUSY predicts the existence of scalar partners fi, Jr to each SM chiral fermion, and spin- 
1/2 partners to the gauge bosons and to the scalar Higgs bosons. So far, the direct search 
for SUSY particles has not been successful. One can only set lower bounds of O(IOO) GeV 
on their masses [3]. The search reach will be extended in various ways in the ongoing 
Run II at the upgraded Fermilab Tevatron [4]. The LHC [5,6] and the e~^e~ International 
Linear Collider (ILC) [7-9] have very good prospects for exploring SUSY at the TeV scale, 
which is favoured from naturalness arguments. From the interplay of both machines detailed 
information on the SUSY spectrum can be expected in this case [10]. 

In the Minimal Supersymmetric extension of the Standard Model (MSSM) two Higgs 
doublets are required resulting in five physical Higgs bosons [11]. The direct search resulted 
in lower limits of about 90 GeV for the neutral Higgs bosons and about 80 GeV for the 
charged ones [12,13]. The Higgs search at the Tevatron will be able to probe significant parts 
of the MSSM parameter space at the 95% C.L. even with rather moderate luminosity [14]. 
The LHC will discover at least one MSSM Higgs boson over most of the MSSM parameter 
space [5,6,15-17]. The ILC will be able to detect any Higgs boson that couples to the Z boson 
in a decay-mode independent way. The properties of all Higgs-bosons which are within the 
kinematic reach of the ILC will be determined with high precision [7-9]. 

Contrary to the SM case, where the mass of the Higgs boson is a free parameter, within the 
MSSM the quartic couplings of the Higgs potential are fixed in terms of the gauge couplings 
as a consequence of SUSY [11]. Thus, at the tree- level, the Higgs sector is determined 
by just two independent parameters besides the SM electroweak gauge couplings g and g', 
conventionally chosen as tan/3 = V2/V1, the ratio of the vacuum expectation values of the 
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two Higgs doublets, and M^, the mass of the CP-odd A boson. As a consequence, the mass 
of the hghtest CP-even MSSM Higgs boson can be predicted in terms of the other model 
parameters. 

Besides the direct detection of SUSY particles and Higgs bosons, SUSY can also be probed 
via the virtual effects of the additional particles to precision observables. This requires a 
very high precision of the experimental results as well as of the theoretical predictions. The 
wealth of high-precision measurements carried out at LEP, SLC and the Tevatron [18] as well 
as the "Muon g — 2 Experiment" (E821) [19] and further low-energy experiments provide 
a powerful tool for testing the electroweak theory and probing indirect effects of SUSY 
particles. The most relevant electroweak precision observables (EWPO) in this context are 
the W boson mass, M^, the effective leptonic weak mixing angle, sin^^eff) the anomalous 
magnetic moment of the muon, = [g — 2)^/2, and the mass of the lightest CP-even MSSM 
Higgs boson, rrih. While the current exclusion bounds on nih already allow to constrain the 
MSSM parameter space, the prospective accuracy for the measurement of the mass of a light 
Higgs boson at the LHC of about 200 MeV [5,6] or at the ILC of even 50 MeV [7-9] could 
promote ruh to a precision observable. Owing to the sensitive dependence of rrih on especially 
the scalar top sector, the measured value of rrih will allow to set stringent constraints on the 
parameters in this sector. 

Since the experimental data — with few exceptions — are well described by the SM [18], 
the electroweak precision tests at present mainly yield constraints on possible extensions of 
the SM, e.g. lower limits on SUSY particle masses. Nevertheless, one can use the available 
data to investigate whether small deviations from the SM predictions could be caused by 
quantum effects of the SUSY particles: sleptons, squarks, gluinos, charginos/neutralinos and 
additional Higgs bosons, and what regions of the SUSY parameter space might be favoured. 

1.2 The structure of the MSSM 

The MSSM constitutes the minimal supersymmetric extension of the SM. The number of 
SUSY generators is = 1, the smallest possible value. In order to keep anomaly cancellation, 
contrary to the SM a second Higgs doublet is needed [20]. One Higgs doublet, TCi, gives 
mass to the d-type fermions (with weak isospin -1/2), the other doublet, TC2, gives mass to 
the M-type fermions (with weak isospin +1/2). All SM multiplets, including the two Higgs 
doublets (2HDM), are extended to supersymmetric multiplets, resulting in scalar partners 
for quarks and leptons ( "squarks" and "sleptons" ) and fermionic partners for the SM gauge 
boson and the Higgs bosons ("gauginos" and "gluinos"). In Tab. II. ll the spectrum of the 
MSSM fields is summarized (family indices are suppressed). In this report we do not consider 
effects of complex phases, i.e. we treat all MSSM parameters as real. 

The mass eigenstates of the gauginos are linear combinations of these fields, denoted as 
"neutralinos" and "charginos". Also the left- and right-handed squarks (and sleptons) can 
mix, yielding the mass eigenstates (denoted by the indices 1, 2 instead of L, R). All physical 
particles of the MSSM are given in Tab. 11.21 
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superfield {SU{3), SU{2),U{1)) 


2HDM particle spin 


SUSY partner spin 


Q (3, 2, I) 
U (3*, 1, -|) 
D (3*, 1, 1) 


iu,d)L 1 
Ur \ 
dR \ 


{u,d)L 

U*R 

4 


L (1, 2, -1) 
^ (1, 1, 2) 


eR \ 


= (z>, e)L 


^1 (1, 2, -1) 
H2 (1, 2, 1) 


{HIH^)l 
{m.Hl)L 




W (1, 3, 0) 
S (1, 1, 0) 


W' 1 
5° 1 


5° \ 


Ga (8, 1, 0) 




9a \ 



Table 1.1: Superfields and particle content of the MSSM. 



1.2.1 The Higgs sector of the MSSM 

The two Higgs doublets form the Higgs potential [11] 

Y _ I |,.|2Maj' |2 I 1^2 I |,.|2 



K + + K + \m\'H2V - mUea^mnl + h.c.) 

+\{9i' + 92') [m^ -m'Y + \92'\n\n2\\ 



(1.1) 



which contains mi,m2,mi2 as soft SUSY breaking parameters and /i as the Higgsino mass 
parameter; g,g' are the SU{2) and U{1) gauge couplings, and ei2 = —1. 
The doublet fields Tii and 0.2 are decomposed in the following way: 



Hi 

n2 



Hi 



V2 

^2 + ^(0^ + ZX^) 



(1.2) 



The potential ()1.1|1 can be described with the help of two independent parameters (besides 
g and g'): tan (3 = V2/V1 and M\ = — m^gl^^^/^ + cot (5 ), where is the mass of the 
CP-odd A boson. 

The diagonalization of the bilinear part of the Higgs potential, i.e. the Higgs mass ma- 
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2HDM particle 


spin 


SUSY particle 


spin 


quarks: 


q 


1 

2 


squarks: 







IfiDtons" 


/ 


1 
2 


slentons" 







gluons: 


9a 


1 


Gluinos: 




1 

2 


gauge bosons: 




1 


Neutralinos: 


Xi) X2-I Xsi Xa 


1 

2 


Higgs bosons: 







Charginos: 


xt, xt 


1 
2 



Table 1.2: The particle content of the MSSM. 



trices, is performed via the orthogonal transformations 

\ / cos a sin a 

/i° / \ — sin a cos a 

_ ( cos/3 sin/3 

A" 
i7± 



— sin /3 cos /3 
cos /3 sin /3 

— sin /? cos /? 



The mixing angle a is determined through 



tan 2q; = tan 2/3 



One gets the following Higgs spectrum: 

2 neutral bosons, CP = +1 
1 neutral boson, CV = — 1 
2 charged bosons 
3 unphysical Goldstone bosons 



0? 
X? 

0f 



■- < a < . 
2 



A" 

G^,G+,G- 



The masses of the gauge bosons are given in analogy to the SM: 



- 2/ 2 I 2\ 



= 0. 



(1.3) 
(1.4) 
(1.5) 



(1.6) 



;i.7) 



(1.8) 



At tree level the mass matrix of the neutral CP-even Higgs bosons is given in the 0i-02- 
basis in terms of Mz, Ma, and tan/3 by 



yi^2,tree 
Higgs 



mz 



mz 



m 



•1<P2 
2 

(j>2 
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Ml sin^ 13 +Ml cos^ {3 -{Ml + M|) sin /? cos/? \ , . 

-(Ml + M|) sin/? cos/? Ml cos2 /3 +M|sin2/3 ^^■^> 

which by diagonahzation according to eq. ()1.3|1 yields the tree-level Higgs boson masses 



^h,tree 



^ ( ""r^^ ^2 ] ■ (1.10) 



The mixing angle a satisfies 



Ml + M| vr 



tan 2a = tan 2/3— — |, < a < 0. (1.11) 



Ml - M| ' 2 



Since we treat all MSSM parameters as real there is no mixing between CP-even and CP-odd 
Higgs bosons. 

The tree-level results for the neutral CP-even Higgs-boson masses of the MSSM read 



m 



H,h 



Ml + M'z± \ (Ml + M|)2 - AM^Ml cos^ 2/? 



1.121 



This implies an upper bound of tree < Mz for the light CP-even Higgs-boson mass of the 
MSSM. For a discussion of large higher-order corrections to this bound, see Sect. 12.71 The 
direct prediction of an upper bound for the mass of the light CP-even Higgs-boson mass is 
one of the most striking phenomenological predictions of the MSSM. The existence of such 
a bound, which does not occur in the case of the SM Higgs boson, can be related to the fact 
that the quartic term in the Higgs potential of the MSSM is given in terms of the gauge 
couphngs, while the quartic couphng is a free parameter in the SM. 



1.2.2 The scalar quark sector of the MSSM 

The squark mass term of the MSSM Lagrangian is given by 

1 




^m, = -^[flPn)Z [ y ] , (1.13) 



where 

M? + M| cos 2/3 (// - Qfsl^r) + mj mf{Af - /i{cot /5 ; tan/? }) 
mf{Af - /i{cot/5;tan/?}) M? + M| cos 2/5 Q/s^ + 



;i.l4) 



and {cot P ;tan/5 } corresponds to {u; dj-type squarks. The soft SUSY breaking term Mg, 
is given by: 

_ ( Mfj for right handed u-type squarks ,^ . 

^' y Mjj for right handed d-type squarks ' \ ■ ) 

In order to diagonalize the mass matrix and to determine the physical mass eigenstates the 
following rotation has to be performed: 



fi \ ^ { cos^^- sin^j \ ( ^ 
h J V -sin^/ cos Of J \ fn 



1.16) 
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The mixing angle is given for tan /3 > 1 by: 



cos 6? 



sin 6? 



JR h 



^ mj (Af - fi{cot ; tan/5 })2 + 



m : 



=F sgn 



Af — yu{cot P ; tan P } 



mj {Af — ;u{cot (3 ; tan (3 }y 



^ mj {Af - fi{cot ; tan/3 }y + {mj^ 



m 



2 \2 



;i.i7) 



(1.18) 



The negative sign in ()1.18p corresponds to w-type squarks, the positive sign to d-tjpe ones. 
m| = M? + M| cos 2/? QfS^r + m'l denotes the lower right entry in the squark mass 

JR 



matrix ()1.14p . The masses are given by the eigenvalues of the mass matrix: 



m 



h,2 



+ -M|cos2/3 J| +m 
2 



;i.i9) 



sgn 



2 



M? - M?, + M| cos 2/3 (// - 2Q /s^^,)] ' + Am} (a^ - /i cot /3 )' 
M? - M?, + M| cos 2f3 (// - 2Qjs2^) ^ + 4mj (a^ - /i tan/? 



M|, + M| cos 2/3 (// - 2QfS 



w) 



for M-type and d-type squarks, respectively. Since the non-diagonal entry of the mass matrix 
eq. ()1.14j) is proportional to the fermion mass, mixing becomes particularly important for 
f = t, in the case of tan (3^1 also for f = b. 

For later purposes it is convenient to express the squark mass matrix in terms of the 
physical masses mf_^,mj^ and the mixing angle dj: 



cos2 ^,m| + sin^ ^.m^- sin % cos ^(ml - m 

sin 9f cos Of{m% — m% ) siv? 6fm% + cos^ Orm 
J J ^ Jl J2' J h J . 



f2' 
2 

h 



;i.2o) 



Af can be written as follows: 



Ai 



sin Of cos 6f{rn?j: — m ^ ] 
J I ^ Jl rr_ 

rrif 



+ /i{cot /3 ; tan /3 }. 



1.21) 



1.2.3 Charginos 

The charginos xt ~ l'^) are four component Dirac fermions. The mass eigenstates are 
obtained from the winos and the charged higgsinos H{ , 



w- 



-iX- 



w~ 



'1.221 
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The chargino masses are defined as mass eigenvalues of tlie diagonalized mass matrix, 



or given in terms of two- component fields 
where X is given by 



;i.24) 



M2 V2Mwsm/3 
X = I ' 1 • (1.25) 

V2Mw cos (3 jj, 

In the mass matrix M2 is the soft SUSY-breaking parameter for the Majorana mass term. 
H is the Higgsino mass parameter from the Higgs potential eq. (jl.lj) . 

The physical (two-component) mass eigenstates are obtained via unitary (2 x 2) matrices 
U and V: 



This results in a four- component Dirac spinor 



z,j = l,2. (1.26) 



Xt=[^-- ] ^ = 1,2, (1.27) 



where U and V are given by 



with 



0+ det X > 
U = 0_ ; V = <^ (1.28) 

a30+ detX<0 



_ , COS0+ sm0-|- , 

0± = 7^ ] ; (1.29) 



— sm (p± cos (p± 
cos 0-1- und sin 0-1- are given by (e = sgn[detX]) 

Mw{smP m-+ + e cos (3 m-+] 

iim& — 



(M2 + e/im^+) 

— u m -+ — e M2 m -+ 

tan0_ = — ^2 . (1.30) 

V 2 Mw (sin /3 m -+ + e cos /3 m -+ ) 

(If 0-1- < it has to be replaced by 0-(- + vr.) m-+ and m-+ are the eigenvalues of the 
diagonalized matrix 

MLg,x+ = vxtxv"' = U*XXt(U*)"' 
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m 



(1.31) 



They are given by 



xr2 



+ 4M^ cos^ 2P + 4M^(M2^ + /^^ + 2 M2 sin 2/5 ) 



;i.32) 



1.2.4 Neutralinos 



Neutralinos x° (i = 1,2,3,4) are four- component Majorana fermions. They are the mass 
eigenstates of the photino, 7, the zino, Z, and the neutral higgsinos, and H2, with 
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-iXz 
iXz 



H2 



;i.33) 



Analogously to the SM, the photino and zino are mixed states from the bino, B, and the 
wino, W, 



with 



7 



SH' + 5 Cm/ 

W^cw - Bsw ■ 



The mass term in the Lagrange density is given by 

with the two-component fermion fields 
The mass matrix Y is given by 



;i.35) 



(1.36) 



(1.37) 



Y = 



/Ml -MzSw cos f3 MzSw sin P \ 

M2 Mzcw cos P -Mzcw sin P 

—MzSw cos P MzCw cos P —/i 

\ MzSw sin P —MzCw sin P —/i / 



(1.38) 



The physical neutralino mass eigenstates are obtained with the unitary transformation ma- 
trix N: 

XH^V'? i,j = l,...,4, (1.39) 
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resulting in the four-component spinor (representing the mass eigenstate) 

X°=(^I) i = l,...,4. (1.40) 

The diagonal mass matrix is then given by 

Mdiag,xO = N*YN-i . (1.41) 

1.2.5 Gluinos 

The gluino, g, is the spin 1/2 superpartner (Majorana fermion) of the gluon. According to 
the 8 generators of SU{3)c (colour octet), there are 8 gluinos, all having the same Majorana 
mass 

rUg^Ms. (1.42) 
In SUSY GUTs Mi, M2 and M3 are not independent but connected via 

m~g^M3^^M2 = — s^M2, Mi = ^4^M2. (1.43) 

1.2.6 Non-minimal flavour violation 

The most general flavour structure of the soft SUSY-breaking sector with flavour non- 
diagonal terms would induce large flavour-changing neutral-currents, contradicting the exper- 
imental results [3]. Attempts to avoid this kind of problem include flavour- diagonal SUSY- 
breaking scenarios, like minimal Supergravity (with universality assumptions) or gauge- 
mediated SUSY-breaking, see the next subsection. In these scenarios, the sfermion-mass 
matrices are flavour diagonal in the same basis as the quark matrices at the SUSY-breaking 
scale. However, a certain amount of flavour mixing is generated due to the renormalization- 
group evolution from the SUSY-breaking scale down to the electroweak scale. Estimates of 
this radiatively induced off-diagonal squark-mass terms indicate that the largest entries are 
those connected to the SUSY partners of the left-handed quarks [21,22], generically denoted 
as A.LL. Those off-diagonal soft SUSY-breaking terms scale with the square of diagonal soft 
SUSY-breaking masses Msusy, whereas the A^r and A^l terms scale linearly, and Arr with 
zero power of Msusy- Therefore, usually the hierarchy ^ Al^j^l ^ A^i^ is realized. It 
was also shown in Refs. [21,22] that mixing between the third and second generation squarks 
can be numerically significant due to the involved third-generation Yukawa couplings. On 
the other hand, there are strong experimental bounds on squark mixing involving the first 
generation, coming from data on K^-K^ and D^-D^ mixing [23,24]. 

Considering the scalar quark sector with non-minimal fiavour violation (NMFV) for the 
second and third generation, the squark mass matrices in the basis of {cL,tL,CR,tR) and 
{sl, bL, sr, Ir) are given by 



/ 


M? 


A* 


rricXc 


A* 

^LR 


\ 




A* 


M? 


A* 


mtXt 






rricXc 


A* 

^RL 




A* 

^RR 




V 




mtXt 


A* 

^RR 







(1.44) 
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with 



^1 



■l-'s 



A^ 



A^ 

^RR 



rribXb 

Kr 



(1.45) 



Mi 

rCq 



Ml 



m 
m 



w) 



COS2/3 M|(T| - QgS 
COS2I3 MlQgsl, {q = t,c) 



+ cos 2(3 M'zQqsiy {q = b, s) 
-//(tan/3)-2^3 



(1.46) 



where rUg, Qg and T| are the mass, electric charge and weak isospin of the quark q. Mq^, Mfj^, 
Mj^^ are the soft SUSY-breaking parameters. The SU{2) structure of the model requires 

Mq to be equal for t and b as well as for c and s. 

In order to diagonalize the two 4x4 squark mass matrices, two 4x4 rotation matrices, 
Ru and R^, are needed. 



h 

cr 

\iRj 



a 



( 5l \ 

bL 
Sr 

\bRj 



(1.47) 



yielding the diagonal mass-squared matrices as follows. 



diag{m|,m|,m|,m|J"'^ 



diag{m?^,m?^,m?3,m?J 

4 

For the numerical analysis we use 

Ai^ = AM^ M^^ , 
Ati = \M-^Mg^ , 



a,l3 



= Rf (Mj). . (i?;?'^y 

d V d'l.'i ^ d ' 



aL = a^^ = a^^ = o. 



(1.48) 
(1.49) 



(1.50) 



Feynman rules that involve two scalar quarks can be obtained from the rules given in the 
fh-i /r basis by applying the corresponding rotation matrix {q — u,d), 



V{XUp) = Rf' Rif V{XM-) . 



;i.5i) 



Thereby V{Xqiqj) denotes a generic vertex in the /l, /r basis, and V{Xqaq'p) is the vertex 
in the NMFV mass-cigcnstatc basis. The Feynman rules for the vertices needed for our 
applications, i.e. the interaction of one and two Higgs or gauge bosons with two squarks, can 
be found in Ref. [25]. 
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1.2.7 Unconstrained MSSM versus specific models for soft SUSY 
breaking 

In the unconstrained MSSM no specific assumptions are made about the underlying SUSY- 
breaking mechanism, and a parametrization of all possible soft SUSY-breaking terms is used 
that do not alter the relation between the dimensionless couplings (which ensures that the 
absence of quadratic divergences is maintained). This parametrization has the advantage of 
being very general, but the disadvantage of introducing more than 100 new parameters in 
addition to the SM. While in principle these parameters (masses, mixing angles, complex 
phases) could be chosen independently of each other, experimental constraints from flavour- 
changing neutral currents, electric dipole moments, etc. seem to favour a certain degree of 
universality among the soft SUSY-breaking parameters. 

Within a specific SUSY-breaking scenario, the soft SUSY-breaking terms can be pre- 
dicted from a small set of input parameters. The most prominent scenarios in the literature 
are minimal Supergravity (mSUGRA) [26,27], minimal Gauge Mediated SUSY Breaking 
(mGMSB) [28] and minimal Anomaly Mediated SUSY Breaking (mAMSB) [29-31]. The 
mSUGRA and mGMSB scenarios have four parameters and a sign, while the mAMSB sce- 
nario can be specified in terms of three parameters and a sign. 

Detailed experimental analyses within the multi-dimensional parameter space of the un- 
constrained MSSM would clearly be very involved. Therefore one often restricts to certain 
benchmark scenarios, see e.g. Refs. [32-35], or relies on underlying assumptions of a specific 
SUSY-breaking scenario. 

The EWPO can be analyzed within the unconstrained MSSM (or extensions of it), which 
allows to set constraints on the SUSY parameter space in a rather general way. In our numer- 
ical anaysis in chapter IHl we discuss the impact of EWPO in the context of the unconstrained 
MSSM, while in chapter |3] we focus on the mSUGRA, mGMSB and mAMSB scenarios as 
special cases. 

1.2.8 Experimental bounds on SUSY particles 

The non-observation of SUSY particles at the collider experiments carried out so far place 
lower bounds on the masses of SUSY particles which are typically of (9(100 GeV) [3]. These 
bounds, however, depend on certain assumptions on the SUSY parameter space, for instance 
on the couplings and decay characteristics of the particles or the validity of a certain SUSY- 
breaking scenario. 

Relaxing some of these assumptions can result in bounds that are much weaker than 
the ones that are usually quoted. As an example, collider experiments do not provide any 
lower bound on the mass of the lightest neutralino if the GUT relation connecting Mi and 
M2, see eq. ()1.43|) . is lifted [36]. It is interesting to investigate in how far the results for 
EWPO can narrow down the parameter space where the bounds from direct searches are 
very weak. Such an analysis has been carried out, for instance, for a scenario with a light 
scalar bottom quark of 0{5 GeV). In Ref. [37] it has been shown that a light scalar bottom 
quark is consistent with the constraints from the EWPO and the LEP Higgs search. 
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1.3 Electroweak precision observables 



In general there are two possibilities for virtual effects of SUSY particles to be large enough 
to be detected at present and (near future) experiments. On the one hand, these are rare 
processes, where SUSY loop contributions do not compete with a large SM tree-level contri- 
bution. Examples are rare b decays like b — > S7, Bg — > lJ''^fJ'~, and electric dipole moments 
(EDMs). For processes of this kind the SUSY prediction for the rates can be much larger 
(sometimes by orders of magnitude) than the SM one. 

On the other hand, EWPO which are known with an accuracy at the per cent level or 
better have the potential to allow a discrimination between quantum effects of the SM and 
SUSY models. Examples are the W boson mass, Mw, and the Z-boson observables, like the 
effective Icptonic weak mixing angle, sin^^eff- 

This distinction between rare processes and EWPO is of course not a completely rigid 
one. The anomalous magnetic moment of the muon, for instance, corresponds both to a 
rare process according to the above definition and to an EWPO which has been measured 
with high accuracy. In view of the prospects for precision measurements of the mass of the 
lightest CP-even Higgs boson, nih, at the next generation of colliders, we also treat nih as 
an EWPO. 

In the present report we concentrate our discussion on EWPO, in particular the observ- 
ables in the W- and Z-boson sector, the anomalous magnetic moment of the muon, and 
the mass of the lightest CP-even Higgs boson. We just briefly comment on rare processes 
in the following section and occasionally in our numerical discussion. For a more thorough 
investigation of the constraints on the SUSY parameter space we refer to the literature. 
For reviews of rare decays see Ref. [38], results for EDMs in the MSSM can be found in 
Refs. [39,40] and in references therein. 

1.3.1 Constraints on the SUSY parameter space from rare pro- 
cesses 

The branching ratio BR(6 57) receives, besides the SM loop contribution involving 
the W boson and the top quark, additional contributions from chargino/stop and charged 
Higgs/stop loops [41]. The SUSY contributions are particularly large for hght charged Higgs 
bosons and large fi or tan/?. The currently available SUSY contributions to BR(b 57) 
include the one-loop result and leading higher-order corrections. The comparison of the the- 
ory prediction with the data imposes important constraints on the parameter space both of 
general two-Higgs-doublet models and of the MSSM. In the latter case it is possible that the 
two kinds of additional contributions are individually large but interfere destructively with 
each other, resulting in only a small deviation of the decay rate from the SM prediction. 

Another interesting channel is the decay Bs fi'^fi~. The SM contribution to this decay 
is tiny, resulting in a BR of about 10~^ [42]. Within SUSY, however, diagrams enhanced by 
tan (3^ can contribute. Thus the decay width can grow with tan^ (3 and the BR can be much 
larger than in the SM [43] , see Ref. [44] for a recent review. The available corrections in the 
MSSM consist of the full one-loop evalution and the leading two-loop QCD corrections. The 
current bound from the Tevatron is BR(B,. H'^li') < 2.7 x 10"^ at the 90% C.L. [45]. A 
substantial improvement in this bound can be expected in the forthcoming years. 
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A different way for probing SUSY is via its contribution to EDMs of heavy quarks [46] , of 
the electron and the neutron (see Refs. [40,47] and references therein), or of deuterium [48]. 
While SM contributions start only at the three-loop level, due to its complex phases the 
MSSM can contribute already at one-loop order. Also the leading two-loop corrections for 
the electron and neutron EDMs are available. Large phases in the first two generations of 
(s)fermions can only be accomodated if these generations are assumed to be very heavy [49] 
or large cancellations occur [50], see however the discussion in Ref. [39]. 

1.3.2 Pseudo-observables versus realistic observables 

The quantities that can be directly measured in experiments are cross sections, line shape 
observables, forward-backward asymmetries etc., deemed "realistic observables" in the lan- 
guage of Rcf. [51]. The obtained results depend on the specific set of experimental cuts 
that have been applied and are infiuenced by detector effects and other details of the ex- 
perimental setup. In order to determine quantities like masses, partial widths or couplings 
from the primarily measured observables, a deconvolutiuon (unfolding) procedure is ap- 
phed. This procedure involves manipulations hke unfolding the QED corrections, subtracting 
photon-exchange and interference terms, subtracting box-diagram contributions, unfolding 
higher-order QCD corrections, etc. These secondary quantities are therefore called "pseudo- 
observables" in Ref. [51]. 

The procedure of going from realistic observables to pseudo-observables results in a slight 
model dependence of the pseudo-observables. As an example, the experimental value of the 
Z-boson mass has a slight dependence on the value of the Higgs-boson mass in the SM, see 
Refs. [18,52]. The EWPO on which we focus in this report arc pseudo-observables in the 
sense outlined above. At the level of electroweak precision physics, it is important to keep 
in mind that in order to obtain the numerical values of the EWPO given in the literature 
the Standard Model has been used in several steps for calculating the subtraction terms. An 
obvious model dependence also occurs if, instead of performing an explicit subtraction of SM 
terms, parameters like M^, as{Mz), etc. are determined directly from a SM fit, containing 
the full set of SM corrections, to the realistic observables. 

Using the same numerical values of the EWPO as input for analyses within the MSSM 
(or other extensions of the SM) is obviously only justified if new physics contributions to 
the subtraction terms and the implemented higher-order corrections are negligible. As an 
example, the experimental value extracted for as{Mz) in the MSSM (for a given SUSY mass 
spectrum) would somewhat differ from the SM value of as{Mz). 

A consistent treatment of the model dependence of the EWPO is necessary in a precision 
analysis of the MSSM. At the current level of experimental precision the shift induced in 
the EWPO from taking into account the full MSSM particle content instead of the SM will 
normally be of minor importance. In some regions of the parameter space, in particular where 
some of the SUSY particles are very light, an explicit verification of the above assumption 
would however be desirable. 

Concerning the determination of the MSSM parameters, additional comphcations arise 
compared to the SM case. In general the model dependence is relatively small for masses, 
since the mass of a particle can closely be related to one particular realistic observable. For 
couplings (with the exception of the electromagnetic coupling in the Thomson limit), mixing 
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angles, etc., on the other hand, the model dependence is relatively large. In contrast to the 
SM, many of the MSSM parameters are not closely related to one particular observable, e.g. 
tan/3, 11, the stop and sbottom mixing angles, complex phases, etc., resulting in a relatively 
large model dependence. Therefore, the approach of extracting pseudo-observables with only 
a fairly small model dependence seems not to be transferable to the case of the MSSM. It 
seems that eventually the MSSM parameters will have to be determined in a global fit of 
the MSSM to a large set of observables, taking into account higher-order corrections within 
the MSSM, see Refs. [53,54] or Ref. [55] for an attempt of a coordinated effort. 

1.3.3 EWPO versus effective parameters 

In this report we focus our discussion on the EWPO, i.e. (pseudo-)observables like the W- 
boson mass, Mw, the effective leptonic weak mixing angle, sin^ 9cs, the the leptonic width of 
the Z boson, F;, the anomalous magnetic moment of the muon, = {g — 2)^/2, the mass of 
the lightest CP-even MSSM Higgs boson, nih, etc. In the literature virtual effects of SUSY 
particles are often discussed in terms of effective parameters instead of the EWPO (see e.g. 
Ref. [56] and references therein). We do not follow this approach, and just briefly comment 
about it in the following. 

Since for the accuracies anticipated at future colliders, see Tab. 11.41 below, it is particu- 
larly important to have a precise understanding of how effects of new physics can be probed 
in a sensible way, the virtues and range of applicability of effective parameters need to be 
assessed. 

A widely uses set of parameters are the S, T, U parameters [57]. They are defined 
such that they describe the effects of new physics contributions that enter only via vacuum- 
polarization effects (i.e. self-energy corrections) to the vector boson propagators of the SM 
(i.e. the new physics contributions are assumed to have negligible couplings to SM fermions). 
The S, T, U parameters can be computed in different models of new physics as certain combi- 
nations of one-loop self-energies. Experimentally, their values are determined by comparing 
the measured values Af^^ of a number of observables with their values predicted by the SM, 
Af^, i.e. Af^^ = Af^ + ff'^{S,T, U). Here A^^ contains all known radiative corrections in 
the SM, while ff^{S, T, U) is a (linear) function of the parameters 5, T, U and describes the 
contributions of new physics. The SM prediction A^^ is evaluated for a reference value of mt 
and Mh- For most precision observables the corrections caused by a variation of mt and Mh 
at one-loop order can also be absorbed into the parameters S", T, and U . A non-zero result 
for S, T, U determined in this way indicates non-vanishing contributions of new physics 
(with respect to the SM reference value). 

From their definition, it is obvious that the S, T, U parameters can only be applied for 
parameterizing effects of physics beyond the SM. Taking into account the full contributions 
within the SM cannot be avoided, as these contributions (containing also vertex and box 
corrections) cannot consistently be absorbed into the S, T, U parameters (for a more detailed 
discussion of this point, see Ref. [58]). 

Examples of new physics contributions that can be described in the framework of the S*, 
T, U parameters are contributions from a fourth generation of heavy fermions or effects from 
scalar quark loops to the W- and Z-boson observables. A counter example going beyond 
the S", T, U framework are SUSY corrections to the anomalous magnetic moment of the 
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muon. According to their definition, the S, T, U parameters are restricted to leading order 
contributions of new physics. They should therefore be apphed only for the description of 
small deviations from the SM predictions, for which a restriction to the leading order is 
permissible. It appears to be questionable, on the other hand, to apply them to cases of very 
large deviations from the SM, like extensions of the SM with a very heavy Higgs boson in 
the range of several TeV. 

Other parameterizations have been suggested (see e.g. Refs. [59,60]) with no reference 
to the SM contribution and which are not restricted in the possible kinds of new physics. 
These parameterizations are defined as certain linear combinations of different observables. 
It is however not in all cases obvious that studying the experimental values and the theory 
predictions for these parameters is of advantage compared to studying the EWPO themselves. 
For a recent discussion of effective parameters, see also Ref. [61]. 

1.3.4 Current experimental status of EWPO 

LEP, SLC, the Tevatron, and low-energy experiments have collected an enormous amount 
of data on EWPO. Examples for the current experimental status of EWPO are given in 
Tab. II. H| including their relative experimental precision. The quantities in the first three 
lines, Mz, Gp, and rut, are usually employed as input parameters for the theoretical pre- 
dictions. The observables M^, sin^6'efT, Tz, on the other hand, are used for testing the 
electroweak theory by comparing the experimental results with the theory predictions. Com- 
paring the typical size of electroweak quantum effects, which is at the per cent level, with the 
relative accuracies in Tab. II. 3( which are at the per mille level, clearly shows the sensitivity 
of the electroweak precision data to loop effects. 





central value 


absolute error 


relative error 


Mz [GeV] 


91.1875 


±0.0021 


±0.002% 


G, [GeV-2] 


1.16637 X 10"^ 


±0.00001 X 10"^ 


±0.0009% 


mt [GeV] 


178.0 


±4.3 


±2.4% 


Mw [GeV] 


80.425 


±0.034 


±0.04% 


sin^ e^s 


0.23150 


±0.00016 


±0.07% 


Tz [GeV] 


2.4952 


±0.0023 


±0.09% 



Table 1.3: Examples of EWPO with their current absolute and relative experimental er- 
rors [3,18]. 

The experimental accuracy of the precision observables will further be improved at the 
currently ongoing Run II of the Tevatron, the LHC and a future ILC, with the possible 
option of a high luminosity low-energy run, GigaZ [7-9,62]. The most significant improve- 
ments among the EWPO can be expected for M]y and sin^ ^efr- If the Higgs boson will be 
detected, a precise measurement of its mass will be important for testing the electroweak 
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theory. Concerning the input parameters, the experimental error of the top-quark mass is 
the dominant source of theoretical uncertainty in electroweak precision tests. This will re- 
main to be the case even with the accuracy on rrit reachable at the LHC [63]. Thus, the 
high-precision measurement of rrit at the ILC will be crucial for an increased sensitivity to 
virtual effects of new physics [63,64]. 

The prospective accuracy for M]y, sin^ 9es, rrit and nih (for a value of ^ 120 GeV) at 
the Tevatron, at the LHC (combined with the data collected at the Tevatron) and the ILC 
(with and witout GigaZ option) are summarized in Tab. 11.41 (see Ref. [65] and references 
therein) . 





now 


Tevatron 


LHC 


ILC 


ILC with GigaZ 


(5sin2 ^eff(xlO^) 


16 




14-20 




1.3 


6Mw [MeV] 


34 


20 


15 


10 


7 


6mt [GeV] 


4.3 


2.5 


1.5 


0.2 


0.1 


6mh [MeV] 






200 


50 


50 



Table 1.4: Current and anticipated future experimental uncertainties for sin^^efr, M^, "m-t, 
and ruh (the latter assuming ruh ~ 115 GeV). Each column represents the combined results 
of all detectors and channels at a given collider, taking into account correlated systematic 
uncertainties, see Ref. [65] for details. Updated Tevatron numbers can be found in Ref. [66]. 

Another EWPO with a high sensitivity to virtual effects of SUSY particles is the anoma- 
lous magnetic moment of the muon, = {g — 2)^/2. The final result of the Brookhaven 
"Muon g-2 Experiment" (E821) for reads [19] 

aj^^P = (11 659 208 ± 5.8) x 10"^° . (1.52) 

The interpretation of this measurement within SUSY strongly depends on the corresponding 
SM evaluation. The SM prediction depends on the evaluation of the hadronic vacuum po- 
larization and light-by- light contributions. The former have been evaluated by Refs. [67-70], 
the latter by Ref. [71], but there is a recent reevaluation [72], describing a possible shift of 
the central value by 5.6 x 10^^°. Depending on which hadronic evaluation is chosen, the 
difference between experiment and the SM prediction lies between the two values (including 
the updated QED result from Ref. [73]) 

"If"-"?'" ( [68]+ [71]) = (31.7 ±9.5) X 10"^° : 3.3 a , (1.53) 
^exp_^theo ^ ^Qj^^ ^72]) = (20.2 ±9.0) X 10^1° : 2.1a . (1.54) 

These evaluations are all obtained with a Aa^ad determination from e~^e~ data. Recent 
analyses concerning r data indicate that uncertainties due to isospin breaking effects may 
have been underestimated earlier [69]. Furthermore new data obtained at KLOE [74], where 
the radiative return is used to obtain data on Aohad, agrees with the older e~^e~ data. 
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This, together with a continuing discussion about the uncertainties inherent in the isospin 
transformation from r decay, has led to the proposal to leave out the r data in the Acthad 
determination, resulting in the estimate [75] 

- af^° = (25.2 ± 9.2) x 10"^° : 2.7(7 . (1.55) 
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Chapter 2 



Theoretical evaluation of precision 
observables 



2.1 Regular isat ion and renormalization of supersym- 
metric theories 

2.1.1 Basic strategy 

In higher-order perturbation theory the relations between the formal parameters and measur- 
able quantities are different from the tree-level relations in general. Moreover, the procedure 
is obscured by the appearance of divergences in the loop integrations. For a mathematically 
consistent treatment one has to regularize the theory, e.g. by dimensional regularization 
(DREG), where the regularization is performed by analytically continuing the space-time 
dimension from 4 to D [76,77]. But then the relations between the physical quantities and 
the parameters become cut-off-dependent. Hence, the parameters of the basic Lagrangian, 
the "bare" parameters, have no physical meaning. On the other hand, the relations between 
measurable physical quantities, where the parameters drop out, are finite and independent 
of the cut-off. It is therefore in principle possible to perform tests of the theory in terms of 
such relations by eliminating the bare parameters. 

Alternatively, one may replace the bare parameters by renormalized ones by multiplica- 
tive renormalization for each bare parameter qq, 

ao = Zad = a + 5a (2-1) 

with renormalization constants Za different from 1 by a higher-order term. The renormalized 
parameters a are finite and fixed by a set of renormalization conditions. The decomposition 
()2.ip is to a large extent arbitrary. Only the divergent parts are determined directly by the 
structure of the divergences of the loop amplitudes. The finite parts depend on the choice 
of the explicit renormalization conditions. These conditions determine the physical meaning 
of the renormalized parameters. 

Before predictions can be made from the theory, a set of independent parameters has to 
be taken from experiment. In practical calculations the free SM parameters are usually fixed 
by using a, G^, Mz, mj, as (and possibly entries of the quark and lepton mass matrices. 
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if the off-diagonal entries are not neglected) as physical input quantities. They have to be 
supplemented by the empirically unknown input parameters for the Higgs sector and the 
SUSY breaking sector. Differences between various schemes are formally of higher order 
than the one under consideration. The study of the scheme dependence of the perturbative 
results, possibly after improvement by resummation of the leading terms, gives an indication 
of the possible size of missing higher-order contributions. 

On the theoretical side, a thorough control of the quantization and the renormalization 
of the MSSM as a supersymmetric gauge theory, with spontaneously broken gauge symme- 
try and softly broken supersymmetry, is required. This is not only a theoretical question 
for establishing a solid and consistent theoretical framework but also a matter of practical 
importance for concrete higher-order calculations, where the quantum contributions to the 
Green functions have to fulfil the symmetry properties of the underlying theory. An increas- 
ing number of phenomenological applications has been carried out in the Wess-Zumino gauge 
where the number of unphysical degrees of freedom is minimal, but where supersymmetry 
is no longer manifest. 

Moreover, a manifestly supersymmetric and gauge-invariant regularization for divergent 
loop integrals is missing. The prescription of DREG preserves the Lorentz and the gauge 
invariance of the theory, apart from problems related to the treatment of 75 in dimensions 
other than 4. In supersymmetric theories, however, a Z?- dimensional treatment of vector 
fields leads to a mismatch between the fermionic and bosonic degrees of freedom, which 
gives rise to a breaking of the supersymmetric relations. This led to the development of 
dimensional reduction (DRED) [78]. In this scheme only the momenta are treated as D- 
dimensional, while the fields and the Dirac algebra are kept 4-dimensional. It leads to 
ambiguities related to the treatment of 75 [79], and therefore cannot be consistently applied 
at all orders (for a review, see Ref. [80]). Hence, renormalization and the structure of coun- 
terterms have to be adapted by exploiting the basic symmetries expressed in terms of the 
supersymmetric BRS transformations [81]. An additional complication in the conventional 
approach assuming an invariant regularization scheme, however, arises from the modification 
of the symmetry transformations themselves by higher-order terms. 

The method of algebraic renormalization, applied in Ref. [82] to the electroweak SM and 
in Ref. [83] for the MSSM, avoids the difficulties of the conventional approach. The theory is 
defined at the classical as well as the quantum level by the particle content and by the basic 
symmetries. The essential feature of the algebraic method is the combination of all symme- 
tries into the BRS transformations leading to the Slavnov- Taylor (ST) identity. In this way, 
the theory is defined by symmetry requirements that have to be satisfied after renormaliza- 
tion in all orders of perturbation theory. In the case of symmetry violation in the course of 
explicitly calculating vertex functions in a given order, additional non-invariant counterterms 
are uniquely determined to restore the symmetry, besides the invariant counterterms needed 
for absorbing the divergences and for the normalization of fields and parameters. Examples 
are given in Ref. [84,85] for supersymmetric QED and QCD and in Ref. [86] for the SM case. 
Explicit evaluations at the one-loop level in supersymmetric models [84,85,87] have shown 
that DRED yields the correct counter terms. 

In the following we discuss the renormalization of several sectors of the MSSM. We focus 
on the sectors that are needed for the one- and two-loop calculations reviewed below and 
restrict ourselves to the order in perturbation theory required there. These sectors are the 



23 



SM gauge bosons, the electric charge, the quark and scalar quark sector as well as the MSSM 
Higgs boson sector. 

As mentioned above, many MSSM parameters are not closely related to one particular 
physical observable, so that no obvious 'best choice' exists for their renormalization. Ex- 
amples treated below are tan/3 and the mixing angles in the scalar quark sector. Various 
definitions for these parameters already exist in the literature (the situation is similar to the 
case of the weak mixing angle of the electroweak theory, where the use of several different def- 
initions in the literature caused some confusion in the early days of electroweak higher-order 
corrections). We will briefly comment on some of them below. In view of the large number 
of MSSM parameters there is clearly a need to establish some common standards in the 
literature in order to allow for a transparent comparison of different results. Requirements 
that a renormalization scheme for the whole MSSM should fulfll are in particular a coherent 
treatment of all sectors, applicability for both QCD and electroweak corrections, and numer- 
ical stability. Furthermore aspects of gauge (in-)dependence need to be addressed. When 
formulating renormalization prescriptions for the MSSM particular care has to be taken in 
order to respect the underlying symmetry relations of the theory. While in the SM all masses 
of the particles can be flxed by independent renormalization conditions, in a supersymmetric 
theory various relations exist between different masses. Therefore only a subset of the mass 
parameters of the theory can be renormalized independently. The counterterms for the other 
masses are then determined in terms of the independent counterterms. For a discussion of 
these issues, see e.g. Ref. [88]. 



2.1.2 Gauge boson mass renormalization 

We discuss here the renormalization of the gauge-boson masses in the on-shell scheme [89] 
at the one-loop level. Writing the W and Z self-energies as 



(2.2) 



where the scalar functions E ' (g ) are the transverse parts of the self-energies, and deflning 
to correspond to {—i) times the loop diagrams by convention, we have for the one-loop 
propagators {V = W, Z) 



flfj 



-E^(g2 



(2.3) 



where terms proportional to q^q^ (see eq. ()2.2j) ) have been omitted (they are suppressed if 
the propagator is attached to a light external fermion). 

Resumming all self-energy terms yields a geometric progression for the dressed propaga- 
tors: 



1 + 



Ml 



+ 



+ 



g2-M2 +E^(g2) 



(2.4) 
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The locations of the poles in the propagators are shifted by the self-energies. Consequently, 
the masses in the Lagrangian can no longer be interpreted as the physical masses of the 
W and Z bosons once loop corrections are taken into account. The mass renormalization 
relates these "bare masses" to the physical masses Myy, Mz by 

= ty| + 5M|, (2.5) 

with counterterms of one-loop order. The propagators corresponding to this prescription are 
given by 

instead of ()2.4|) . The renormalization conditions which ensure that M^/^z are the physical 
masses fix the mass counterterms to be 

5M^ = ReS'^(M^), 

(5M| = ReE^(M|). (2.7) 



These are the on-shell renormalization conditions. In an MS (or DR) renormalization, on 
the other hand, the counterterms SM^^, 6Mz are defined such that they essentially only 
contain the divergent (in the limit D — ^ 4) contribution. The renormalized mass parameters 
in this case do not directly correspond to the physical masses. They explicitly depend on 
the renormalization scale. 

While the Z-boson mass is commonly used as an input parameter, Mw is normally traded 
as an input parameter for the Fermi constant G^, which is precisely measured in muon decay. 
The prediction for M\y in terms of G^, Mz, a and the parameters of the theory that enter via 
loop corrections can therefore be compared to the experimental value of Mw, constituting a 
sensitive test of the theory (see below). 

Extending the above on-shell definition to higher orders requires to take into account 
that the pole of the propagator of an unstable particle is located in the complex plane 
rather than on the real axis (which is the case for stable particles). A gauge-invariant 
mass parameter is obtained if the mass is defined according to the real part of the complex 
pole. The expansion around the complex pole leads to a Breit-Wigner shape with a fixed 
width. The experimental determination of the gauge-boson masses, on the other hand, uses 
a Breit-Wigner parametrization with running width for historical reasons. This needs to be 
corrected for by a finite shift in Mw and Mz- (For a more detailed discussion, see Ref. [90] 
and references therein.) 

2.1.3 Charge renormalization 

The electroweak charge renormalization is very similar to that in pure QED. In the on-shell 
scheme, the definition of e as the classical charge in the Thomson cross-section 

^Th = (2.8) 



25 



is maintained. Accordingly, the Lagrangian carries the bare charge Cq = e+5e with the charge 
counterterm 5e of one-loop order. The charge counterterm 5e has to absorb the electroweak 
loop contributions to the ee7 vertex in the Thomson limit. This charge renormalization 
condition is simplified by the validity of a generalization of the QED Ward identity which 
implies that those corrections related to the external particles cancel each other. Hence, for 
5e only two universal contributions are left, 

'-^ = \m^-'-^^, n^(o)-^EV) , . (2.9) 

e 2 cw oq^ q^=o 

The first contribution is given by the photon vacuum polarization, 11'^, for real photons, 
= 0. Besides the charged- fermion loops, it contains also bosonic loop diagrams from 
W^W~ virtual states and the corresponding ghosts, as well as from extra charged particles 
in extensions of the SM. The second term contains the mixing between photon and Z boson, 
in general described as a mixing propagator, A'''^, with S'^^ normalized according to 

All loop contributions to T,'^^ vanish at = 0, except the non-Abelian bosonic loops yield 
E'^^(O) 7^ 0. They are the same in the standard model and in supersymmetric exten- 
sions. S'^^(O) completely vanishes in the background-field quantization of the electroweak 
theory [91]. 

The fermion-loop contributions to the photon vacuum polarization in (|2.9p are analogous 
to the electron loop in standard QED and do not depend on the details of the electroweak 
theory. They give rise to a logarithmic dependence on the fermion masses. While for the 
leptonic contributions the known lepton masses can be inserted, perturbative QCD is not 
applicable in this regime, and quark masses are no reasonable input parameters. 

In order to evaluate the contribution of light fermions, i.e. the leptons and the quark 
flavours except the top quark, it is convenient to add and subtract the photon vacuum 
polarization at = M| and to consider the finite quantity (for the top quark and other 
heavy fermions 11''' (0) can be evaluated directly) 

Refl^(M|) = Rem(M|) - n^(0) . (2.11) 

Splitting it into the contribution of the leptons and the five light quarks yields the quantity 

Aa = A«iept + Atthad = - Refll^,{Ml) - ReIi;[,,(M|) , (2.12) 

which represents a QED-induced shift in the electromagnetic fine structure constant 

a a{l + Aa). (2.13) 

The evaluation of the leptonic content of Aa in terms of the known lepton masses yields at 
three-loop order [92] 

Aaiept = 314.97687- 10"^ (2.14) 
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The 5-flavour contribution of the hght quarks to the shift in the fine structure constant can 
be derived from experimental data with the help of a dispersion relation 



— M|Re / ds'— ^ . , 

37r ^ s'{s'-Ml.-te) 



Achad = - ^ M| Re / ds' : ' (2.15) 



where 

\ — c^(e"^e~ ^ 7* ^ hadrons) 
(T(e+e~ — > 7* — > iJriJi~) 

is an experimental input quantity for the low energy range. Recent compilations yield 
the values Aa = 0.02761 ± 0.00036 [93], Act = 0.02769 ± 0.00035 [94], Aa = 0.02755 ± 
0.00023 [68]. 



2.1.4 Renormalization of the quark and scalar quark sector 

The renormalization of the quark sector can differ from the SM case, since the quark masses 
also appear in the scalar quark sector. Therefore the renormalization of the quark and of 
the scalar quark sector cannot be discussed separately. Since the scalar top and bottom 
quarks are most relevant for the evaluation of EWPO, we will focus on their renormalization 
here. Concerning the EWPO calculation reviewed below, only a renormalization in 0{aa) 
is necessary. 



The top and scaleir top sector 

The t/i sector contains four independent parameters: the top-quark mass mj, the stop masses 
mi^ and m^^j ^"^^ either the squark mixing angle 9i or, equivalently, the trilinear coupling At. 
Accordingly, the renormalisation of this sector is performed by introducing four counterterms 
that are determined by four independent renormalisation conditions. 

In an on-shell scheme, the following renormalisation conditions are imposed (the proce- 
dure is equivalent to that of Ref. [95], although there no reference is made to the mixing 
angle) . 

(i) On-shell renormalization of the top-quark mass yields the top mass counterterm, 

Smt = ^mt[ReEtJm?) + ReEt^(m^) + 2ReEt^(m,')] , (2.16) 

with the scalar coefficients of the unrenormalized top-quark self-energy, Et(p), in the 
Lorentz decomposition 

Et(p) = j^u;_EtJp2) +^uj+T.t^{p'')+mtT.t^{p'') . (2.17) 



(ii) On-shell renormalization of the stop masses determines the mass counterterms 

5ml = ReE,~^^(m|) , 5m| = ReE,~,,(m|) , (2.18) 
in terms of the diagonal squark self-energies. 
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(iii) The counterterm for the mixing angle, 6'^, (entering eq. ()1.20|) ) can be fixed in the 
following way, 



Se, = ^i-ii^ ^^^^^ , (2.19) 

2[mi — mi ^ 



involving the non-diagonal squark self-energy. (This is a convenient choice for the treat- 
ment of 0{as) corrections. If electroweak contributions were included, a manifestly 
gauge-independent definition would be more appropriate.) 

In renormalized vertices with squark and Higgs fields, the counterterm of the trilinear cou- 
pling At appears. Having already specified 66^, the At counterterm cannot be defined inde- 
pendently but follows from the relation 

. „^ 2mt( At — a cot 3) 

sin 2ei = — PAIZ ^ 2.20 

mi — mf 

tl 12 



yielding 



5At = — 
nit 



- sin 26i {5m'f_^ — Sm^^ ) + cos 2^j (m?^ 



—sm29f(m'j-m^r)6mt]. (2.21) 

This relation is valid at 0{as) since both fi and tan/3 do not receive one-loop contributions 
from the strong interaction. 

The bottom and scalar bottom sector 

Because of SU(2)-invariance the soft-breaking parameters for the left-handed up- and down- 
type squarks are identical, and thus the squark masses of a given generation are not inde- 
pendent. The stop and sbottom masses are connected via the relation 

cos^ 01 mf^ + sin^ % m?^ = cos^ 6'^ m| + sin^ m| + ml - mf - cos(2/5) , (2.22) 

with the entries of the rotation matrix in eq. p.l6|) . Since the stop masses have already 
been renormalized on-shell, only one of the sbottom mass counterterms can be determined 
independently. Following Ref. [96], the 62 mass is chosen as the pole mass yielding the 
counterterm from an on-shell renormalization condition, i.e. 

Sml = ReEj^Jml) , (2.23) 

whereas the counterterm for m^^ is determined as a combination of other counterterms, 
according to 

^"^t = — (cos^ OiSml^ + sin^ 6'j5m| - sin^ O^^mf^ - sin 26't(m| - mf^)59i 

COS t/^ V 

-f sin 29i{ml^ - ml^)59i - 2mt 6mt + 2mi, 6mb) . (2.24) 
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Accordingly, the numerical value of m^^ does not correspond to the pole mass. The pole 
mass can be obtained from m^^ via a finite shift of 0{as) (see e.g. Ref. [97]). 

There are three more parameters with counterterms to be determined: the 6-quark mass 
nib, the mixing angle 9i, and the trilinear coupling Ab. They are connected via 

su.2e,^ '''^'^i'-^'T^\ (2.25) 
mi —mi 

01 02 

which reads in terms of counterterms 

2 cos 29i 69-, = sin 26-,-^ + , " - sin 29-,^ ^ . (2.26) 

mb mi — mi mi — mi 

01 02 01 02 

Only two of the three counterterms, Smb, 69^, SAb can be treated as independent, which 
offers a variety of choices. 

As discussed in Ref. [96] a convenient choice is the "m(, DR" scheme, whereas a scheme 
analogous to the one in the t/i sector, involving a bottom pole mass, can lead to artificially 
enhanced higher-order corrections. 

Concerning the renormalization of the top and the bottom mass, there are important 
differences. The top-quark pole mass can be directly extracted from experiment and, due to 
its large numerical value as compared to other quark masses and the fact that the present 
experimental error is much larger than the QCD scale, it can be used as input for theory 
predictions in a well-defined way. For the mass of the bottom quark, on the other hand, 
problems related to non-perturbative effects are much more severe. Therefore the parameter 
extracted from the comparison of theory and experiment [3] is not the bottom pole mass. 
Usually the value of the bottom mass is given in the MS renormalisation scheme, with the 
renormalisation scale fj,^^ chosen as the bottom-quark mass, i.e. m^^(m^^) [3]. 

Another important difference to the top/stop sector is the replacement of cot (3 — > tan /3. 
As a consequence, very large effects can occur in this scheme for large values of /i and 
tan/3 [98]. 

Potential problems with the bottom pole mass can be avoided by adopting a renormal- 
isation scheme with a running bottom-quark mass. In the context of the MSSM it seems 
appropriate to use the DR scheme [78] and to include the SUSY contributions at 0{as) into 
the running. We denote this running bottom mass as m^^''^^^^(/i°^). 

The "m^ DR" scheme mentioned above uses a DR renormalization for both and Ab. 
In the DR scheme, the 6-quark mass counterterm can be determined by the expression 

5mb = ^mb[ReJ:t^(ml) + ReJ:t^(ml) + 2ReE^^(m^)] , (2.27) 

where J]'^" means replacing the one- and two-point integrals A and Bq in the quark self- 
energies by their divergent parts. The counterterm for the trilinear coupling Ab in the DR 
scheme reads [96] 



SAb^ — 

mb 



tanesReE£(m?J + ^(ReEg;(m?J + ReEg;(m?J) 
+ tan 9i (cos' ^,-Re E^;;(m? ) + sin^ ^,-Re Egj(m|) 
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- (Re E^^(m?) + Re Ef^{m',) + 2Re E^^(m^) 



+ - (2 tan 6'f,m;, 



1 



[mr — mr 



jsin2^^) 



X 



2V 2mb'"~^^ ^2 

Re (m^) + Re ^^-(m^) + 2Re 11^7(^1 



(2.28) 



The counterterms for the mixing angle, 66^, and the bi mass, (5m? , are dependent quantities 
and can be determined as combinations of the independent counterterms, invoking (j2.24j) 
and (jT^ . 



59: 



1 1 
mi — mi 

01 02 



mi,6Ab + tan d^Sm^^ + — (m?^ — m?_^) sin 26^^ — 2 tan 

- tan 6'^ (cos^ 0i5m\ + sin^ 6'f:5m| - sin 26lt-(m| - mfJdOi - 2mt5mt^ , (2.29) 

5m?^ = tan^ 9i5m^^ + 2tan^^m{,5Afc + 2^^ — sin^ 6i{m^^ — fn^^ + (1 — tan^ 9i)mhjSmi, 
+ (1 — tan^ 9^) (^cos^ 9^6m^^ + sin^ 9^6mf^ — sin 29f{mj^ — m^^)59i — 2mt6m 



[2.30) 



DR 



The renormahzed quantities in this scheme depend on the DR renormahzation scale fi 

In order to determine the value of mj^^'^^^^(/i^^) from the value m^^(/i^^) that is 
extracted from the experimental data one has to note that by definition 77^^^'^^^^ contains 
all MSSM contributions at 0{as), while m^^ contains only the 0{as) SM correction, i.e. the 
gluon-exchange contribution. Furthermore, a finite shift arises from the transition between 
the MS and the DR scheme. 



The expression for m^^'^^^^(/i°^) is most easily derived by formally relating 



m. 



DR.MSSM 



to the bottom pole mass first and then expressing the bottom pole mass in terms of the MS 
mass (the large non-perturbative contributions affecting the bottom pole mass drop out in 



the relation of ^j^^'^^^^ iq m^^). Using the equality m^^ + (5mP^ = 777,^r,mssm _|_ ^^dr,mssm 
and the expressions for the on-shell counterterm and the DR counterterm one finds 



m 



DR.MSSM/, DR 



m 



OS 



b vr- / -^0 ■ 2 

fin 



+ -m,(E^°(m,2) + S^°(m,2)) + m, S^^"(m^) . (2.31) 



Here the S are the UV-finite parts of the bottom quark self-energy coefficients. They 
depend on the DR scale ijP^ and are evaluated for on-shell momenta, = m^. Inserting 



m 



OS 



m^^{Mz)V^'^\ where 



.shift 



(A [m 



MS\2 
b ) 



(2.32) 
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one finds the desired expression for mf^, 

^DR,MSSM^^DR) ^ ^fs zW"^'^'' + Im^^E^.M + J^^M"")) + ^^^tsi^D ■ (2-33) 

2.1.5 MSSM Higgs boson sector renormalization 

In order to perform liigfier-order calculations in the Higgs boson sector, the renormalized 
Higgs boson self-energies are needed (see Sect. I2.7|l . The parameters appearing in the Higgs 
potential, see eq. (jl.lj) . are renormalized as follows: 

Ml^Ml + 6Ml Th-^n + Sn, (2.34) 

M^igg, M^igg, + 5M^igg, , tan /5 tan /3 ( 1 + 5tan/? ) . 

mjj± — > mjj± + 5m%± 

The renormalization of My/ and Mz has been described in Sect. 12.1.21 M^jggs denotes the 
tree-level Higgs boson mass matrix given in eq. ()1.9j) . Th and Th are the tree-level tadpoles, 
i.e. the terms linear in h and H in the Higgs potential. 

The field renormalization matrices of both Higgs multiplets can be written symmetrically. 




2 



6ZhH 1 + \^Zhh 




(2.35) 



and for the charged Higgs boson 

H^(l + 5Zh-h+) ■ (2.36) 
For the mass counterterm matrices we use the definitions 

The renormalized self-energies, £(p^), can now be expressed through the unrenormalized 
self-energies, S(j9^), the field renormalization constants and the mass counterterms. This 
reads for the CP-even part, 

thhip") = ^hhip') + SZhhip" - ml,^J - 5ml (2.38a) 

^hHip^) = ^huip^) + 5ZhH{V^ - |("^h,tree + "^/f.tree)) " ^'^hH^ (2.38b) 
^Hh{P^) = ^Hh{P^) + 5Zhh{P^ - "^H,tree) " ^^^'h ) (2.38c) 

and for the charged Higgs boson 

t.H-H+{p^) = T.H-H+{p^) + SZu-H+iv^ - m]j±) - 5m]j± . (2.39) 

Inserting the renormalization transformation into the Higgs mass terms leads to ex- 
pressions for their counter terms which consequently depend on the other counter terms 
introduced in ()2.34|) . 
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For the CP-even part of the Higgs sectors, these counter terms are: 

6ml = 6 Ml cos2(a - f3) + 6 Ml sm^{a + (3) (2.40a) 

+ 2MZSWCW ^^'^^ ~ ^) sin^(a - P) + STh sin(a - + cos^(a - /?))) 
+ StanP sin /3 cos /3 {Ml sin 2(a - /?) + M| sin 2(a + /?)), 

^m^^ = i(5Mlsin2(a-/5) -5M|sin2(a + /?)) (2.40b) 

- Stanp sin /? cos (3 {Ml cos 2 (a - /?) + M| cos 2(a + /?)), 

= 5Ml sin2(a - l3) + 6Ml cos^{a + (3) (2.40c) 
~ 2MZSWCW ^^^H - + sin2(a - /?)) + sin(a - /3) cos^(a - /?)) 

— (5tan/3 sin (3 cos /5 (M^ sin 2 (a — /5) + M|sin2(a + /5)) . 

For the charged Higgs boson it reads 

5m]j± = 6Ml + 5Ml, . (2.41) 

For the field renormahzation it is sufficient to give each Higgs doublet one renormalization 

constant, 

7^1 ^ (1 + \5Zn,)ni, 7^2 ^ (1 + \5Zn,)n2 • (2.42) 

This leads to the following expressions for the various field renormalization constants in 
eq. (ESSD: 

5Zhh = sin^a SZt-i-^ + cos^a 5.^7^2 5 (2.43a) 
SZhH = sin a cos a {SZ-^^ — SZ-i^J, (2.43b) 
SZhh = cos^aSZ-Hi + sin^a 6Z'H2, (2.43c) 
5Zh-h+ = sin^P SZn, + cos^P SZn^ . (2.43d) 

The counter term for tan P can be expressed in terms of the vaccuum expectation values as 

X s s 

6tanP = - {6Zn, - SZnJ + — - — , (2.44) 

2 V2 Vi 

where the Svi are the renormalization constants of the vf. 

v,^{l + 6ZnJ {vi + 6v,) , V2^{1 + 6Zn,) {v2 + 6V2) . (2.45) 

The renormalization conditions are fixed by an appropriate renormalization scheme. For 
the mass counterterms besides the on-shell conditions for M]y and Mz (see eq. (|2.7p ) also 
Ma can be renormalized on-shell: 

6Ml = ReJ:AA{Ml). (2.46) 
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Since the tadpole coefficients are chosen to vanish in all orders, their counter terms follow 
from T{h,H} + ST{h,H} = 0: 

6n = -Th, 6Th = -Th . (2.47) 

For the remaining renormalization constants for 6 tan f3, dZ-j-i^ and 5Z-}i^ several choices are 
possible, see e.g. Ref. [99,100]. A convenient choice is a DR renormalization of (5tan/5, 5Z-yi^ 
and 5Z-}-i^, 

5tan/? = 5tan/5^ = -j^^ [Re ^^.(m^,,,,,,) - Re SW«,tree)]''\ (2.48a) 
5Zn, =5Z^ = - [Re S'^^ |„=o] , (2.48b) 
5Zn, = 5Z^ = - [Re E',, .^^J . (2.48c) 



2.2 Sources of large SUSY corrections 
2.2.1 Possible sources 

Besides the known sources of sizable higher-order corrections in the SM, e.g. contributions 
enhanced by powers of mt or logarithms of light fermions, there are additional sources of 
possibly large corrections within the MSSM: 

• Large corrections can arise not only from loops containing the top quark, but also its 
scalar superpartners. In the MSSM Higgs sector, Yukawa corrections from the top and 
scalar top quark sector can be especially large. The one-loop corrections, for instance 
to the upper bound on the mass of the lightest CP-even Higgs boson, can reach the 
level of 100%. The leading one- loop term from the top and scalar top sector entering 
the predictions in the Higgs sector is given by [101] 

^G.mtlogl'^) . (2.49) 



m 



While the Higgs sector of the MSSM is CP-conserving at tree level, large CP-violating 
effects can be induced by the loop corrections. 

Effects from the h/h sector of the MSSM can also be very important for large values of 
tan/3 and /i. 

The h Yukawa coupling can receive large SUSY corrections, yielding a shift in the 
relation between the b quark mass and the corresponding Yukawa coupling [98], 

Vb = TTTT^ • 2.50 

V cos p 1 + Arrih 

The quantity Arrih contains in particular a contribution involving a gluino in the loop, 
which gives rise to a correction proportional to (os/im^ tan/?), which can be large. 
For Anif, ^ — 1 the b Yukawa coupling even becomes non-perturbative. This issue is 
discussed in Sect. 12.2.21 
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• Besides the scalar quark sector, SUSY theories have further possible sources of large 
isospin sphtting, which can give large contributions to the p parameter [102,103]. 

• Soft SUSY-breaking masses can induce splittings in the supersymmetric coupling rela- 
tions [104,105] (i.e. the equality of a SM coupling gi with the corresponding supersym- 
metric coupling hi). If scalar superpartners have masses at a high scale M, and all the 
other masses are light with mass m ~ M^eak, the resulting corrections are given by 

^-l-^A6aog^, (2.51) 
gi[m) 16 TT^ m 

where A6j is the one-loop beta function coefficient contribution from all light particles 
whose superpartners are heavy. If M ^ m these corrections to the SUSY coupling 
relation can be sizeable. 

• Another type of possibly large corrections in supersymmetric theories are the so-called 
Sudakov logs (see Ref. [106] and references therein). They appear in the form of 
log(g^/M|ugY) (where q is the momentum transfer) in the production cross sections of 
SUSY particles at e+e~ colliders. 

• In general, SUSY loop contributions can become large if some of the SUSY particles 
are relatively light. 

2.2.2 Resummation in the 6/6 sector 

The relation between the bottom-quark mass and the Yukawa coupling yb, which in lowest 
order reads nib = ybVi/\/^, receives radiative corrections proportional to ybV2 = ybtan/Svi. 
Thus, large tan /3-enhanced contributions can occur, which need to be properly taken into 
account. As shown in Refs. [98,107] the leading terms of C(a6(as tan/3)") can be resummed 
by using an appropriate effective bottom Yukawa coupling. 

Accordingly, an effective bottom-quark mass is obtained by extracting the UV-finite 
tan /3-enhanced term Am^ from eq. (j2.33p (which enters through Sf,.,) and writing it as 
1/(1 + Am;,) into the denominator. In this way the leading powers of (a^ tan /5)" are correctly 
resummed [98,107]. This yields 

DRMSSM m rnf^iMz)b^''''' + |m,(Efi-(m,2) + ^IiM) + mb E^-(m2) 
mb ' (/. ) = ^^^^ , (2.52) 

where = + Arrib denotes the non-enhanced remainder of the scalar 6-quark self- 
energy at 0{as), and b^^^^^ is given in ()2.32j) . The tan /3-enhanced scalar part of the 6-quark 
self-energy, Am^, is given at 0{as) hj^ 

2 

Arrib = — astBxi(3 prrig I{m\ ,m? ,m?), (2.53) 
37r ^ ^ ^ 



^There are also corrections of 0{at) to Amh that can be resummed [107]. These effects usually amount 
up to 5-10% of the 0{as) corrections. Since in this report we only consider 0{abas) contributions, these 
corrections have been omitted. Further corrections from subleading resummation terms can be found in 
Ref. [108]. 
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with 



Tf 2 2 2\ 




02 02 01 ' 01 tf 01 ' y y bo ^ y 02' 




(2.54) 



and Amf, > for > 0. 

In the "m;, DR" defined above, the effective bottom-quark mass as given in eq. (j2.52j) 
should be used everywhere instead of the DR bottom quark mass. This also applies to the 
bottom mass in the sbottom-mass matrix squared, eq. ()1.14|) . from which the sbottom mass 
eigenvalues are determined. The effects of Am?,, i.e. the leading effects of C(as), can be 
incorporated into a lowest-order result (e.g. the one-loop results for the renormalized Higgs 
boson self-energies, see Sect. 12. 7|) by using the effective bottom-quark mass of eq. ()2.52|1 (or 
the correspondingly shifted value in other renormalization schemes). 

2.3 Electroweak precision observables in the MSSM 

In this section we briefly introduce the electroweak precision observables that are discussed 
in this report. A description of the current status of their theoretical evaluation within the 
MSSM will be given in the following sections and the remaining theoretical uncertainties will 
be discussed. 

The current experimental status of the EWPO and prospective improvements of their 
precision in the future have been summarized in Sect. 11.3.41 In order to fully exploit the 
experimental precision of the EWPO the theoretical uncertainties should be reduced signif- 
icantly below the level of the experimental errors. 

Concerning the theoretical predictions, two kinds of uncertainties need to be taken into 
account: the theoretical uncertainties from unknown higher-order corrections ("intrinsic" 
theoretical uncertainties) and the uncertainties induced by the experimental errors of the 
input parameters ("parametric" theoretical uncertainties). The parametric uncertainty in- 
duced by the known input parameters (in the SM case in particular mt and Acthad) needs to 
be reduced in order to increase the sensitivity to the unknown parameters of the model (in 
the SM case Mh). 

The EWPO discussed in the following sections are: 

• The W boson mass can be evaluated from 



where a is the flne structure constant and the Fermi constant. This relation arises 
from comparing the prediction for muon decay with the experimentally precisely known 
Fermi constant. The radiative corrections are summarized in the quantity Ar, derived 
first for the SM in Ref. [109]. The prediction for Mw within the SM or the MSSM 
is obtained from evaluating Ar in these models and solving eq. (j2.55p in an iterative 
way. The theory status of the prediction for Mw is reviewed in Sect. 12.51 




(2.55) 
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Another important group of EWPO are the Z boson observables, among which we 
mostly concentrate on the effective leptonic weak mixing angle at the Z boson reso- 
nance, sin^ ^eff- It can be defined through the form factors at the Z boson pole of the 
vertx coupling of the Z to leptons (/). If this vertex is written as il'y'^{gy — gx'^^)lZ^ 
then 

2n 1/^1 T.^9\ 



sin' = - 1 - Re^ . (2.56) 

At the tree level this amounts to the sine of the weak mixing angle, sin^ O^j = 1 — 
My^/M"^, in the on-shell scheme. Loop corrections enter through the form factors Qy 
and ^^A- The theoretical evaluation is reviewed in Sect. 12.61 

The quantity Ap, 

E^_E-(0) 
^ " M| ' ^^-^^^ 

parameterizes the leading universal corrections to the electroweak precision observables 
induced by the mass splitting between fields in an isospin doublet [102]. S^'^(O) 
denote the transverse parts of the unrenormalized Z and W boson self-energies at 
zero momentum transfer, respectively. The induced shifts in the two above described 
observables are given in leading order by 



5Mv,^^^^Ap, ^sin^^eff^-^^Ap. (2.58) 

2 Cy^ S^r C^r S^r 



The theoretical evaluation of Ap is discussed in Sect. 12.41 

Another very powerful observable for constraining the parameter space of the MSSM 
is the mass of the lightest CP-even Higgs boson, rrih. If the Higgs boson will be found 
at the next generation of colliders, its mass will be measured with high precision. We 
therefore refer to rrih also as an EWPO. While rrih is bounded from above at tree- level 
by < Mz, it receives large radiative corrections. The leading one- loop contribution, 
arising from the t/t sector, reads [101] 

= ^ . . < log { . (2.59) 



v^7r2 sin2/5 



ml 



The loop corrections, entering via Higgs-boson propagator corrections, can shift vrih by 
50-100%. The theoretical status is reviewed in Sect. 12.71 

As a further precision observable that we investigate in detail in this report we consider 
the anomalous magnetic moment of the muon, = {g — 2)^. It is related to the 
photon-muon vertex function F^^^p as follows: 

u{p')T ^,fiAp{p, -p\ q)u{p) = u{p') [7pFy(g2) + {p + p')pFM{q^) + ■ ■ ■] u{p) , 

= -2m^FM(0) , (2.60) 

where F/v/(g^) = at tree-level. Non-zero values are induced via loop corrections. The 
theoretical evaluation is discussed in Sect. 
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2.4 The p parameter 



We start our discussion with the quantity Ap, see eq. ()2.57|) . which parametrizes in particular 
the leading contributions from loops of scalar quarks and leptons to the W-hoson mass and 
the Z-boson observables. 



2.4.1 One-loop results 

In the SM the dominant contribution to Ap at the one-loop level arises from the t/h doublet 
due to its large mass splitting. With both fermion masses non-zero it reads 

Apo'^ = ^^FoK,m^), (2.61) 

with 

Fo(x,i/) =x + 2/-^^log-. (2.62) 
x-y y 

Fo has the properties Fo(m^,m^) = Fo(m^,m^), Fo(m^,m^) = 0, Fo(m^,0) = m^. Therefore 
for rrit ^ nib eq. ()2.61|) reduces to the well known quadratic correction 

ApSM = rnl (2.63) 

Within the MSSM the dominant SUSY correction at the one-loop level arises from the 
scalar top and bottom contribution to eq. ()2.57|) . see Fig. 12.11 



V / V 
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Figure 2.1: Feynman diagrams for the contribution of scalar quark loops to the gauge boson 
self-energies at one-loop order. 



For rrtfc 7^ it is given by 
3G„ 



Ap; 



SUSY 

'0 



8^2 



- sin^ 6'i cos^ 614^0 (m|, m|) - sin^ 6^ cos^ 6'^Fo(m|^, m?J 
+ cos^ 9^ cos^ OiFq (m| , m?^ ) + cos^ 9f sin^ Oj^Fq (m|^ , m?^ ) 
+ sin^ Of cos^ e^Fo (m| , ) + sin^ 9^ sin^ O-^Fq (m| , m- 



(2.64) 



The size of the SUSY one-loop contributions are shown for an exemplary case in Fig. 12.21 as 
a function of MgusY- The parameter Mgusy is defined by setting the soft SUSY-breaking 
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parameters in the diagonal entries of the stop and sbottom mass matrices equal to each other 
for simplicity, 

MsusY = Mq = M^ = M^, (2.65) 
see eq. fll.l4|l . We furthermore use the shorthands 

Xt = At — n I tan [3^ Xfo = — /itan/? . (2.66) 

The other parameters in Fig. 12.21 are tan /3 = 3 and Xt = 0, 2Msusy- In this case can 
reach values of up to 2 x 10~^. The line for Xt = 2Msusy starts only at Msusy ~ 300 GeV. 
For lower values of Msusy one of the scalar top mass squares is below zero. 

0.0020 I — , — , — , — , — I — , — , — , — , — I — , — , — , — , — I — , — , — , — , — I 




Figure 2.2: One-loop contribution of the (t, b) doublet to Ap as a function of the common 
squark mass Msusy for tan /? = 3, and X5 = and = or 2Msusy- 



2.4.2 Results beyond the one-loop level 

SM results 

Within the SM the one-loop 0{a) result from the contribution of the t/b doublet has been 
extended in several ways. The dominant two-loop corrections arise at 0{aas) and are given 
by [110] 

^pSM,aa. = -ApSM^^ (l + ^2/3) _ (2.67) 
O TT 

These corrections screen the one-loop result by approximately 10%. Also the three-loop 
result at 0{aa1) is known. Numerically it reads [111] 



Apr' ^ = < — ■ 14.594... . (2.68) 
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Furthermore the leading electroweak two-loop contributions of (9(G^m^) have been calcu- 
lated. First the result in the approximation = had been evaluated [112], 



SM,G2 



^Pl\M„=0 



SM 
|AfH=0 



1287r 

19 - 277^ 



4 

4 



SM 
\M„=0 



fSM 



(2.69) 

extends 



Later the full 0{Gji^mf) result for arbitrary Mh became available [113], where 5i|a/j:^=o 
to 

^iIMh^o = 19 - 2n^ + fct(m„ Mh). (2.70) 

The leading two-loop contribution to Ap in an asymptotic expansion for large of 
0(G^M|fM^) was obtained in Ref. [114]. It turned out to be numerically small. 

Leading electroweak three-loop results of 0{G^^mf) and (9(G^asm^) became available 
more recently [115,116]. Numerically they read in the case Mh = 0: 

G,, 



Ap. 



SM,G^ 
2\M„=0 



Ap[ 



SM.G^a^ 
'2\Mh=0 



8^2 

Gn 



• 249.74 , 



8^2 



— ) ■ 2.9394 . 

7r 



(2.71) 
(2.72) 



For the case Mh 7^ the result has been obtained in several limits, allowing a smooth inter- 
polation, see Ref. [116] for details. Most recently also the leading O(G^M^M^) contribution 
was obtained [117]. Besides for very large values of Mh it is numerically insignificant. 

The SUSY corrections at C»(aaJ 

The leading two-loop corrections arising in the MSSM (beyond the SM part) have been 
evaluated at 0{aas) [97] and 0{a1, atab, a^) [118,119] (the latter in the limit of large Mgusv)- 
The leading 0{aas) corrections to the scalar quark loops consist of the diagrams shown in 
Fig. 12.31 (supplemented with the corresponding diagrams for the subloop renormalization, 
see Ref. [97]). The diagrams can be divided into three groups: the pure scalar contribution 
(diagrams a-c), the gluonic correction (diagrams d-j, where the gluon-loop contribution, 
diagrams i-j, is zero) and the gluino exchange correction (diagrams k-n). 

The pure scalar quark diagrams give a vanishing contribution. The gluonic correction 
can be cast into a compact formula [97]: 



A SUSY 
r Ijgluon 
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4V2 7r2 71 



sin^ 9); cos^ 9fFi (m: 



sin^ 6r cos^ 9rFi{m 



+ cos^ Oi cos^ e^Fi (m| , m?^ ) cos^ 6^ sin^ O^Fi 
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+ sin^ 9t cos^ e-^Fi (m| , J sin^ 9^ sin^ (m| , m| 



^2.73) 



with 



Fi{x,y) = x + y-2 log - 
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Figure 2.3: Feynman diagrams for the contribution of scalar quark loops to the gauge-boson 
self-energies at two-loop order. 



+7 log' --2{x- y)U2 1 - - , 2.74 

[x-yr y \ yJ 

where Fi has the properties Fi(m^,m^) = Fi(m^,m^), Fi(m^,m^) = 0, Fi(m^,0) = m^(l + 
7r^/3). The gluino exchange correction results in a lengthy formula, see Ref. [97], and is not 
given here. It decouples for rrig — >• oo. 

The analytical formula for the 0{aas) corrections given in eq. (j2.73|) is expressed in terms 
of the physical squark masses, i.e. an on-shell renormalization has been carried out for all 
four squark masses. As discussed in Sect. 12.1.41 SU(2) invariance leads to a relation between 
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the stop and sbottom masses, so that not all four masses can be renormalized independently. 
This results in a finite mass shift of 0{as) that is given, if expressed in terms of m^^^, as 
the difference between the counterterm of eq. ()2.24|1 and the on-shell counterterm. If the 
two-loop result is expressed in terms of the on-shell masses, this mass shift appears in the 
relation between the physical squark masses and the (unphysical) soft SUSY-breaking mass 
parameters in the squark mass matrices, see eq. ()1.14j) . While this shift is formally of higher 
order in the evaluation of the masses that are inserted in the two-loop result, it needs to be 
taken into account in the one-loop result. This gives rise to an extra contribution compared 
to the results discussed in Sect. 12.4.11 see Ref. [97] for a more detailed discussion. 

The SUSY corrections at 0{af), 0{atab), 0{al) 

Furthermore the leading O^a^), 0{atab), 0{al) corrections to Ap have been evaluated 
in the limit Msusy oo [118,119]. The rrit dependence of Ap differs between the pure 
SM contribution and the additional SUSY corrections. Within the SM, the corrections are 
~ mf for the one-loop and ~ mf for the two-loop correction, leading to sizable shifts in the 
precision observables. The additional SUSY corrections at the one-loop level (from scalar 
quark loops), on the other hand, do not contain a prefactor ~ m^. In the electroweak two- 
loop corrections it is no longer possible to separate out the pure SM contribution because 
of the extended Higgs sector of the MSSM. The leading electroweak two-loop corrections in 
the MSSM are therefore of 0{G'^mf) (as in the SM case) and potentially sizable. 

The leading contributions of Oi^a^), 0{atab) and 0{al) have been derived by extracting 
the contributions proportional to yf, ytUb and yl, where 

y/2mt V^rub . . 

Vt = — ^-T, Vb = -T ■ (2.75) 

V smp V cos p 

The coefficients of these terms could then be evaluated in the gauge-less limit, i.e. for 
Mw, Mz ^0 (keeping cw = Mw/Mz fixed). 

For the Higgs masses appearing in the two-loop diagrams the following relations have 
been used, arising from the gauge-less limit 

m^± = Ml , ml = , ml± = . (2.76) 

Applying the corresponding limit also in the neutral CP-even Higgs sector would yield for the 
lightest CP-even Higgs-boson mass ml = (and furthermore rnjj = M\, sin a = — cos (3 , 
cos a = sin/3 ). Since within the SM the limit turned out to be only a poor 

approximation of the result for arbitrary M^, ml has been kept non-zero (which formally 
is a higher-order effect). Keeping m^ as a free parameter is also relevant in view of the 
fact that the lightest MSSM Higgs boson receives large higher order corrections, which shift 
its upper bound up to 135 GeV (for Msusy < 1 TeV and mt = 175 GeV), see Sect. 12. 7( 
These corrections can easily be taken into account in this way (in the Higgs contributions at 
one-loop order, however, the tree-level value of mh should be used). Keeping a arbitrary is 
necessary in order to incorporate non SM-like couplings of the lightest CP-even Higgs boson 
to fermions and gauge bosons. 

On the other hand, keeping all Higgs-sector parameters completely arbitrary is not pos- 
sible, as the underlying symmetry of the MSSM Lagrangian has to be exploited in order 
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to ensure the UV-finiteness of the two-loop corrections to Ap. Thus only those symmetry 
relations have been enforced in the neutral CP-even Higgs sector which are explicitly needed 
in order to obtain a complete cancellation of the UV-divergences. 

It is convenient to discuss the 0{a^ oc G^^mf) SUSY contributions to Ap separately, i.e. 
the case where i/b = 0. The Ola'f) corrections are by far the dominant subset within the 
SM, i.e. the 0{atab) and C(a^) corrections can safely be neglected within the SM. The same 
is true within the MSSM for not too large values of tan/5. It is well known [120] that the 
SUSY sector of the MSSM decouples if the general soft SUSY-breaking scale goes to infinity 
(corresponding to Msusy ^ cxd in the one-loop result given above). The leading contributions 
of (9((j^mj) in the case where the scalar quarks are heavy is therefore obtained in the limit 
where only the two Higgs doublet sector of the MSSM is active [118,119], corresponding to 
the limit Msusy — ^ oo. 

In Ref. [118] the result has been obtained in the simplified case with tree-level Higgs 
boson masses. In the limit Mw,Mz the neutral CP-even Higgs boson masses at the 
tree-level reduce to 



0, 



m 



H 



Mi 



[2.77) 



In this limit also the relation between the angles a and (3, see eq. fll.6|) . becomes very 
simple, a = j3 — 7r/2, i.e. sin a = — cos/5 , cos a = sin/5 . The only remaining scales left 
are the top quark mass, m^, the CP-odd Higgs boson mass. Ma, and tan/5 (or sin/? = 

In the limit of large tan/? (i.e. (1 — sin^/5 



tan/5/A/l + tan2/5 
particularly simple form. One obtains 



<^ 1) the result takes a 
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(2.78) 



Thus, for large tan (3 the SM limit with MfJ^ (see eq. (|2.69j) ) is reached. 
Keeping tan /5 arbitrary but expanding for large values of Ma yields 
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where A = m1/M\. In the limit A Q one obtains 

Q2 
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(2.80) 
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i.e. exactly the SM limit for — > is reached. This constitutes an important consistency 
check: in the limit A ^ the heavy Higgs bosons are decoupled from the theory. Thus only 
the lightest CP-even Higgs boson should remain, which has in the (9(G'^m^) approximation 
(neglecting higher-order corrections) the mass m/j = 0, see eq. ()2.77|) . As already observed 
in Ref. [97], the decoupling of the non-SM contributions in the limit where the new scale 
(i.e. in the present case Ma) is made large is explicitly seen here at the two-loop level. 

Now we turn to the full 0{af) corrections. As discussed in Ref. [119], an UV-finite result 
could only be obtained if the relations in eq. ()2.7(i|l are taken into account. The masses of the 
neutral Higgs bosons as well as the mixing angle could be kept as 'independent' parameters, 
i.e. they can be obtained taking into account higher order corrections. The full result without 
the tree-level relations is rather lengthy and can be found in Ref. [119]. 

Now also the 0{atah), 0{al) SUSY corrections are considered. The structure of the 
fermion doublet requires that further symmetry relations are taken into account. Within the 
Higgs boson sector it is necessary, besides using eq. ()2.7(i|l . also to use the relations for the 
heavy CP-even Higgs boson mass and the Higgs mixing angle, 

= , sin a = — cos/5 , cos a = sin/3 . (2.81) 

On the other hand, rrih can be kept as a free parameter. The couplings of the lightest 
CP-even Higgs boson to gauge bosons and SM fermions, however, become SM-like, once 
the mixing angle relations, eq. ()2.8H1 . are used. Furthermore, the Yukawa couplings can no 
longer be treated as free parameters, i.e. eq. ()2.75|1 has to be employed, which ensures that 
the Higgs mechanism governs the Yukawa couplings. Corrections enhanced by tan/5 thus 
arise only from the heavy Higgs bosons, while the contribution from the lightest CP-even 
Higgs boson resembles the SM one. 



2.4.3 Results in the NMFV MSSM 

The existing corrections to Ap within the NMFV MSSM [25] consist of squark contributions 
based on the general 4x4 mass matrix for both the t/c and the b/s sector, see Sect. I1.2.(j[ 
These corrections are visualized by the Feynman diagrams in Fig. 12.41 They are denoted 
as Ap«. 

The squark contribution Ap^ can be decomposed according to 

Ap^' = Ez + Qz + ^w + Qw, (2.82) 

where S and correspond to different diagram topologies, i.e. to diagrams with trilinear 
and quartic couplings, respectively (see Fig. 12. 4p . The explicit expressions read as follows, 



^ a,b,c4 a, 13 

= -^^^EE{(^r)%(<) + (^r)'^o(<)}, 

^ a a 

Sg' 



43 




Figure 2.4: Feynman diagrams for the squark contributions to the gauge boson self-energies. 
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E {i^i(l)^R7fMrnlj + ^,{l)\RT?A,{mlj] .(2.83) 



Here the indices run from 1 to 2 for Latin letters, and from 1 to 4 for Greek letters. The 
expressions contain the one-point integral Aq and the two-point integral Bqq in B^y{k) = 
dfivBoQ + k^kyBii in the convention of Ref. [121]. The remaining constants and are 
defined as follows, 
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The CKM matrix only affects Sty. Corrections from the first-generation squarks are 
negligible due to their very small mass splitting. Non-minimal flavor mixing of the first 
generation with the other ones has been set to zero, but conventional CKM mixing is basically 
present. Although it is required for a UV finite result, it yields only negligibly small effects. 
Therefore, for simplification, we drop the first generation and restore the cancellation of UV 
divergences by a unitary 2x2 matrix replacing the {23}-submatrix of the CKM matrix. 



CKM 
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cose sme 
— sin e cos e 



(2.85) 



with |e| ~ 0.04 close to the experimental entries [3] of the conventional CKM matrix. 

Since I\p^ is a finite quantity, and the CKM matrix effects (and therefore, the e depen- 
dence) only appear in Siy, it has been shown [25] that Sv(/ (and thus Ap) is symmetric under 
the simultaneous reversal of signs e — e, A —A (see eq. ()1.50p ). i.e. only the relative 
sign has a physical consequence, affecting the results for Ap significantly. In physical terms, 
non-minimal squark mixing can either strengthen or partially compensate the CKM mixing. 
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2.5 Evaluation of Mw 



One of the most important quantities for testing the SM or its extensions is the relation 
between the massive gauge boson masses, Mw and M^, in terms of the Fermi constant, G^, 
and the fine structure constant, a. This relation can be derived from muon decay, where the 
Fermi constant enters the muon lifetime, r^, via the expression 



T,. 



1927r3 



ml 



1 + 1 



m. 



(2.86) 



with F{x) = 1 — 8a; — 12x^ Inx + 8a;^ — x'^. By convention, this defining equation is supple- 
mented with the QED corrections within the Fermi Model, Ag. Results for Ag have been 
available for a long time at the one-loop [122] and, more recently, at the two-loop level [123] 
(the error in the two- loop term is from the hadronic uncertainty). 



Ag = 1.810 — + (6.701 ± 0.002) ( — 

47r V47r 



^2.87) 



Commonly, tree-level W propagator effects giving rise to the (numerically insignificant) term 
3m^/(5M^) in eq. ()2.86|) are also included in the definition of G^, although they are not 
part of the Fermi Model prediction. With the second order term of eq. ()2.87p the defining 
equation for in terms of the experimental muon lifetime, eq. ()2.86|) . yields the value of 



G. 



given m 
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Within a given model, G^ can be calculated in terms of the model parameters. The 
Fermi constant is given by the expression 



G. 



1 + 



E^(0) 



(VB) 



(2.88) 



This equation contains the bare parameters with the bare mixing angle. The term [VB) 
schematically summarizes the vertex corrections and box diagrams in the decay amplitude. 
A set of infrared-divergent "QED correction" graphs has been removed from this class of 
diagrams. These left-out diagrams, together with the real bremsstrahlung contributions, 
reproduce the QED correction factor of the Fermi model result in eqs. (j2.86p . (j2.87p and 
therefore have no influence on the relation between G^ and the model parameters. 

Equation ()2.88p contains the bare parameters eo, M{y, s^r. Expanding the bare parame- 
ters and keeping only terms of one-loop order yields the expression. 
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(2.89) 



which is equivalent to eq. ()2.55|) . The quantity Ar is the finite combination of loop diagrams 
and counterterms in ()2.89p . The prediction for My/ within the SM or the MSSM is obtained 
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from evaluating Ar in these models and solving eq. ()2.89|) . 




Ar{Mw,Mz,mt,...) }. (2.90) 



In practice, this can be done by an iterative procedure since Ar itself depends on M^y. 
The one-loop contributions to Ar can be written as 



Ar = Aa - ^Ap + {Ar\,^, (2.91) 

where Aa is the shift in the fine structure constant due to the light fermions of the SM, 
Aa oc logrrif (see the discussion in Sect. 12. l!^ . and Ap is the leading contribution to the 
p parameter from fermion and sfermion loops. The remainder part, (Ar)i.em, contains in 
particular the contributions from the Higgs sector. 

In the following we will discuss the status of the theoretical evaluation of Mw- After a 
brief review of the SM contribution, the additional MSSM corrections are described in more 
detail. 



2.5.1 SM corrections 

In the SM, the resuh for (VB) in eq. ^TE^ is 



The singular part of this equation involving the divergence A = 2/(4 — D)— 7 + log 4-7? (see 
App. El coincides, up to a factor, with the non-Abelian bosonic contribution to the charge 
counterterm in eq. ()2.9|) : 

" /A_,og^fU^^. (2.93) 



Extra non-standard vertex and box diagrams do not change the singular part; they contribute 
another finite term (^-B) non-standard- Together with eq. ()2.9p and eq. ()2.92|) . we obtain from 
eq. ()2.89|1 the following expression. 



Ar ^ n^(0)-tf^-^V^"^°^;^^^ 



+ 2 ^ ^ + f 6 + l-g^^H. ) + (^i?)non-standard , (2.94) 



where in the on-shell renormalization the renormalization constants are given by the on-shell 
self-energies, as specified in eq. (j2.7p . 

Beyond the complete one- loop result [109], resummations of the leading one- loop contri- 
butions Aa and Ap are known [124]. They correctly take into account the terms of the form 



46 



(Aa)^, (Ap)^, (AaAp), and (AaArrem) at the two-loop level and the leading powers in Aa 
to all orders. 

Higher-order QCD corrections to Ar are known at 0{aas) [110,125,126] and 0{aal) [111, 
127] since about 10 years. Recently the full electroweak two-loop result for Ar has been 
completed. It consists of the fermionic contribution [90,128,129], which involves diagrams 
with one or two closed fermion loops, and the purely bosonic two-loop contribution [130]. 

Beyond two- loop order, besides higher-order contributions to Ap (see Sect. 12.4. 2j) the 
results for the pure fermion-loop corrections (i.e. contributions containing n fermion loops at 
n-loop order) are known up to four-loop order [131]. They contain in particular the leading 
contributions in Aa and Ap. 

Since the full result for is rather lenghty and contains numerical integrations of 
integrals appearing in the electroweak two-loop contributions, a simple parametrization is 
given in Ref. [132]. It approximates the full result for Mw to better than 0.5 MeV for 
10 GeV < Mh < 1 TeV if the other parameters are varied within their combined 2a region 
around their experimental central values. 

The expected size of the unknown higher-order corrections, i.e. the estimated theory 
uncertainties [132] (for Mh ^ 300 GeV) are summarized in Tab. 12.11 (see Ref. [65,132,133] 
for further details). 



2-loop 


3-loop 


C(a^, ferm) 


C(a^bos) 


0{aal, ferm) 


0(G2aXM|) 




compl. [90,128,129] 


compl. [130] 


compl. [111,127] 


3.0 


1.5 



4-loop 


0(G,afm?) 




1.3 


1.4 



Table 2.1: Estimated uncertainties from unknown higher-order corrections to Mw in 
MeV [132]. 

Currently these intrinsic uncertainties result in [132] 

^^SM,intr (p^^^^gn^) _ 4 MeV . (2.95) 

It seems reasonable that the evaluation of further higher-order corrections will lead to a 
reduction of this uncertainty by a factor of two or more on the timescale of 5-10 years. We 
therefore estimate as future intrinsic uncertainty 

^^SM,intr (f^t^^g) ^ 2 MeV . (2.96) 

The dominant theoretical uncertainty at present is the uncertainty induced by the ex- 
perimental errors of the input parameters. The most important uncertainties arise from the 
experimental error of the top-quark mass and the hadronic contribution to the shift in the 
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fine structure constant. The current errors for rrit [134] and Aahad [93] induce tlie following 
parametric uncertainties 

5m™* = 4.3 GeV AM^""'""' (current) ^ 26 MeV, (2.97) 
5(Aa™*) = 36 X 10-^ ^^para,Aa,.d (curreut) ^ 6.5 MeV. (2.98) 

At the ILC, the top-quark mass will be measured with an accuracy of about 100 MeV [7-9]. 
The parametric uncertainties induced by the future experimental errors of rrit and Acthad [135] 
will then be [63] 

^^futurc _ ^ AM^""'""' (future) ^ 1 MeV, (2.99) 

(5(Aaf",7'=) = 5 X 10-^ ^^para,Aa,.d (future) ^ 1 MeV. (2.100) 

Thus, the precision measurement of the top-quark mass at the ILC and prospective improve- 
ments in the determination of Aahad (see the discussion in Ref. [135]) will reduce the para- 
metric uncertainties to the same level as the prospective intrinsic uncertainties, eq. ()2.96|) . 
allowing a very sensitive test of the electroweak theory. 



2.5.2 SUSY corrections 

In this subsection we review the current status of the SUSY corrections to M^. The intrinsic 
uncertainties from missing higher-order SUSY corrections will be discussed in Sect. 13.1.21 



One-loop corrections 

The complete one-loop corrections to ArsusY ^ were evaluated independently by two groups [136, 
137]. The main part of the contributions stems from the t/b doublet that enters at the one- 
loop level only via gauge-boson self-energies. Therefore, only Feynman diagrams as depicted 
in Fig. 12. II have to be evaluated, but contrary to Sect. 12.4.21 also with no n- vanishing external 
momentum. In general, all scalar-quark contributions (yielding S'''^(0) = 0, according to 
the comment after eq. ()2.10|) ) are contained in 

^r^-^^m-'f(^-m^^^^^- (-01) 

In the approximation of neglecting the external momenta in the self-energies the second term 
in eq. ()2.1()H) reduces to Ap, leading to a decomposition as in eq. ()2.91|1 . For loops of scalar 
quarks the corrections mainly arise from the contribution to Ap (see eq. fl2.64|l ). so that Ar^ 
can be approximated as 

Ar« ~ — ^ Ap ^ -3.5Ap. (2.102) 



The full one- loop result from the t/h sector is compared with this approximation in Figs 
and 12.61 The case of no-mixing in the h sector is shown in Fig. 12.51 for tan/5 = 1.6 and 
Xt = 0, 200 GeV. The full result is reproduced by the Ap approximation within a few per 



^From here on we drop the subscript 
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Figure 2.5: The t/b corrections to Ar at the one-loop level, eq. ()2. 10111 . are compared with 
the approximation, eq. (j2.1(J2|) . The results are shown as a function of mq{= Msusy) for 
tan/3 = 1.6, Xb = and Xt = 0, 200 GeV. 
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Figure 2.6: The tjh corrections to Ar at the one-loop level, eq. ()2.101|1 . are compared with 
the approximation, eq. ()2.102p . The results are shown as a function of mq{= Msusy) for 
tan /3 = 40, X^ = 2500 GeV and Xt = 0, 200 GeV. 



49 



cent. The same applies for large mixing in the b sector, see Fig. 12. 6t with Xj, = 2500 GeV, 
tan /3 = 40 and = 0, 200 GeV. 

As investigated in detail in Refs. [136,137], the full SUSY one-loop contribution to Ar 
does not exceed (9(0.0015) (explicit formulas for the self-energy contributions are given 
in Ref. [138], see also Refs. [103,139].) The main contribution is given by the universal 
corrections, see eq. ()2.94|) . The corrections beyond the t/b sector arise from the other scalar 
quarks, entering only in the universal corrections, and the sleptons and gauginos, entering 
in the universal as well as in the non-universal contributions [136]. 

The full one-loop corrections from third and second generation squarks in the NMFV 
MSSM, using eq. (I2!l(l1jl . have been derived in Ref. [25]. 



Corrections beyond one-loop 

Since the dominant one- loop corrections are given by the t/b contributions, the existing 
two-loop calculations have focused on this sector. The only existing two-loop calculation for 
Ar, going beyond the Ap approximation as presented in Sect. 12.4.21 are the gluon-exchange 
corrections of 0{aas) [140]. This is the only result in the t/b sector beyond one- loop order 
that can be obtained as an analytical formula due to the presence of the massless gluon in 
the two-loop two-point function. The gluino-exchange corrections, on the other hand, have 
been shown to decouple for large nig [97], see Sect. 12.4.21 

The 0{aas) gluonic corrections are evaluated from the Feynman diagrams as shown in 
Fig. 12. 3t but taking into account the momentum dependence. Furthermore the derivative of 
the photon self-energy is needed. It is given by (D = 4 — 26) 
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77 TT 4s?, 



W 



f=l,2 a=l,2 
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The most complicated parts are the gauge boson self-energies with non-zero external mo- 
mentum. The general case is given by 
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a as 3 
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with 
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The functions ^40, Bq [141], T'234/, Ti23'4, 7ii23'4 and Ti23'45 [126,142] can be found in the 
Appendix. 
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2.6 Evaluation of Z boson observables 



The measurement of the Z-boson mass from the Z hneshapc at LEP has provided us with 
an additional precise input parameter besides a and G^. Other observable quantities from 
the Z peak allow us to perform precision tests of the electroweak theory by comparison 
with the theoretical predictions given by specific models. At the Z boson resonance in e'^e" 
annihilation, two classes of precision observables are available: 

a) inclusive quantities: 

• the partial leptonic and hadronic decay width Fy, 

• the total decay width F^, 

• the hadronic peak cross section (T/j, 

• the ratio of the hadronic to the electronic decay width of the Z boson, R^., 

• the ratio of the partial decay width for Z — > cc (66) to the hadronic width, Rcip) ■ 

b) asymmetries and weak mixing angles: 

• the forward-backward asymmetries Ap^, 

• the left-right asymmetries A(^^, 

• the r polarization P^-, 

• the effective weak mixing angle sin^ ^eflf- 

All these quantities can be written in a transparent way with the help of effective vector 
and axial vector couplings, which comprise the genuine electroweak loop contributions, be- 
sides those from the QED virtual-photon corrections, which are the same in the SM and in 
supersymmetric extensions. 



2.6.1 The effective Zff couplings 

The structure of the resonating Z amplitude allows us to define neutral-current (NC) vertices 
at the Z peak with effective coupling constants Qy j^, equivalently to the use oi pf,Kf: 



-.NC 



N 1/2 



Ml 



1/2 



(2.105) 



The complete expressions for the effective couplings read as follows: 

1 - Ar 



^ m^(M|) 

Vf -\- zswCwQf: 



1/2 



1 + nT(M|) 



1 - Ar 
+ ri^(M|) 



1/2 



(2.106) 
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Besides Ar, the building blocks are the following finite combinations of two-point functions 
evaluated at s = M|: 

z 

si, V Ml cw Ml ) ' 

(2.107) 

and the finite form factors -FV,a at s = M| from the vertex corrections a|i^'^°°^^ (including 
the external- fermion wave function renormalizations), 

^l--'""' - 2^ (-"^v-(^) - 7„7.F-(.) . - 0.) ^) . (2.108) 



For the explicit expressions for the self-energies and the vertex corrections including the 
MSSM contributions see Refs. [138,143-145]. Owing to the imaginary parts of the self- 
energies and vertices, the form factors and the effective couplings, respectively, are complex 
quantities. 

Effective Mixing Angles. We can define effective mixing angles for a given fermion species / 
according to 

sin^g. = . , i , 1 1 -Re^ I , (2.109) 

from the effective coupling constants in ()2.106p . They are of particular interest since they 
determine the on-resonance asymmetries. A special case is the effective mixing angle for the 
light leptons (/ = i), which is commonly denoted as the effective mixing angle (assuming 
lepton universality), 

sin^ ^eff = sin^ , (2.110) 
as, e.g., in the analysis of experimental data from LEP and SLC [146]. 

2.6.2 Z boson observables 

From lineshape measurements one obtains the parameters Mz, Tz, ctq, or the partial widths. 
Here Mz will be used as a precise input parameter, together with a and G^; the width and 
partial widths are specific model predictions. 

The total Z width Tz can be calculated as the sum over the partial decay widths 

Tz = J2^f, Tf = T{Z^ff) (2.111) 
/ 
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(other decay channels are not significant). The fermionic partial widths, when expressed in 
terms of the eff'ective coupling constants defined in (j2.1(J6|) . read 
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The photonic QED correction, given at one-loop order by 

3a 



QED 



(2.113) 



(2.114) 



is small, at most 0.17% for charged leptons. 

The factorizable SM (i.e. gluonic) QCD corrections for hadronic final states can be written 
as follows: 

Ar^cD = To (l gi r + I gi T) ^qcd , (2.115) 

where [147] 
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for the light quarks with ~ 0, with = as(M|). 

For h quarks the QCD corrections are different owing to finite h mass terms and to 
top- quark- dependent two-loop diagrams for the axial part: 
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(2.117) 



The coefficients in the perturbative expansions 

Ry 
Ra 



V"s , V f^s 

71 V vr 



as 

-) + 

vr 



c 



TT V TT 



a. 



+ 



depending on rrib and rrit, have been calculated up to third order in as, except for the rrib- 
dependent singlet terms, which are known to 0{a1) [148,149]. For a review of the QCD 
corrections to the width, with the explicit expressions for Ry^A, see Ref. [150]. 

The partial decay rate into b quarks, in particular the ratio Rb = Tb/T^ad, is an observable 
with special sensitivity to the top quark mass. Therefore, beyond the pure QCD corrections. 
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the two-loop contributions of the mixed QCD-electroweak type are also important. The QCD 
corrections were first derived for the leading term of 0{asG^mJ) [151] and were subsequently 
completed by the 0{as) correction to the Xogmt/ M-^y term [152] and the residual terms of 
0{aas) [153]. 

At the same time, the complete two- loop 0{aas) corrections to the partial widths for 
decay into the light quarks have also been obtained, beyond those that are already contained 
in the factorized expression ()2.115|1 with the electroweak one- loop couplings [154]. These 
"non-factorizable" corrections yield an extra negative contribution of —0.55(3) MeV to the 
total hadronic width. 

Besides the standard gluonic QCD corrections, there are supersymmetric QCD correc- 
tions involving virtual gluinos and squarks, which turned out to be very small [155,156], 
for masses of the SUSY partners in accordance with the bounds from direct experimental 
searches. 

From the partial widths and the total width 1)2.1111) the following set of combinations can 
be formed, 

the hadronic peak cross section, with the hadronic width Fhad = 'Yliq 
the ratio of the hadronic to the electronic decay width. 



i?e = ^ , (2.119) 



e 



the ratio of the partial decay width for Z ^ bb (cc) to the total hadronic decay width, 

Rho = ^ ■ (2.120) 



had 



The various asymmetries depend on the ratios of the vector to the axial vector coupling and 
thus on the effective mixing angles defined in eq. ()2.109|) . in terms of the combinations 

Af - 2(l-4|g,|sin^g,) 

-l + (l-4|Q,|sin^^,)2' ^'-'"'^ 

yielding 

the left-right asymmetry and the r polarization, 

Alr = ^% Pr = A\ (2.122) 
the forward-backward asymmetries, 

4B = l^''Af- (2.123) 
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Final-state QCD corrections, in the case of quark pair production, are important for the 
forward-backward asymmetries, at the one-loop level given by 

^ (l - ^) , (2.124) 

in the absence of cuts. Finite-mass effects have to be considered for b quarks only; they are 
discussed in Ref. [157]. Two-loop QCD corrections in the massless approximation are also 
available [158]. The SUSY-QCD corrections again turn out to be small for realistic values 
for squark and gluino masses [156]. 

2.6.3 The effective leptonic mixing angle 

Since sin^ 9^^ is a precision observable with high sensitivity for testing the electroweak theory, 
we discuss in this section the status of the theoretical predictions for sin^ 6^^ in the SM and 
the MSSM. 

SM corrections 

Recently the complete result for the fermionic two- loop corrections has been obtained [159], 
improving the prediction compared to the previously known (9(G^m^M|) term [160]. Con- 
trary to the case of the VT-boson mass, see Sect. 12.5. H the purely bosonic two- loop corrections 
are not yet completely known. 

Beyond two- loop order, the same kind of corrections are known as for Mw, i-e. QCD 
corrections of 0{aas) [110,125,126] and 0{aal) [111,127], pure fermion-loop corrections up 
to four-loop order [131], and three-loop corrections entering via Ap (see Sect. I2.4.2|) . 

A simple parametrization of the SM result for sin^ ^gfr containing all relevant higher-order 
corrections can be found in Ref. [159]. It reproduces the exact calculation with a maximal 
deviation of 4.5 x 10^^ for 10 GeV < Mh < 1 TeV if the other parameters are varied within 
their combined 2a region around their experimental central values. 

The estimated theory uncertainties for different parts of the unknown higher-order cor- 
rections are summarized in Tab. 12.21 fsee Refs. [159,65] for further details). 



2-loop 


3-loop 


0(0^, ferm) 


0{a^, hos) 


ferm) 


0(G>,m2M|) 




compl. [159] 


1.2 


compl. [111,127] 


2.3 


2.5 



4-loop 






1.1 


2.4 



Table 2.2: Estimated uncertainties from unknown higher-order corrections to sin^^efr 
in [10-5] [159,65] 
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Currently these intrinsic uncertainties result in [159] 

6 sin^ ^SM,intr (^^^^gj^^) ^ 5 X 10-5 . (2.125) 

In the future, an improvement down to about 

6 sin^ 0^SM,intr (f^^^^g) ^ 2 X 10"^ (2.126) 

seems achievable. 

Concerning the parametric uncertainties, the current errors for rrit [134] and Aa^ad [93] 
give rise to 

5m™* = 4.3 GeV A sin^ ^p^'^'^''"* (current) ^ 14 x 10-^ (2.127) 
5(Aa™*) = 36 X 10-^ A sin^ ^P^'^"'^"'^-^ (current) ^ 13 X 10-^ (2.128) 

The parametric uncertainties induced by the future experimental errors of rrit and Aahad 
are [63] 

^^future _ ^ Asin^^P™' (future) ^ 0.4 x 10-^ (2.129) 

5(A4"^*7^) = 5 X 10"^ A sin^^P^''^'^^''^'^ (future) ^ 1.8 X 10"^ (2.130) 

Compared to the GigaZ accuracy (see Tab. II. 4p on sin^ O^s also the parametric uncertainty 
induced by the experimental error of Mz is non-negligible [63] 

6Mz = 2.1 MeV =^ A sin^ ^p^'^'^^^ ^ 1.4 x 10"^ (2.131) 

As in the case of Mw, the precision measurement of the top-quark mass at the ILC and 
prospective improvements in the determination of Aa^ad will reduce the parametric uncer- 
tainties to the same level as the prospective intrinsic uncertainties, eq. ()2.126|) . 



MSSM corrections 

As for M]y, the largest correction to sin^ O^s in the MSSM can be expected from scalar quark 
contributions. The shift in sin^ O^s is then given by 

A sin^ els = ^2^'^ Ar'" - swCwfV^{Ml), (2.132) 

with 



fp-(Ml) = _ * f ^ _ (2.133) 

^ ^' Ml a„ V Ml Ml, J' ^ ' 

and Ar^ from eq. ()2.101|) . 

In the MSSM the complete one-loop corrections to sin^ 9es have been evaluated as de- 
scribed in Sect. 12.6.21 Beyond one-loop order the leading term can be included via the p pa- 
rameter approximation, eq. ()2.58p . where Ap at the two-loop level is given in Sect. 12. 4T^ 
The intrinsic uncertainties from missing higher-order SUSY corrections will be discussed in 
Sect. Em 

The full one-loop corrections from third and second generation squarks in the NMFV 
MSSM, using eq. (I2!i;-i2ll . have been derived in Ref. [25]. 
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2.7 The lightest Higgs boson mass as a precision ob- 
servable 



A striking prediction of the MSSM is the existence of at least one hght Higgs boson. The 
search for this particle is one of the main goals at the present and the next generation 
of colliders. Direct searches at LEP have already ruled out a considerable fraction of the 
MSSM parameter space [12,13]. With the forthcoming data from the Tevatron, the LHC 
and the ILC either a light Higgs boson will be discovered or the MSSM will be ruled out as 
a viable theory for physics at the weak scale. Furthermore, if one or more Higgs bosons are 
discovered, their masses and couplings will be determined with high accuracy at the ILC. 
Thus, a precise knowledge of the dependence of the masses and mixing angles of the MSSM 
Higgs sector on the relevant supersymmetric parameters is of utmost importance to reliably 
compare the predictions of the MSSM with the (present and future) experimental results. 

The Higgs sector of the MSSM has been described in Sect. II. 2. H at tree-level, leading to 
the prediction for the lightest MSSM Higgs boson, m/i_trcc < ^z, see eq. fjl.l2|) . However, 
this mass bound, which arises from the gauge sector of the theory, is subject to large radiative 
corrections in particular from the Yukawa sector of the theory [101]. Because of the impor- 
tance of the higher-order corrections, a lot of effort has been devoted to obtain higher-order 
results in the MSSM Higgs sector. Results for the complete one-loop contributions are avail- 
able [161,144]. Corrections beyond one-loop order have been obtained with different methods. 
Leading and subleading two-loop corrections have been obtained in the Effective Potential 
(EP) approach [162,163], the Renormalization Group (RG) improved EP approach [164], and 
with the Feynman-diagrammatic method [165-167]. Detailed comparisons of the different 
methods have been performed [168,169]. The higher-order corrections shift the upper bound 
on ruh to TJih < 136 GeV [166,170] (for mt = 178 GeV and MgusY < 1 TeV, neglecting 
uncertainties from unknown higher-order corrections). 

In the case that the MSSM parameters possess non- vanishing complex phases (cMSSM), 
the upper bound on remains the same as for the MSSM with real parameters, but the 
Higgs-boson couplings can vary significantly compared to the case with real parameters. 
Complex phases are possible for the trilinear couplings, Af,f = t,b,T, . . ., for the Higgsino 
mass parameter, n, and for the gaugino mass terms. Mi, M2, M3 = rrig (where one of the 
latter ones can be rotated away by a redefinition of the corresponding fields). Recently 
the different methods for the evaluation of higher-order corrections in the MSSM Higgs 
sector have even been extended to the cMSSM, reaching nearly the precision as in the real 
MSSM [171-175]. In the following, however, we will focus on the real case. 

2.7.1 Higher-order corrections to ruh 

The tree- level bound on rrih, being obtained from the gauge couplings, receives large correc- 
tions from SUSY-breaking effects in the Yukawa sector of the theory. The leading one-loop 
correction is proportional to mf. The leading logarithmic one-loop term (for vanishing mix- 
ing between the scalar top quarks) reads [101] 




(2.134) 
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Corrections of this kind have drastic effects on the predicted value of rrih and many other 
observables in the MSSM Higgs sector. The one-loop corrections can shift rrih by 50-100%. 

In the Feynman diagrammatic (FD) approach the higher-order corrected Higgs boson 
masses are derived by finding the poles of the h, if-propagator matrix. Its inverse is given 
by 

^^"^^'^ -~\ t^p') p'-rnl... + ^..ip') )^ ^'-'''^ 

where the S(p^) denote the renormalized Higgs-boson self-energies (see eq. ()2.38|l ). and p 
is the external momentum. Determining the poles of the matrix Aniggs in eq. ()2.135p is 
equivalent to solving the equation 

^ = 0. (2.136) 

The status of the available results for the self-energy contributions to eq. (|2.135|) can 
be summarized as follows. For the one-loop part, the complete result within the MSSM is 
known [101,161]. The by far dominant one-loop contribution is the 0{at) term due to top 
and stop loops {at = y^/{4:n), where yt has been defined in eq. ()2.75p ). 

The evaluation of two-loop corrections is quite advanced and it has now reached a stage 
where all the presumably dominant contributions are known. They include the strong cor- 
rections, usually indicated as 0{atas), and Yukawa corrections, 0{af), to the dominant one- 
loop 0{at) term, as well as the strong corrections to the bottom/sbottom one-loop 0{ab) 
term {at = ?/^/(47r)), i.e. the 0{ab<ys) contribution. The latter can be relevant for large values 
of tan/?. Presently, the (!?(«*«,) [162,164-166], 0{af) [162,164,176,177], 0{abas) [178,96], 
0{ataf,), 0{al) [179] contributions to the self-energies are known for vanishing external mo- 
menta. In the (s)bottom corrections the all-order resummation of the tan /^-enhanced terms, 
C(ab(Q;s tan/?)"), is also performed [98,107] (see Sect. 12.2."^ . The above results have been 
implemented into the program FeynHiggs [180,99,181], which evaluates observables in the 
MSSM Higgs sector (including also results with complex phases). 

Recently, also the full electroweak two-loop corrections in the approximation of vanishing 
external momentum [163] and the leading two-loop momentum dependent effects [182] have 
been published. For these corrections no public code is available yet. In order to apply this 
result for expressing nih in terms of physical masses, a transition of the parameters Mz and 
Ma in Refs. [163,182] to their on-shell values will be required at the two- loop level. 

Besides the masses of the Higgs bosons, also their couplings are affected by large higher- 
order corrections. For the MSSM with real parameters, leading corrections can conveniently 
be absorbed into the couplings by introducing an effective mixing angle aes- It is obtained 
from the higher-order corrected Higgs-boson mass matrix in the approximation where the 
momentum dependence of the Higgs-boson self-energies is neglected. 

The Higgs-boson mass matrix in the 0i-02 basis reads in this case 

/M-2 ^ ( - (0) ml^^^ - ±^,^M \ 
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P^ -"^Ktree + ^hhiP^ 



P 



m 



H.trec 



+ f^HH{p^ 



^hH{P^ 



where the Es(0) {s = 0i,0i02;02) denote the renormahzed Higgs-boson self-energies (in 
the 01, 02 basis), including one- and two- loop (and possibly higher-order) corrections. These 
self-energies (at zero external momentum) have to be inserted into eq. (j2.137p . Diagonalizing 
this higher-order corrected Higgs-boson mass matrix 



Higgs 



yields the effective mixing angle aes- 



"off 



m 



H 







(2.13^ 



fteff = arctan 



-(Ml + M|) sin/3 cos/3 - S^,^, 
M| cos2 p +Ml sin^ ^ -t^,- ml 



TT TT 
— < Cteff < — 
2 2 



(2.139) 



Replacing in the Higgs-boson couplings the tree-level mixing angle a by the higher- 
order corrected effective mixing angle Oefr leads to the inclusion of the leading higher-order 
corrections that enter via Higgs-boson propagator corrections [183,184]. 



2.7.2 Remaining intrinsic and parametric uncertainties 

If the MSSM is realised in nature, the light CP-even Higgs-boson mass will be measured 
with high precision at the next generation of colliders. The prospective accuracies for a light 
SM-like Higgs boson that can be obtained in the experimental determination of mh at the 
LHC [185] and at the ILC [7-9] are: 

5ni^^ ^ 200 MeV (LHC) , (2.140) 
(5m^"P ^ 50 MeV (ILC) . (2.141) 

Since ruh depends sensitively on the other sectors of the MSSM, in particular on the t sector 
(see eq. ()2.134|l ). the light CP-even Higgs-boson mass will be very important for precision 
tests of the MSSM. 

The remaining theoretical uncertainties in the prediction for rrih have been discussed in 
Refs. [170,63,227,186]. For recent reviews on the current status of the theoretical prediction 
see also Refs. [175,187]. 

We begin with the discussion of the parametric uncertainties. Since the leading one- 
loop corrections to ruh are proportional to the fourth power of the top quark mass, the 
predictions for rrih and many other observables in the MSSM Higgs sector sensitively depend 
on the numerical value of m^. As a rule of thumb [188], a shift of 6mt = 1 GeV induces a 
parametric theoretical uncertainty of rrih of also about 1 GeV, i.e. 

^^para,m, _ ^ ^2.142) 

The uncertainties induced by the experimental error of rrit at the LHC [189] and the 
ILC [7-9], 

(JmfP ^ 1-2 GeV (LHC) , (2.143) 

^m^^ ^ 0.1 GeV (ILC) , (2.144) 
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can be compared with the parametric uncertainties induced by the other SM input parame- 
ters. Besides m^, the other SM input parameters whose experimental errors can be relevant 
for the prediction of rrih are M^y, a^, and rrib- The W boson mass enters only in higher orders 
through the quantum corrections to muon decay (since Gp is used for the parametrization, 
see eq. (EHHl))- 

The present experimental error of SM^^ = 34 MeV leads to a parametric theoretical 
uncertainty of rrih below 0.1 GeV. In view of the prospective improvements in the exper- 
imental accuracy of Mw, the parametric uncertainty induced by Mw will be smaller than 
the one induced by m^, even for Srrit = 0.1 GeV. 

The current experimental error of the strong coupling constant, 6as{Mz) = 0.002 [3], 
induces a parametric theoretical uncertainty of rrih of about 0.3 GeV. Since a future im- 
provement of the error of as{Mz) by about a factor of 2 can be envisaged [3,62,190], the 
parametric uncertainty induced by rrit will dominate over the one induced by as{Mz) down 
to the level of 5mt = 0.1-0.2 GeV. 

The mass of the bottom quark currently has an experimental error of about 6mb = 
0.1 GeV [3,191]. A future improvement of this error by about a factor of 2 seems to be feasi- 
ble [191]. The influence of the bottom and sbottom loops on rrih depends on the parameter 
region, in particular on the values of tan/5 and /i (the Higgsino mass parameter). For small 
tan/? and/or /i the contribution from bottom and sbottom loops to rrih is typically below 
1 GeV, in which case the uncertainty induced by the current experimental error on rrif, is 
completely negligible. For large values of tan/3 and fi, the effect of bottom/sbottom loops 
can exceed 10 GeV in rrih [178,170,96]. Even in these cases we find that the uncertainty 
in ruh induced by 6mb = 0.1 GeV rarely exceeds the level of 0.1 GeV, since higher-order 
QCD corrections effectively reduce the bottom quark contributions. Thus, the parametric 
uncertainty induced by rut will in general dominate over the one induced by m^, even for 
Srrit^ 0.1 GeV. 

The comparison of the parametric uncertainties of rrih induced by the experimental errors 
of Mw, cts{Mz) and rrib with the one induced by the experimental error of the top quark mass 
shows that an uncertainty of 6mt ~ 1 GeV, corresponding to the accuracy achievable at the 
LHC, will be the dominant parametric uncertainty of ruh. The accuracy of 6mt ~ 0.1 GeV 
achievable at the ILC, on the other hand, will allow a reduction of the parametric theoretical 
uncertainty induced by 6mt to about the same level as the uncertainty induced by the other 
SM input parameters. 

We now turn to the intrinsic theoretical uncertainties in the prediction for rrih from 
unknown higher-order corrections. Even if all the available higher-order corrections described 
above are taken into account, the intrinsic uncertainty in from unknown higher-order 
corrections is still estimated to be quite substantial [99,170]^. The numerical relevance of 
the unknown higher-order corrections depends on the region of MSSM parameter space that 
one considers. An overall estimate of the intrinsic uncertainty can therefore be only a rough 
guidance for "typical" MSSM parameter regions. In regions where higher-order corrections 
are particularly enhanced (for instance very large mixing in the stop sector or regions where 

■^For codes that do not include all the existing higher-order corrections the intrinsic theoretical uncertain- 
ties can be much larger. 
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the bottom Yukawa coupling is close to being non-perturbative) the theoretical uncertainties 
can be significantly larger. 

At the two-loop level, various genuine electroweak two-loop corrections from different 
sectors of the MSSM are not yet included in the publicly available codes. A rough estimate 
of their numerical impact can be obtained from the relative importance of the correspond- 
ing contributions at the one-loop level. This has been performed in Ref. [170] and yielded 
an estimate of the remaining uncertainty of unknown two-loop corrections of about 2 GeV. 
Another way of estimating the effect of unknown two-loop corrections is to apply differ- 
ent renormalization schemes at the one-loop level and to vary the renormalization scale of 
quantities that are renormalized according to the DR scheme [99]. As an example for the 
latter approach, Fig. 12.71 shows the effect of varying the renormalization scale that enters via 
the renormalization of tan/5 and the Higgs field renormalization constants at the one-loop 
level for "typical" MSSM parameters (see caption). The corresponding shift in the one- loop 
result for rrih, which is of the order of genuine two-loop corrections that are not included in 
the current prediction for ruh, is indicated by the grey areas. The uncertainty in ruh from 
varying /ijjj^ from mt/2 to 2mt is in accordance with the above estimate of the uncertainty 
from missing two-loop corrections of about ±2 GeV. 
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Figure 2.7: The renormalization scale dependence of m/j introduced via a DR definition 
of tan/3 and the Higgs field renormalization constants is shown as a function of Ma (left 
plot) and tan/3 (right). The lower curves correspond to tan/? = 2 (left) and Ma = 100 GeV 
(right). For the upper curves we have set tan/? = 20 (left) and Ma = 500 GeV (right), 
/ig^ has been varied from mt/2 to 2mt. The other parameters are Msusy = 500(1000) GeV, 
Xt = 2 Msusy, M2 = = 500 GeV. The dotted line corresponds to a full on-shell scheme, 
for more details see Ref. [99]. 

Beyond two-loop order, corrections that effectively shift the value of the top-quark mass 
entering the calculation are particularly important because of the sensitive dependence of rrih 
on rrit. Corrections of this kind of 0{ata'^) can be estimated by varying the renormalization 
scheme of the top-quark mass at the two-loop level. Another possibility for estimating the size 
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of three-loop corrections is to analyze the numerical impact of the leading logarithmic three- 
loop term that can easily be obtained with renormalization group methods [170,192]. Both 
possibilities have been investigated in detail in Ref. [170], yielding an estimate of the intrinsic 
theoretical uncertainty beyond the two-loop level of about 1.5 GeV. Similar strategies in 
the case of the O(afea^) correction [96] lead to an intrinsic uncertainty of up to 2 GeV in 
the case of /x < (in regions where the effects of the bottom/sbottom sector are strongly 
enhanced), and of about ~ 100 MeV for /i > 0. 

As an overall estimate for the current intrinsic uncertainties in the prediction of ruh we 
obtain 

Am^;^^' (current) = 3 GeV. (2.145) 

On the timescale of 5-10 years it seems reasonable to expect that the complete two-loop 
calculation (which is already technically feasible with the currently existing tools) can be 
incorporated into efficient codes and that the higher-order uncertainties can be reduced by 
at least a factor of two, leading to the estimate 

Am^;^^' (future) = 0.5 GeV. (2.146) 

2.7.3 Higgs sector corrections in the NMFV MSSM 

Within the MSSM with MFV, the dominant one-loop contributions to the self-energies 
in ()2.137j) result from the Yukawa part of the theory (i.e. neglecting the gauge couplings); 
they are described by loop diagrams involving third-generation quarks and squarks. Within 
the MSSM with NMFV, the squark loops have to be modified by introducing the generation- 
mixed squarks, as given in Sect. 11.2.61 The leading terms are obtained by evaluating the 
contributions to the renormalized Higgs-boson self-energies at zero external momentum, 
Ss(0),s = hh,hH,HH. The evaluation has been restricted to the dominant Yukawa con- 
tributions resulting from the top and t/t (and c/c) sector. Corrections from b and b/b (and 
s/s) could only be important for very large values of tan/5, tan/5 ^ mt/rrib and have not 
been considered so far. The analytical result of the renormalized Higgs boson self-energies, 
based on the general 4x4 structure of the t/c mass matrix, has been derived in Ref. [25]. 
However, as has also been shown in Ref. [25], the corrections for are not significant for 
moderate generation mixing. 



2.8 The anomalous magnetic moment of the muon 

Another observable which is important in the context of precision tests of the electroweak 
theory is the anomalous magnetic moment of the muon, a^ = {g — 2)^/2. For the interpre- 
tation of the results in the context of Supersymmetry (or other models of new physics) 
the current status of the comparison of the SM prediction with the experimental result is 
crucial, see Refs. [193,194] for reviews and the discussion in Sect. 11.3."^ It currently results 
in a deviation of [75] 

- = (25.2 ± 9.2) x 10"^° : 2.1 a . (2.147) 
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2.8.1 MSSM one-loop calculation 

The anomalous magnetic moment of the muon is related to the photon-muon vertex 
function F^^^p as follows: 

u{p')T,Mp^-p'^^Xp) = uip')bpFviQ') + iP + P')pFMiq') + ---]uip), (2.148) 

= -2m^FM{0). (2.149) 

It can be extracted from the regularized vertex function using the projector [195,196] 

1 „ iD-2 



a/1 



2{D - 1){D - 2)ml 



Til -—- [ml^p - Dp^^ -{D- l)m^pp] V 



TO 

+ ^ {i + m,) (7.7p - 7p7.) (# + m^) T"^ }>, (2-150) 



Vp = r^M-(p,-p,0), (2.151 
d 

Tpu = -Q^T^-^A-{P - (g/2), -p - (g/2), q) 



(2.152) 

5=0 



Here the muon momentum is on-shell, p^ = m^, and D is the dimension of space-time. For 
more details see Refs. [195-197]. 

The complete one-loop contribution to can be devided into contributions from dia- 
grams with a smuon-neutralino loop and with a sneutrino-chargino loop, see Fig. 12. 8[ leading 
to 

^^susY,iL _ Aaf + Aaff^ (2.153) 




Figure 2.8: The generic one-loop diagrams for the MSSM contribution to a^: diagram with 
a sneutrino-chargino loop (left) and the diagram with a smuon-neutralino loop (right). 

The full one- loop expression can be found in [198], see Ref. [199] for earlier evaluations. 
If all SUSY mass scales are set to a common value, 

-^susY = fn^± = "^x° =fTT'p = "^i>M (2.154) 
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the result is given by 

^susY.iL ^ ^3 ^ f tan/5 sign(/.) . (2.155) 

V ^SUSY / 

Obviously, supersymmetric effects can easily account for a (20 . . . 30) x 10~^° deviation, if 
is positive and Msusy lies roughly between 100 GeV (for small tan j3) and 600 GeV (for 
large tan/?). Eq. (j2.155|) also shows that for certain parameter choices the supersymmetric 
contributions could have values of either a^^^^ ^ 55 x 10"^'' or a^^^'*' ^ — 5 x 10~^°, both 
outside the Sex band of the allowed range according to eq. ()2.147p . This means that the 
[g — 2)n measurement places strong bounds on the supersymmetric parameter space. 

2.8.2 MSSM two-loop calculation 

In order to fully exploit the precision of the {g — 2)^ experiment within SUSY, see e.g. 
Refs. [200-203] for discussions of the resulting constraints on the parameter space, the the- 
oretical uncertainty of the SUSY loop contributions from unknown higher-order corrections 
needs to be under control. It should be significantly lower than the experimental error given 
in eq. (|1.52j) and the hadronic uncertainties in the SM prediction, leading to the combined 
uncertainty given in eq. ()2.147|) . 

For the electroweak part of the SM prediction the desired level of accuracy has been 
reached with the computation of the complete two-loop result [195,196], which reduced the 
intrinsic uncertainty from QED and electroweak effects below the level of about 1 x 10~^° [75]. 
For the SUSY contributions, a similar level of accuracy has not been reached yet, since 
the corresponding two-loop corrections are partially unknown. Four parts of the two-loop 
contribution have been evaluated up to now that will be reviewed in the next subsections. 

Two-loop QED corrections 

The first part are the leading log (m^/MgusY)-terms of supersymmetric one-loop diagrams 
with a photon in the second loop. They are given by [204] 

^^SUSY,2L,QED ^ ^^SUSY.IL ^ (^^) ) " (2-156) 

They amount to about —8% of the supersymmetric one-loop contribution (for a SUSY mass 
scale of Msusy = 500 GeV). 

Two-loop Two-Higgs-doublet contributions 

In the MSSM, the bosonic electroweak two-loop contributions differ from the SM because 
of the extended MSSM Higgs sector. This class is defined by selecting all MSSM two-loop 
diagrams without a closed loop of fermions or sfermions and without pure QED-diagrams, 
see the first line in Fig. 12.91 The results presented in this section have been obtained in 
Ref. [197]. 
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Figure 2.9: Some MSSM two-loop diagrams for with (depending on the diagram) F — 
li,u„ /, /' = t, 6, T, (t>^KH, A, H^, G, G^; ^ = G^, H^; /, /' = i, 6, f, i>,; 1/ = 7, ^, W; 

X\,2i Xl,2,3,4- 
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The result a^°^'^^MSSM ^^^^g 



a 



bos,2L,MSSM 



5 Gf,ml a 



bos,2L,MSSM 



Cr 



where the coefficient of the logarithm is given by 



log 
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fj. I bos,2L,MSSM 
+ Co ' 



bos,2L,MSSM 



30 

C/3 



[98 + 9c^ + 23(1 - Aslrf] 



m 



H 



mi 



(2.157) 



(2.158) 
(2.159) 



Here Cq, = cos a, etc. Using the tree-level relations in the Higgs sector, it can be shown that 
1, and thus the logarithms in the SM and the MSSM are identical. The coefficient 



c, 



bos,2L,MSSM 



is more complicated and not given here, see Ref. [197]. 



Two-loop corrections with a closed SM fermion/sfermion loop 

The third known part are the diagrams with a closed loop of SM fermions or scalar fermions 
calculated in Ref. [205], extending previous results of Refs. [206,207]. 

The two-loop diagrams discussed in this subsection can be subdivided into three classes 
(all diagrams are understood to include the corresponding subloop renormalization): 
(fVcp) diagrams with a sfermion {t, b, f, 9^) loop, where at least one gauge and one Higgs 
boson is exchanged, see the second line of Fig. 12.91 

{fVV) diagrams with a sfermion loop, where only gauge bosons appear in the second loop, 
see the third line of Fig. 12. 9| 

(/1/0) diagrams with a fermion (t, b, r, z/,-) loop, where at least one gauge and one Higgs 
boson are present in the other loop, see the fourth line of Fig. 12.91 The corresponding 
diagrams with only gauge bosons are identical to the SM diagrams and give no genuine SUSY 
contribution. The difference between the SM and the MSSM originates from the extended 
Higgs sector of the MSSM. Diagrams where two Higgs bosons couple to the external muon 
are suppressed by an extra factor of m^/M^ and hence negligible. 

The counterterm diagrams contain the renormalization constants 6M^^, SZ^., Sth^n 
corresponding to mass, charge and tadpole renormalization and can be easily evaluated. 
For the evaluation the on-shell renormalization scheme has been chosen. This leads to 
SM^rz = R-eS'^'^(M^^), where S'^'^ denote the transverse parts of the gauge-boson self- 
energies, see Sect. Em The charge renormalization is given by SZe = —1/2 H'''(0), see 
eq. ()2.9|) . The tadpoles are renormalized such that the sum of the tadpole contribution T 
and the counterterm vanishes, i.e. Sth^n = —Th,H, see Sect. 12.1.51 

Numerically the most important contribution comes from the diagrams with a Higgs 
boson and photon exchange. This type of contributions can be particularly enhanced by 
the ratio of the mass scale of the dimensionful Higgs-Sfermion coupling divided by the mass 

2 

scale of the particles running in the loop, i.e. by ratios of the form {/i. A, -^j}/{my, nih^H}, 
which can be much larger than one. For large tan (3 and large sfermion mixing, the leading 
terms are typically given by the parts of the couplings with the highest power of tan f3 and 
by the loop with the lightest sfermion. These contributions involve only if-exchange, since 
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the /i-couplings approach the SM-Higgs couphng for not too small Ma- They can be well 
approximated by the formulas [205] 



AafL = -0.013 X 10-^0 m^/itan/^ ^-^^^^^^^ ^2.160) 



mb Afotan /? 



Aa°'^^ = -0.0032 x lO"^'^ sign(/i), (2.161) 

where and are the masses of the lighter t and b, respectively, and tuh is the mass 
of the heavy CP-even Higgs boson. The formulas holds up to few percent if the respective 
sfermion mass fulfils m^i ^ mn- Since the heavier sfermions also contribute and tend to 
cancel the contributions of the lighter sfermions, these formulas do not approximate the full 
result very precisely, but they do provide the right sign and order of magnitude. 



Two-loop contributions with a closed chargino/neutralino loop 

The 2-loop contributions to containing a closed chargino/neutralino loop constitute a 
separately UV-finite and gauge-independent class and have been evaluated in Ref. [197]. 
Corresponding diagrams are shown in the last line of Fig. 12.91 

The chargino/neutralino two-loop contributions, a^'^^, depend on the mass parameters 
for the charginos and neutralinos /i. Mi 2, the CP-odd Higgs mass Ma, and tan/3. It is 
interesting to note that, contrary to Ref. [119], no tree- level relations in the Higgs sector 
were needed in order to find a UV-finite result. This is due to the fact that each two- loop 
diagram contributing to {g — 2)^ together with its corresponding subloop renormalization is 
finite. 

The parameter dependence of a^'^'" is quite straightforward [197]. If all supersymmetric 
mass scales are set equal, fi = M2 = Ma = Msusy (with the only exception that Mi = 
5/3 s^/c^ M2), the approximate leading behaviour of a^'^^ is simply given by tan/5/M|usY, 
and the following relation holds, 

aJ^-«llxlO-»(^)(^^) sig„(,). (2.162) 

As shown in Ref. [197], the approximation is very good except for very small Msusy and 
small tan /?, where the leading term is suppressed by the small /i, and subleading terms begin 
to dominate. 



Remaining intrinsic uncertainties 

So far, at the two-loop level, the MSSM corrections to the Two-Higgs-Doublet model (THDM) 
one-loop diagrams have been evaluated. The only exception here are the diagrams that con- 
tain as a second loop an additional closed smuon-neutrino or muon-sneutrino loop. However, 
these corrections are expected to be small. 

The remaining two-loop corrections that are not yet available are 
— the contributions with a mixed SM fermion/sfermion loop attached to a SUSY one-loop 
diagram. 
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— the full THDM corrections to the SUSY one-loop diagrams. This will include as a subset 
also the QED corrections evaluated in Ref. [204], where, however, all SUSY masses had been 
set equal to Msusy- 

The first missing class of mixed SM fermion/sfermion contributions might in principle be 
as large as the SM fermion or scalar fermion corrections obtained in Ref. [205], see above. 
This leaves an intrinsic uncertainty of about ~ 3 x 10~^°. The second class gives corrections 
smaller than 10% to the MSSM one-loop result. Assuming that the corresponding intrinsic 
uncertainties are less than half of the evaluated corrections, the combined effect of the 
unknown two-loop corrections can be estimated to be about 

Aa'/^^' (current) = 6 x 10"^° . (2.163) 

After a full two-loop calculation will be available, the intrinsic theoretical uncertainty from 
unknown QED and electroweak higher-order corrections should be at the level of 

AaJ*'' (future) = 1 x 10"^° . (2.164) 

2.9 Tools and codes for the evaluation of electroweak 
precision observables 

The large number of different fields in the MSSM gives rise to a plethora of possible interac- 
tion vertices. Calculations at the one-loop level and beyond therefore usually involve a lot 
of Feynman diagrams. The diagrams in general contain several mass scales, making their 
evaluation (in particular beyond one- loop order) increasingly difficult. Since the necessary 
steps can be structured in a strictly algorithmic way, they can be facilitated with the help 
of computer algebra tools and numerical programs. 

Computer algebra tools have heavily been used in deriving the results discussed above. 
Because of the multitude of scales involved in SUSY higher-order corrections, in most cases 
the result cannot be expressed in a compact form. Instead, the results presented above have 
been transformed into public computer codes (also being used for the numerical evaluation 
in Sects. ElandlH). 

2.9.1 Tools for the calculation of EWPO 

The calculation of higher-order SUSY Feynman- diagrams consists of several steps. First the 
topologically different diagrams for the given loop order and the number of external legs need 
to be generated. Inserting the fields of the model under consideration into the topologies 
in all possible ways leads to the Feynman diagrams. The Feynman rules translate these 
graphical representations into mathematical expressions. Since the loop integrals in general 
lead to divergences, the expressions need to be regularized and renormalized. The evaluation 
of the Feynman amplitudes involves a treatment of the Lorentz structure of the amplitude, 
calculation of Dirac traces etc. At the one-loop level it is possible to reduce all tensor inte- 
grals to a set of standard scalar integrals, which can be expressed in terms of known analytic 
functions. In contrast to the one-loop case, no general algorithm exists so far for the eval- 
uation of two-loop corrections in the electroweak theory. The main obstacle in two-loop 
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calculations in massive gauge theories is the complicated structure of the two-loop integrals, 
which makes both the tensor integral reduction and the evaluation of scalar integrals very 
difficult. In general the occurring integrals are not expressible in terms of polylogarithmic 
functions [208]. For the evaluation of some types of integrals that do not permit an analytic 
solution numerical methods and expansions in their kinematical variables have been devel- 
oped. Computer-algebraic methods can facilitate most of the above-mentioned steps. There 
are computer algebra packages available based on FORM [209], Mathematica [210] or both. 

A package for the generation of SUSY amplitudes and drawing the corresponding dia- 
grams is FeynArts [121,211]. As a feature of particular importance for higher-order calcula- 
tions in the electroweak theory, FeynArts generates not only the unrenormalized diagrams 
at a given loop order but also the counterterm contributions at this order and the countert- 
erm diagrams needed for the subloop renormalization. For one-loop calculations with up 
to four external legs (the inclusion of five external legs is currently under way) the package 
FormCalc [212] can be used, where for numerics the LoopTools [213] package can easily be 
linked. For the evaluation of two-loop diagrams with up to two external legs the program 
TwoCalc [214] can be used. It is based on an algorithm for the tensor reduction of gen- 
eral two-loop 2-point functions and can be used for an automatic reduction of Feynman 
amplitudes for two-loop self-energies with arbitrary masses, external momenta, and gauge 
parameters to a set of standard scalar integrals. The above computer algebra codes evaluate 
the multi-loop diagrams analytically without performing expansions for small parameters 
etc. 

The program QGRAPH [215] is an efficient generator for Feynman diagrams (so far 
restricted to the SM, see however Ref. [216]). As output the diagrams are encoded in a 
symbolic notation. Being optimized for high speed, QGRAPH is particularly useful for 
applications involving a very large number (i.e. (9(10^)) of diagrams. Its output, depending 
on the number of scales and external legs can then be passed to MATAD [217], MINCER [218] 
or EXP [219], where expansions for small parameters are performed. 

An alternative package for SM and SUSY one- loop calculations is the GRACE sys- 
tem [220]. 

Overviews about codes for higher- loop and -leg calculations can be found in Refs. [221-223]. 

2.9.2 Public codes for the numerical evaluation of EWPO 

The results presented in Sects. 12.41 12.5l and l2.6l have been implemented in the code POMSSM"^, 
which has been used for the numerical evaluation in Sects. El and 0] The Higgs boson sector 
evaluations have been done with the code FeynHiggs [180,99,175,181], including the correc- 
tions described in Sect. 12.71 This code also performs an evaluation of all Higgs boson decay 
widths as well as production cross sections for photon colliders. Also the results for Ap as 
described in Sect. 12.41 are included as a subroutine. Other codes for evaluations of Higgs 
sector observables are Hdecay [224] and CPsuperH [225]. The results for the anomalous 
magnetic moment of the muon, described in Sect. 12. 8[ are available as a subroutine for the 
code FeynHiggs. 



new version of POMSSM is currently prepared and will be available from the authors. 
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Chapter 3 

MSSM predictions versus 
experimental data 

Now we study the impact of the higher-order corrections to the electroweak precision observ- 
ables discussed above. The MSSM predictions are compared with the current experimental 
results and constraints on the parameter space of the unconstrained MSSM are discussed. 
We furthermore investigate how the improved electroweak precision measurements at the 
next generation of colliders enhance the sensitivity of testing the electroweak theory. 

3.1 MSSM predictions for Mw and sin^ 6^^ 

3.1.1 Numerical analysis in the MSSM 
Results for Ap 

We start our discussion of the numerical results with the quantity Ap, which parametrizes 
leading SUSY contributions to the VT-boson mass and the Z-boson observables, see Sect. 12.41 
The effect of the gluonic SUSY two-loop contributions as given in eq. ()2.73|) (the four squark 
masses are renormalized on-shell; the mass shift arising from the SU(2) relation is understood 
to be absorbed into the one- loop result, see Sect. I2.4.2j) is shown for an exemplary case in 
Fig. 13.11 as a function of MgusY- The other parameters are tan/5 = 3 and Xt = 0,2Msusy- 
The line for Xt = 2MgusY starts only at MgusY ~ 300 GeV. For lower values of Msusy one of 
the scalar top mass squares is below zero. Ap^^p^^^^ can reach values of up to 0.2 x 10~^. The 
results for the gluino-exchange contribution are shown in Fig. 13.21 fX^ = 0) and Fig. l3.3l fXt = 
2Msusy) for m-g = 0, 10, 200, 500 GeV (and m-g = 800 GeV in the latter) as a function of 
Msusy- The results for rrig = and 10 GeV are indistinguishable for Xt = 0. The decoupling 
for large nig is visible already for nig = 500 GeV. In the case of Xt = 2Msusy5 see Fig. 13.31 
^Pi^fufno is general positive and can reach values up to 0.5 x 10~^ for nig = 200 GeV. 
As can be seen in the figure, for larger values of rrig the contribution to Ap decouples as 
expected. Contrary to the SM case where the strong two-loop corrections screen the one-loop 
result, the 0{aas) corrections in the MSSM increase the one-loop contributions by up to 
35%, thus enhancing the sensitivity to scalar quark effects. 

In Fig. 13.41 the numerical result of the leading 0{at) MSSM corrections in the limit of 
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Figure 3.1: Ap^ gf^^^^ as a function of the common squark mass Msusy for tan P = 3, Xt, = 
andXt = 0,2MsusY. 
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Figure 3.2: Ap^^^f^fj^^^ as a function of the common squark mass Msusy for tan P 



3, Xf, = 0, 



Xt = and rrig 



0, 10 (the curves are indistinguishable), 200, 500 GeV [97]. 



large MgusY (see Sect. I2.4.2j) is shown. It is compared with the other contributions to Ap: the 
0(q!j) SM correction (with Mfl^ = rrih) and the SUSY contributions from the scalar quark 
sector at 0{a) and 0{aas)- The results are shown as a function of MgusY? which enters the 
0{al) corrections indirectly via its effect on rrih- For small tan/5 and Ma = 300 GeV, see 
the left plot of Fig. 13.41 the effective change arising from the new genuine MSSM corrections 
compared to the 0{af) SM contribution with Mf^^ = rrih is sizable. While the full 0{a'^) 
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Figure 3.3: Ap^^fjf^^ as a function of the common squark mass Msusy for tan (3 = 3, Xb = 0, 
Xt = 2MsusY and m-g = 0, 10, 200, 500, 800 GeV. 



result is larger than the 0{a) corrections for Msusy ^ 600 GeV, it is larger than the 0{aas) 
corrections for all Msusy- However, the genuine MSSM corrections are always smaller than 
the MSSM 0{a) contributions but they are of equal size as the 0{aas) corrections for 
Msusy ~ 600 GeV. (Note that for smaller Msusy the approximation of neglecting the 
scalar-quark contributions in the 0{af) result may no longer be valid.) Since they enter 
with a different sign into Ap, they can compensate each other. Similar results are found in 
the no-mixing case, which is not shown here. 

The case of large tan j3 and Ma = 300 GeV is shown in the right plot of Fig. 13.41 The 
curve for the 0{af) MSSM corrections in the limit Msusy ^ cxd is indistinguishable in the 
plot from the 0{af) SM contribution with Mf^ = rrih- The difference between these two 
corrections is approximately 1.5 x 10~^, while the 0{aas) corrections are about 10~^ even 
for Msusy = 1000 GeV. The purely electroweak corrections decouple much faster for large 
tan/5 than the 0{aas) corrections, see also eq. ()2.78p . 

In Fig. 13.51 we show the result for the 0{af), 0{atab), and 0{al) MSSM contributions 
to Ap in the m^^^ and the no-mixing scenarios [32,33] (see also App. |B)) compared with the 
corresponding SM result with = rrih. In the left plot tan/3 is fixed to tan/3 = 40 and 
Ma is varied from 50 GeV to 1000 GeV. In the right plot Ma is fixed to Ma = 300 GeV, 
while tan [3 is varied. 

For large tan/5 the 0{atab) and 0{al) contributions yield a significant effect caused by 
the heavy Higgs bosons in the loops, entering with the other sign than the 0{af) corrections, 
while the contribution of the lightest Higgs boson is SM-like. As one can see in Fig. 13.51 for 
large tan/5 the MSSM contribution to Ap is smaller than the SM value. For large values of 
Ma, the SM result is recovered. 
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Figure 3.4: The contribution of the leading 0{af) MSSM corrections in the hmit of large 
MgusY, '^Pi^^^s) is shown as a function of MgusY for Ma = 300 GeV and tan/? = 3 (left 
plot) or tan /? = 40 (right plot) in the case of the m™'^'' scenario, see App. IHI Apf^^^^ is 
compared with the leading 0{a^) SM contribution and with the leading MSSM corrections 
originating from the i/b sector of 0{a) and 0{aas)- Both 0{a^) contributions are negative 
and are for comparison shown with reversed sign. In the right plot the 0{a^) corrections 
differ by about 1.5 x 10^'', which is indistinguishable in the plot. 
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Figure 3.5: The 0{a^^), 0{atab), and 0{al) MSSM contributions to Ap in the m^^^ and 
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Quality of the Ap approximation 

We now turn to the numerical effects on Mw and sin^ 9es induced by Ap. As a first step 
the quality of the Ap approximation, using eq. ()2.58p . is analyzed [140]. We show the 
comparison of the Ap approximation with the full evaluation at the two-loop level, where such 
a calculation is available. As described in Sect. 12.5.21 only the two-loop gluonic corrections 
to Ar have been calculated so far. In Figs. IH.fi^ IH.7I we show the full gluonic two-loop 
contribution to Ar together with the corresponding Ap approximation. The no-mixing case 
in the b sector is presented in Fig. 13.61 with tan/3 = 1.6 and Xf = 0,200 GeV. The case 
with Xb = 2500 GeV is shown in Fig. O with tan/3 = 40 and Xt = 0, 200 GeV. As for the 
one-loop case, see Figs. 12.51 12.61 also in the two-loop case the Ap approximation reproduces 
the full result to better than 10%. 

Corrections to Mw and sin^ O^s induced by Ap 

We illustrate the effects of the corrections to Ap discussed above on the observables Mw 
and sin^ ^efr for the example of the O(a^) MSSM contributions in the limit of large Msusy- 

Fig. ESI shows the shift 6Mw induced by the 0{a^) MSSM contribution for MgusY = 
1000 GeV in the m™'^^ scenario, see App. El The other parameters are p = 200 GeV, = 
At. rrih is obtained in the left (right) plot from varying Ma from 50 GeV to 1000 GeV, while 
keeping tan (3 fixed at tan /? = 3, 40 (from varying tan j3 from 2 to 40, while keeping Ma fixed 
at Ma = 100,300 GeV). Besides the absolute 0{a'^) MSSM contribution (solid and short- 
dashed lines) also the "effective change" compared to the SM is shown, i.e. the difference 
between the 0{a^) MSSM contribution and the 0{af) SM contribution with M^^ = mt 
(long-dashed and dot-dashed lines). While the full result shows contributions to Mw of up 
to 11 MeV, the effective change is much smaller, mostly below the level of 2 MeV. 

For large tan/3 the 0{atah) and 0{al) contributions yield a significant effect from the 
heavy Higgs bosons in the loops, entering with the other sign than the 0{af) corrections, 
while the contribution of the lightest Higgs boson is SM-like, see Sect. 12.4.21 The effective 
change in the predictions for the precision observables from the 0{atab) and 0{al) correc- 
tions can exceed the one from the 0{a^) corrections. It can amount up to 6Mw ~ +5 MeV 
for tan/3 = 40. 

Fig. 13.91 shows the shift 5sin^6'efr induced by the absolute O^a^) MSSM contribution 
(solid and short-dashed lines) and the effective change (long-dashed and dot-dashed lines) for 
-^susY = 1000 GeV in the m'^^^ scenario. The other parameters are p = 200 GeV, Af, = At. 
rrth is obtained in the left (right) plot from varying Ma from 50 GeV to 1000 GeV, while 
keeping tan (3 fixed at tan /3 = 3, 40 (from varying tan (3 from 2 to 40, while keeping Ma 
fixed at Ma = 100,300 GeV). While the full result shows contributions to ?,ir? 6^^ of up to 
6 X 10^^, the effective change is much smaller, mostly below the level of 1 x 10^^. 

For large tan /3, the effective change in the predictions for the precision observables from 
the 0{atai,) and 0{al) corrections can exceed the one from the 0{a'^) corrections. It can 
amount up 5sin^ 6^^ ^ —3 x 10~^ for tan/3 = 40. 
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change compared with the SM result are shown for Msusy = 1000 GeV in the m™*^^ scenario. 
The other parameters are /i — 200 GeV, Ai, = A^. rrih is obtained in the left (right) plot 
from varying Ma from 50 GeV to 1000 GeV, while keeping tan P fixed at tan /3 = 3, 40 (from 
varying tan/3 from 2 to 40, while keeping Ma fixed at Ma = 100, 300 GeV). 



0.00006 
0.00005 
0.00004 

(D 
i 

- 0.00003 

5 
:> 

0.00002 
0.00001 
0.00000 









max 




tanp = 3 '- 




tanp = 40 - 




tanp = 3 (MSSM - SM) 




tanp = 40 (MSSM - SM) - 







110 120 

m, [GeV] 



0.00002 
0.00001 
0.00000 
-0.00001 



1 1 , 1 




















M^= 100 GeV 






- 300 GeV 






Mj- 100 GeV (MSSM 


SM) - 




- 300 GeV (MSSM 


SM) ; 








-, , , , 1 


. . < 1 . < < . 1 < . . < 1 . 



10 



20 30 

tanp 



40 



50 



Figure 3.9: The shift S sin^ ^eff induced by the 0{a^) MSSM contribution and the effective 
change compared with the SM result are shown for Msusy = 1000 GcV in the m™^^ scenario. 
The other parameters are /j, = 200 GeV, Ai, = A^. nih is obtained in the left (right) plot 
from varying Ma from 50 GeV to 1000 GeV, while keeping tan (3 fixed at tan /3 = 3, 40 (from 
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77 



MSSM predictions for M^r and sin^ in comparison with present and future 
experimental precisions 

Now we focus on the comparison of the Mw and sin^ 6es prediction with the present data 
and the prospective experimental precision at the next generation of colliders. 

In Fig. l3.1Ul we compare the SM and the MSSM predictions for as a function of m^. 
The predictions within the two models give rise to two bands in the rrit-Mw plane with 
only a relatively small overlap region (indicated by a dark-shaded (blue) area in Fig. IH.lOl) . 
The allowed parameter region in the SM (the medium-shaded (red) and dark-shaded (blue) 
bands) arises from varying the only free parameter of the model, the mass of the SM Higgs 
boson, from = 113 GeV (upper edge of the dark-shaded (blue) area) to 400 GeV (lower 
edge of the medium-shaded (red) area). The light-shaded (green) and the dark-shaded 
(blue) areas indicate allowed regions for the unconstrained MSSM. SUSY masses close to 
their experimental lower limit are assumed for the upper edge of the light-shaded (green) 
area, while the decoupling limit with SUSY masses of 0{2 TeV) yields the lower edge of 
the dark-shaded (blue) area. Thus, the overlap region between the predictions of the two 
models corresponds in the SM to the region where the Higgs boson is light, i.e. in the 
MSSM allowed region (m/j < 140 GeV). In the MSSM it corresponds to the case where 
all superpartners are heavy, i.e. the decoupling region of the MSSM. The current 68% C.L. 
experimental results^ for rrit and M^y slightly favor the MSSM over the SM. The prospective 
accuracies for the LHC and the ILC with GigaZ option, see Tab. 11.41 are also shown in the 
plot (using the current central values), indicating the potential for a significant improvement 
of the sensitivity of the electroweak precision tests [62]. 

In Fig. 13. Ill the comparison between the SM and the MSSM is shown in the Mwsin^ 9eH 
plane (see also Refs. [226,227]). As above, the predictions in the SM (medium-shaded and 
dark-shaded (red and blue) bands) and possible MSSM regions (light-shaded and dark- 
shaded (green and blue) bands) are shown together with the current 68% C.L. experimental 
results and the prospective accuracies for the LHC and the ILC with GigaZ option. Again the 
MSSM is slightly favored over the SM. It should be noted that the prospective improvements 
in the experimental accuracies, in particular at the ILC with GigaZ option, will provide a 
high sensitivity to deviations both from the SM and the MSSM. 

The central value for the experimental value of sin^ 9es in Fig. 13.111 is based on both 
leptonic and hadronic data. The two most precise measurements, Alr from SLD and Ap^ 
from LEP, differ from each other by about 3a (see Ref. [18]). This, together with the NuTeV 
anomaly (see below), gave rise to a relatively low fit probability of the SM global fit in the 
past years, and had caused considerable attention in the literature. In particular, several 
analyses have been performed where the hadronic data on A-pB have been excluded from 
the global fit (see e.g. Refs. [228,229]). It has been noted that in this case the SM global 
fit, possessing a significantly higher fit probability, yields an upper bound on Mh which 
is rather low in view of the experimental lower bound on Mh of Mh > 114.4 GeV [13]. 
The value of sin^ 9^^ corresponding to the measurement of y4LR(SLD) alone is sin^ ^^cfr = 
0.23098 ± 0.00026 [18]. Fig. EHU shows that adopting the latter value of sin"^ 9es makes the 
agreement between the data and the SM prediction much worse, while the MSSM provides 
a very good description of the data. In accordance with this result, in Ref. [229] it has 

^The plot shown here is an update of Refs. [136,226,227]. 
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Figure 3.10: The current experimental results for M^y and rrit and the prospective accuracies 
at the next generation of colliders are shown in comparison with the SM prediction (medium- 
shaded and dark-shaded (red and blue) bands) and the MSSM prediction (light-shaded and 
dark-shaded (green and blue) bands). 



been found that the contribution of light gauginos and scalar leptons in the MSSM (in a 
scenario with vanishing SUSY contribution to Ap) gives rise to a shift in M\y and sin^ 6'eff 
as compared to the SM case which brings the MSSM prediction in better agreement with 
the experimental values of and y4LR(SLD). 

On the other hand, it has also been investigated whether the discrepancy between Alr 
and could be explained in terms of contributions of some kind of new physics. The 
(loop-induced) contributions from SUSY particles in the MSSM are however too small to 
account for the 3a difference between the two observables (see e.g. Ref. [229]). Thus, the 
quality of the fit to Alr and Ap^ in the MSSM is similar to the one in the SM. 

With the latest experimental values of the precision observables and the most up-to-date 
theory predictions the probability of the global fit in the SM is about 26% [18] (if the NuTeV 
result is not included). In particular, the most recent experimental value of the top-quark 
mass, rrit = 178.0 ± 4.3 GeV [134], and a slight shift in the experimental value of M^r have 
led to an improvement of the overall fit quality. Although the discrepancy between Alr 
from SLD and Ap^ from LEP remains, it seems not well motivated to discard any of the two 
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accuracies at the next generation of colliders are shown in comparison with the SM prediction 
(medium-shaded and dark-shaded (red and blue) bands) and the MSSM prediction (light- 
shaded and dark-shaded (green and blue) bands). 



measurements. 

As mentioned above, another observable for which the SM prediction shows a large 
deviation by about 3a from the experimental value is the neutrino-nucleon cross section 
measured at NuTeV [230]. Also in this case loop effects of SUSY particles in the MSSM are 
too small to account for a sizable fraction of the discrepancy (see e.g. Refs. [231,232]). 

3.1.2 Intrinsic uncertainty in Mw and sin^^eff from SUSY correc- 
tions 

The remaining theoretical uncertainties in the prediction for Mw and sin^ ^eff from unknown 
higher-order corrections in the MSSM (i.e. loop corrections from SM particles and superpart- 
ners) are considerably larger than in the SM, since the results for higher-order corrections in 
the MSSM are not quite as advanced yet as in the SM. The current intrinsic uncertainties 
in the MSSM can roughly be estimated by comparing the size of the known corrections in 
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the MSSM (see above) to the corresponding corrections in the SM and by assuming that 
the unknown higher-order corrections in the MSSM enter with the same relative weight as 
the corresponding corrections in the SM, whose numerical effects are known. This kind of 
estimate does not take into account specific enhancement factors in the MSSM, like for in- 
stance corrections that grow with powers of tan f3. In general, the additional contributions 
from superpartners in the loops will be bigger the smaller the SUSY mass scale is. As in the 
case of nih, the estimate for the intrinsic uncertainty of Mw and sin^ 9es should be under- 
stood to refer to "typical" regions of the MSSM parameter space. In parts of the parameter 
space where certain corrections are particularly enhanced (see the discussion in Sect. 12. 2|) 
the intrinsic uncertainties can be larger. 

Taking the above considerations into account, a crude estimate of the current intrinsic 
uncertainties yields [227] 

MSSM: 6M'^'' (current) = 10 MeV , ^sin^^;"^*"^ (current) = 12 x 10"^ , (3.1) 

i.e. uncertainties that are roughly twice as large as the current uncertainties in the SM. 

With sufficient effort on higher-order calculations in the MSSM, it should be possible in 
the future to reduce the intrinsic uncertainties to the same level as we had estimated for the 
SM (see eqs. (ITT??)|l . ^UT^ ): 

MSSM: 6M'^'' (future) = 2 MeV , ^sin^ 9'^^' (future) = 2 x 10"^ . (3.2) 
3.1.3 Results in the NMFV MSSM 

The analytical results obtained for the EWPO in the NMFV MSSM have been derived for 
the general case of mixing between the third and second generation of squarks, i.e. all NMFV 
contributions, ^ll,lr,rl,rr^ can be chosen independently in the t/c and in the h/s sector, 
see Sect. 11.2.61 Corrections from the first-generation squarks are not considered, for reasons 
discussed in Sect. 11.2.61 The numerical analysis of NMFV effects for the EWPO, however, 
have been performed for the simpler, but well motivated, scenario where only mixing between 
II and cl as well as between hi and sl is considered. The only fiavor off-diagonal entries in the 
squark-mass matrices are normalized according to A^^ = A'-^Mg^Mg^, following [21,23,24] ^, 
where Mg^ g^ are the soft SUSY-breaking masses for the SU (2) squark doublet in the third 
and second generation. NMFV is thus parametrized in terms of the dimensionless quantities 
A* and A^ (see [23,24,233,234] for experimentally allowed ranges). The case of A* = A^ = 
corresponds to the MSSM with minimal flavor violation (MFV). In detail, it has been set 

= A*M^ M^^ , Ai^ = A*,^ = A*^^ = , 

A^^ = \'M-^M~^^ , A^^ = A^^ = A^^ = , (3.3) 

for the entries in the matrices ()1.44|) and ()1.45|) . 

For the sake of simplicity, the same flavour mixing parameter has been assumed in the 
numerical analysis for the i — c and b — s sectors, A = A* = A''. It should be noted in 
this context that the LL blocks of the up-squark and down-squark mass matrices are not 

^The parameters A' and introduced here are denoted by ((S]5i)23 and {Sf 1^)23 in [21,23,24]. 
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independent because of the SU (2) gauge invariance; they are related trough the CKM mass 
matrix [24], which also implies that a large difference between these two parameters is not 
allowed. 

Results for Ap 

For the numerical evaluation [25], the m™^ and the no-mixing scenario have been used [33], 
but with a free scale Msusy, see App. |B| The results are independent of Ma- The numerical 
values of the SUSY parameters are 

Msusy = 1 TeV and 2 TeV, tan p = 30, fi = 200 GeV, e = 0.04, (3.4) 

if not explicitly stated otherwise. The variation with fi and tan/? is very weak, since they 
do not enter the squark couplings to the vector bosons. 

The behaviour with the sign of e is shown in Fig. 13.121 for the corrections to Ap^ as 
a function of A(= A* = A**). The results are shown for different relative signs of e and 
A, choosing A > 0, and fixing |e| = 0.04. Msusy has been set to Msusy = 2 TeV. For 
the m™''^ scenario the effect is small, but in the no-mixing scenario the results are affected 
significantly by the sign of e. The squark contribution to Ap^ can become of O{10^^) for 
A > 0.5. 
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Figure 3.12: The variation of Ap^ with A(= A* = A'') in the m^'^^ and no-mixing scenarios 
for different relative signs of e and A [25]). Msusy = 2 TeV, the other SUSY parameters are 
given in eq. (j3.4|) . 

In Fig. 13.131 we show the dependence of Ap^ on A(= A* = A^) for both the m™^^ and no- 
mixing scenario and for two values of the SUSY mass scale, Msusy = 1 TeV and Msusy = 
2 TeV. It is clear that Ap^ grows with the A parameter, being close to zero for A = and 
Msusy = 2 TeV. One can also see that the effects on Ap^ are in general larger for the 
no- mixing scenario (see also the results shown in Ref. [97]). For large values of Msusy the 
correction increases with increasing A since the splitting in the squark sector increases. 
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Figure 3.13: The variation of Ap^ with A = A* = in the m™^^ and no-mixing scenarios. 
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Figure 3.14: The variation of Ap^ with Msusy in the m^^^ scenario (left panel) and no- 
mixing scenario (right panel), for different values of A [25]. 



The behavior of the corrections with the SUSY mass scale is shown in Fig. 13.141 for 
different values of A in the m™**^ scenario (left panel) and in the no-mixing scenario (right 
panel). The region below Msusy ~ 400 GeV (depending on the scenario) implies too low and 
hence forbidden values for the squark masses. The curves are only for the allowed regions. 
For A = 0, Ap^ decreases, being zero for large Msusy values, in agreement with the results 
shown in Ref. [97]. We have also found that, for A 7^ and small values of Msusy? Ap^ 
decreases until it reaches a minimum and then increases for largest values of the SUSY scale. 
This increasing behavior is more pronounced for larger A values, reaching the level of a few 
per mill. The reason lies once again in the increasing mass splitting. 
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Numerical evaluation for M-^y and sin^ 9^^ 

Here the numerical effects of the NMFV contributions on the electroweak precision observ- 
ables, 5Mw and 5 sin^ ^eff, are briefly analyzed [25]. The shifts in My/ and sin^ 6^^ have been 
evaluated both from the complete expressions for the scalar quark contributions, eq. ()2.10ip 
and eq. ()2.132j) . and using the Ap^ approximation ()2.58p . The corrections to these two ob- 
servables based on eq. ()2.58|1 as a function of A(= A* = A**) are presented in Fig. IH. 151 with the 
other parameters chosen according to ()H.4j) . The m™^^ scenario and no- mixing scenario are 
selected for both plots, with two values of Msusy, as before. The induced shifts in Mw can 
become as large as 0.14 GeV for the extreme case, i.e. when Msusy = 2 TeV, A = 0.6 and 
the case of no-mixing is considered. In the ut^^^ scenario 5My/ is smaller, 5M\y ^ 0.05 GeV, 
but still sizeable. Using the complete expressions, eq. ()2.101|) and eq. ()2.132|1 . yields results 
practically indistinguishable from those shown in Fig. 13.151 [25] . Thus eq. ()2.58|1 is a suffi- 
ciently accurate, simple approximation for squark-mixing effects in the electroweak precision 
observables. 




A, X 

Figure 3.15: The variation of 5My/ and 5sin^ O^q as a function of A = A* = A^, for the m™^^ 
and no- mixing scenarios and different choices of Msusy obtained with eq. (j2.58p [25] . Using 
the complete expressions, eq. (|2.101|) and eq. (|2.132|) . yields practically indistinguishable 
results. 

The shifts 5 sin^ ^eff, shown in the right plot of Fig. 13.151 can reach values up to 7 x 10^"^ 
for MgusY = 2 TeV and A = 0.6 in the no-mixing scenario, being smaller (but still sizeable) 
for the other scenarios considered here. 

Extreme parts of the NMFV parameter space (especially for A* 7^ A^) can be excluded 
already with today's precision. But even small values of A = A* = A** could be probed with 
the future precision on sin^ 6'efr, provided that theoretical uncertainties will be sufficiently 
under control [227]. 
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3.2 The lightest MSSM Higgs boson mass 



The light CP-even MSSM Higgs boson mass, rrih, depends at tree-level on Ma and tan (3. Via 
loop corrections, see Sect. 12.7.11 it depends most strongly on the top quark mass and on the 
parameters of the scalar top sector. As an example, in Fig. I3.16l we show irih as a function of 
tan (3 in two benchmark scenarios, the m^^^ and the no- mixing scenario [33], see App. El 
is shown for a central value of mt = 178.0 GeV (dashed curves), and the variation with rrit 
by ±4.3 GeV is shown as the shaded (green) band. Higher rrih values are obtained for larger 
mt- (All results in this section have been obtained with FeynHiggs2.1 [99,175,180,181].) 

From the result for the m™*^^ scenario in Fig. 13.161 the upper bound of rrih ^ 136 GeV 
for mt = 178 GeV and Msusy = 1 TeV (neglecting the intrinsic theoretical uncertainties) 
can be read off that was mentioned in Sect. 12.71 Allowing a 1 a variation of mt shifts the 
upper bound on mh to about 140 GeV. The variation of the mh prediction with mt is even 
larger in the region of small tan j3. Fig. 13. 161 shows that a 1 a upward fluctuation of mt shifts 
the minimum of mh in the m™'^^ scenario to a value above 114 GeV. Thus, in this case the 
exclusion bound from LEP does not rule out any value of tan /3. The comparison of the 
MSSM prediction with the LEP exclusion bound is shown in more detail in Fig. l3.18l below. 
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Figure 3.16: m^ is shown as a function of tan/5 in the m™'^^ and the no-mixing scenario. 
mt has been varied in the interval mt = 178.0 ± 4.3 GeV. 

The relevance of the parametric uncertainty in m^ induced by different experimental 
errors on mt is emphasized in Fig. 13.171 [63], where the prediction for m/j is shown as a 
function of in the m™'^^ benchmark scenario. The evaluation of m^ has been done for a 
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central value of the top quark mass of = 175 GeV and for tan/3 = 5. The figure shows 
that a reduction of the experimental error from 6mf^^ = 1-2 GeV (LHC) to 5m^^^ = 0.1 GeV 
(ILC) has a drastic effect on the prediction for ruh. 

The prospective experimental error on m/j is also shown in Fig. I3.17[ while no intrinsic 
theoretical uncertainty from unknown higher-order corrections is included. If this intrinsic 
uncertainty can be reduced to a level of ^m)^*'^'^^*"'^'' 0.1 GeV, its effect in the plot would 
be roughly as big as the one induced by 6mf^'^ = 0.1 GeV. An intrinsic uncertainty of 
Srn^^^'^^^^'^'^ ~ 1 GeV, on the other hand, would lead to a significant widening of the band 
of predicted rrih values (similar to the effect of Smf^^ = 1 GeV). In this case the intrinsic 
uncertainty would dominate, implying that a reduction of 6mf^^ = 1 GeV to Smf'^ = 
0.1 GeV would lead to an only moderate improvement of the overall theoretical uncertainty 
of rrih. 
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Figure 3.17: rrih is shown as a function of Ma in the m™^'' scenario for tan/? = 5 [63]. 
Three different precisions for rrit are indicated (with a central value of rrit = 175 GeV). The 
anticipated experimental error on ruh at the ILC of 0.05 GeV is indicated by a horizontal 
band. 

Confronting the theoretical prediction for ruh with a precise measurement of the Higgs- 
boson mass constitutes a very sensitive test of the MSSM, which allows one to obtain con- 
straints on the model parameters. The sensitivity of the ruh prediction to Ma shown in 
Fig. 13.171 cannot directly be translated into a prospective indirect determination of Ma, 
however, since Fig. 13.171 shows the situation in a particular benchmark scenario [33] where, 
by definition, certain fixed values of all other SUSY parameters are assumed. In a realistic 
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situation the anticipated experimental errors of the other SUSY parameters and possible 
effects of intrinsic theoretical uncertainties have to be taken into account. In Sect. 13.51 the 
prospects for an indirect determination of SUSY parameters from precision physics in the 
MSSM Higgs sector will be discussed. 

As another example we demonstrate the impact of the current theory uncertainty of 
^^mtr ^ 2 GeV [170] on the exclusion bound of tan/3, see Ref. [188] for a detailed discussion. 
The m™^^ benchmark scenario [33] has been designed such that for fixed values of rrit and 
-^susY the predicted value of the lightest CP-even Higgs boson mass is maximized for each 
value of Ma and tan/5. In Fig. 13.1^ we show again rrih as a function of tan/5, together 
with the LEP exclusion bound for the mass of a SM-like Higgs [13], Mf^ > 114.4 GeV, as 
a vertical long-dashed line. The solid thick line shows the result in the m™*^^ scenario for 
mt = 178.0 GeV. 



10 - 



CO 



m^""^", FH2.1, mj = 178 GeV 
+ 3 GeV theory unc. 



, FH2.1 , Mg^gY = 2 TeV 
mj = (178+4.3) GeV 
+ 3 GeV theory unc. 



SM exclusion bound I 




80 90 100 110 120 130 

[GeV] 



140 



150 



Figure 3.18: rrih is shown as a function of tan/3 in the m™'^^ scenario (solid) and for 
mt = 178.0 + 4.3 GeV, Msusy = 2 TeV (dot-dashed). A theory uncertainty from unknown 
higher-order uncertainties of Sm^i^^^ = 3 GeV has been neglected (thick lines) or included 
(thin). The SM exclusion bound of Mh = 114.4 GeV, which for small values of tan/3 also 
roughly apphes for the MSSM, is indicated by a dashed fine. 
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While in general a detailed investigation of a variety of different possible production and 
decay modes is necessary in order to determine whether a particular point of the MSSM 
parameter space can be excluded via the Higgs searches or not, the situation simplifies 
considerably in the region of small tan (3 values. In this parameter region the lightest CP-even 
Higgs boson of the MSSM couples to the Z boson with SM-like strength, and its decay into a 
bb pair is not significantly suppressed. Thus, within good approximation, constraints on tan [3 
can be obtained in this parameter region by confronting the exclusion bound on the SM Higgs 
boson with the upper limit on nth within the MSSM. From the intersection of the theoretical 
upper bound in the m^^^ scenario (solid thick line) with the experimentally excluded region 
for nih the experimentally excluded region for tan/5 can be read off. For comparison we 
also show the same upper bound including the theory uncertainty from unknown higher 
order corrections, 5im}^^^ ^ 3 GeV [170] (solid thin line). Taking the theory uncertainty into 
account, the bound on tan/? is considerably weakened (see also Ref. [187]). Furthermore we 
show the m™*^^ scenario with the top-quark mass shifted upwards by one standard deviation, 
mt = 182.3 GeV and with Msusy = 2 TeV (dot-dashed thick line), also including the 3 GeV 
intrinsic theoretical uncertainty (dot-dashed thin line). Even without taking into account 
the intrinsic theoretical uncertainty, in this case no region of tan (3 can be excluded from the 
Higgs search at LEP. This example shows that both a reduction of the experimental error 
on rrit and of the intrinsic theoretical uncertainty will be crucial in order to obtain reliable 
bounds on the SUSY parameters from measurements in the Higgs sector (see also Sect. 13. 5|) . 

3.3 MSSM predictions for {g - 2)^ 

In our numerical discussion of SUSY contributions to the anomalous magnetic moment of the 
muon we first analyze the one-loop results from a scan over the MSSM parameter space [202] 
and then focus on two recently obtained two-loop corrections: the corrections involving a 
closed SM fermion/sfermion loop [205], and the ones involving closed chargino/neutralino 
loops [197]. 

3.3.1 One-loop results from a MSSM parameter scan 

The possible size of the MSSM one-loop contributions to can be assessed by a parameter 
scan. In Fig. 13.191 (from Ref. [202]) the possible MSSM contributions to are shown as a 
function of the lightest observable SUSY particle (LOSP). The lighter (green) dots corre- 
spond to a /i LOSP, darker (red) dots represent charginos/neutralinos as LOSP. The dashed 
lines show the allowed contours if {A^l is allowed to vary up to 100 TeV. The shaded bands 
correspond to the one/two a allowed ranges in the year 2001. One can see that the MSSM 
can easily explain the discrepancy in eq. (j2.147|) . On the other hand, can place stringent 
constraints on the allowed MSSM parameter space. In order to set reliable bounds in the 
MSSM the theoretical uncertainties have to be under control. This requires the evaluation 
of higher-order contributions. The existing two-loop corrections are reviewd in the following 
subsections. 
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Figure 3.19: MSSM one-loop contributions to are shown as a function of the mass of the 
hghtest observable SUSY particle (LOSP), obtained from a scan over the MSSM parameter 
space [202]. The lighter (green) dots correspond to a. ft LOSP, darker (red) dots represent 
charginos/neutralinos as LOSP. The dashed lines show the allowed contours if |v4^| is allowed 
to vary up to 100 TeV. The shaded bands correspond to the one/two a allowed ranges in 
the year 2001. 



3.3.2 Contributions from closed SM fermion/sfermion loops 

The two-loop corrections to {g — 2)^ involving a closed SM fermion/sfermion loop, corre- 
sponding to the diagrams in lines 2-4 of Fig. 12.91 have been evaluated in Ref . [205] , ex- 
tending earlier analyses of Refs. [206,207]. These two-loop corrections have a complicated 
parameter dependence. Therefore a parameter scan has been performed, tan/3 was set to 
tan /3 = 50, and universal soft SUSY-breaking parameters in the scalar fermion mass matri- 
ces were assumed. It turned out to be crucial to take experimental constraints from rrih, Ap, 
BR(6 57) and BR.{Bs — > fi^p^) into account (for details see Ref. [205]). It was shown 
that the diagrams involving a photon and a Higgs boson (diagram no. 12 in Fig. 12. 9p give 
the by far largest contribution. 

The whole contribution of this set of diagrams is shown in Fig. 13.201 The results shown 
in the figure are the following (see Ref. [205] for further details): 

• The outer lines show the largest possible results if all experimental constraints are 
ignored. They show a steep rise of Aa^^ for decreasing m^^; for mj_^ < 150 GeV 
contributions larger than 15 x 10"^°, corresponding to two standard deviations of the 
experimental error on a^, are possible. 

• The next two lines show the possible results if the bound ruh > 106.4 GeV (it results 
from the experimental bound of 114.4 GeV by taking into account a 5 GeV parametric 
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uncertainty from the experimental error of rrit and a 3 GeV intrinsic uncertainty, see 
Sect. 12.7. 2|) and then in addition the bound on Ap are satisfied. The maximum contri- 
butions are very much reduced already by the nih bound, and the Ap bound reduces 
further the positive region for small sfermion masses. If both bounds are taken into 
account, AaJ}^ > 5 x 10"^° and Aaf,^ < -10 x 10"^° is excluded for m? > 100 GeV. 

• The two innermost lines correspond to taking into account in addition the bound on 
BR(i?s and finally also on BR(6 — > 57), resulting in the shaded area. In 

particular taking into account the BR(6 — ^ 57) bound eliminates most data points 
with TUf^ 150 GeV and thus leads to a strong reduction of the possible size of the 
contributions (see however the discussion in Ref. [235]). The largest contributions of 
±4 X 10~^° to Aa^, corresponding to ~ 0.7a of the experimental error, are possible 
for m? ^ 150. . .200 GeV. 



15 



o 

'o 



05 
< 



I ' ' '''-^ ' ' ' ' I 
s 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

all data 

bound 



Ap bound 

-> (1(1 bound 

B->X7bound 



iMiiiiiMiiiiiMiiiiii allowed 
after all contraints 

I I I I I I I I I I I I I I I I 



50 100 150 200 250 300 350 400 450 500 

Min. sfermion mass [GeV] 



Figure 3.20: Maximum contributions of the diagrams with a closed SM fermion/sfermion 
loop to Aa^'" as a function of the lightest squark mass, minjm^^, m^^, "nij,^, ttli^}. The 
constraints from ruh, Ap, BR(6 — >■ 57) and BR(Ss f^~^f^~) have been taken into account 
(for details see Ref. [205]). 



It should be kept in mind that the size of the corrections shown in Fig. 13.201 depend 
on the assumption of the univesality of the soft SUSY-breaking parameters. It has been 
shown in Ref. [205] that lifting this universality assumption can lead to substantially larger 
contributions. As an example, for Mj^/M^ = 10 (see eq. ()1.15|) ). Aa^^ > 15 x 10"^° could 
be achieved without violating any experimental constraint. 
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3.3.3 Contributions from closed chargino/neutralino loops 



The 2-loop contributions to containing a closed chargino/neutralino loop [197] constitute a 
separately UV-finite and gauge-independent class. The corresponding diagrams were shown 
in the last line of Fig. 12.91 The chargino/neutralino two-loop contributions, a^'^^, depend 
on the mass parameters for the charginos and neutralinos /i, Mi 2, the CP-odd Higgs mass 
M4, and tan /3. 
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Figure 3.21: Comparison of the supersymmetric one-loop result aSUSY.iL (^(;Jag]-^g(j) with the 
two-loop chargino/neutralino contributions a^'^^ (dash-dotted) and the sum (full line) [197]. 
The parameters are = M2 = = M^usy, tan/5 = 50, and the sfermion mass parameters 
are set to ITeV. 



The chargino/neutralino sector does not only contribute to a^'"^^ but already to a^^^^'^'", 
so it is interesting to compare the one- and two-loop contributions. For the case that all 
masses, including the smuon and sneutrino masses, are set equal to Msusy, the one- loop and 
two- loop contributions can be trivially compared using eqs. ()2.155|1 . ()2.1(j2j) . showing that 
the two-loop contribution shifts the one-loop result by about 2%. 

However, the chargino/neutralino sector might very well be significantly lighter than the 
slepton sector of the second generation, in particular in the light of FCNC and CP-violating 
constraints, which are more easily satisfied for heavy 1st and 2nd generation sfermions. In 
Fig. 13.211 the chargino/neutralino two-loop contributions are therefore compared with the 
supersymmetric one-loop contribution a^^^'^'^'" at fixed high smuon and sneutrino masses 
Mj = 1 TeV. The other masses are again set equal, /i = M2 = Ma = MgusY- Furthermore, 
we use a large tan/3 value, tan/? = 50, which enhances the SUSY contributions to a^. 

It has been found that for Msusy ~ 400 GeV the two-loop contributions become more and 
more important. For Msusy ~ 100 GeV they amount to 50% of the one- loop contributions, 
which are suppressed by the large smuon and sneutrino masses. 
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Figure 3.22: Constraints on the MSSM parameter space in the /i-M2-plane for Ma = 
200 GeV from comparing the MSSM prediction with the data. The different regions re- 
suhing from the MSSM prediction based on qSUSY.il _|_ (jX-sl (contours with sohd border) 
and from the prediction based on a^^^^'^'" alone (dashed contours) are shown. The slepton 
mass scale (which enters only the one-loop prediction) and tan f5 are indicated for each plot. 
The contours are at (24.5, 15.5, 6.5, —2.5, —11.5, —20.5) x 10^^° corresponding to the central 
value of a;fP - a^^ecSM ^ ^24.5 ± 9.0) x 10"^° and intervals of l-5a [197]). 



The two-loop corrections have an important impact on constraints on the MSSM pa- 
rameter space obtained from confronting the MSSM prediction with the experimental value. 
This is shown in Fig. 13.221 where the regions in the /i-M2-plane resulting from the MSSM 
prediction including the two-loop correction, qSUSY.il _j_ qX.sl^ ^-^.^ compared with the cor- 
responding regions obtained by neglecting the two-loop correction, i.e. with a^^^^'^'" alone. 
The different panels correspond to different values of tan/? and the common smuon and 
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sneutrino mass Mj (the latter has an impact only on the one-loop contribution), while 
has been fixed to = 200 GeV. These parameter choices are allowed essentially in the 
entire /x-M2-plane by the current experimental constraints mentioned above, provided the 
t and b mass parameters are of 0{1 TeV). The contours drawn in Fig. 13.221 correspond to 
the Icr, 2a, ...regions around the value a^^^ — aJ^^eo.sM ^ (24.5 ± 9.0) x 10~^°, based on 
Refs. [68,72]. We find that for the investigated parameter space the SUSY prediction for 
lies mostly in the — 2 a region if fi is positive. However, the new two-loop corrections shift 
the 1 a and 2 a contours considerably. This effect is more pronounced for smaller tan j3 and 
larger M^. 

3.4 MSSM fits and constraints from existing data 

There have been many studies of the sensitivity of low-energy observables to the scale of 
super symmetry, including the precision electroweak observables [62,203,229,236-240]. Such 
analyses face the problem of the large dimensionality of the MSSM parameter space. In this 
section we discuss global fits in the unconstrained MSSM (for real parameters and using 
certain universality assumptions). Analyses in specific soft SUSY-breaking scenarios, such 
as mSUGRA, will be discussed in Sect. 01 An overview of non-super symmetric analyses of 
precision observables and resulting constraints can be found in Ref. [241]. 

The most recent global fit of the MSSM to the electroweak precison data has been 
performed in Ref. [239] (for previous analyses, see Refs. [56,236-238]). The results are shown 
in Fig. I3.23t where the predictions in the SM, the MSSM and the constrained MSSM (i.e. the 
mSUGRA scenario) are compared with the experimental data (the SUSY predictions are for 
tan/3 = 35). Fig. 13.231 shows the features discussed above: the MSSM predictions for Mw 
and (for large tan/3) {g — 2)^ are in better agreement with the data than in the SM (slight 
improvements also occur for the total width of the Z boson, Tz, and for B Xg'j). On the 
other hand, for the observables with the largest deviations between theory and experiment, 
namely Ap^ and the neutrino-nucleon cross section measured at NuTeV (the latter is not 
shown in Fig. 13.231) . the MSSM does not yield a significant improvement compared to the 
SM. The global fit in the MSSM has a lower value than in the SM. Since the MSSM fit 
has less degrees of freedom than the SM one, the overall fit probability in the MSSM is only 
slightly better than in the SM. 

3.5 Future expectations 

In this section we give a few examples of the possible physics gain obtainable with the 
anticipated improvements of the accuracies of the experimental results and the theoretical 
predictions for the precision observables (see Tab. 11.41 and the discussion in chapter 2). We 
focus here on the effects from Mw, sin^ 6'eff and nih- For a discussion of {g — 2)^ in the 
framework of the mSUGRA scenario, see chapter 4 below. 

Two examples of future prospects were already presented in Sect. 13.1.11 In Fig. I3.10l the 
SM and MSSM predictions in the mt-Mw plane are shown and compared with the current 
and future experimental precisions. Likewise, in Fig. 13.111 the results for the Mwsin^ O^s 
plane are given. It becomes aparent that the prospective improvements in the experimental 
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Figure 3.23: The predictions in the SM, the MSSM and the mSUGRA scenario (CMSSM) 
are compared with the data [239]. Deviations between theory and experiment are indicated 
in units of one standard deviation of the experimental results. 

accuracies, in particular at the ILC with GigaZ option, will provide a high sensitivity to 
deviations both from the SM and the MSSM. 

The indirect constraints on supersymmetric models from electroweak precision tests, in 
particular with GigaZ accuracy, will yield complementary information to that obtained from 
the direct observation of supersymmetric particles at the Tevatron, the LHC or the ILC (for 
a comprehensive overview on the prospects of the LHC and the ILC and the potential for 
combined analyses using LHC and ILC data, see Ref. [10]). As an example we present an 
analysis in the scalar top sector [226]. Direct information on the stop sector parameters 
and 6j: can be obtained at the ILC from the process e~^e~ titi, yielding a precision 
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of (9(1%) [242]. These direct measurements can be combined with the indirect information 
from requiring consistency of the MSSM with a precise measurement of rrih, and sin^ ^cfr- 
This is shown in Fig. 13.241 where the allowed parameter space according to measurements of 
uih, My/ and sin^^efr is displayed in the plane of the heavier stop mass, m^^^ ^ind | cos6'j| for 
the accuracies at the ILC with and without the GigaZ option and at the LHC (see Tab. II. 4j) . 
For (with an assumed central value of 180 GeV) a precision at the ILC of 1.25 GeV is 
taken [242], while for the LHC an (optimistic) uncertainty of 10% in m^^ is assumed. For 
the other parameters the following central values and prospective experimental errors have 
been used: Ma = 257 ± 10 GeV, = 263 ± 1 GeV, M2 = 150 ± 1 GeV, m-g = 496 ± 10 GeV. 
For the top-quark mass an error of 0.2 GeV has been used for GigaZ/ILC and of 2 GeV 
for the LHC. For tan f5 a lower bound of tan jS > 10 has been taken. For the future theory 
uncertainty of rrih from unknown higher-order corrections an error of 0.5 GeV has been 
assumed. The central values for Mw and sin^ ^cfr have been chosen in accordance with a 
non-zero contribution to the precision observables from SUSY loops. For the experimental 
errors at the different colliders the values given in sect. Tab. 11.41 have been used. For the 
future intrinsic theoretical uncertainties the estimates of eq. ()3.2j) have been taken. 

As one can see in Fig. I3.24[ the allowed parameter space in the m^^-\ cos9^\ plane is 
significantly reduced from the LHC to the ILC, in particular in the GigaZ scenario (i.e. 
precision measurements of Mw and sin^^efr)- Using the information on |cos%| from the 
direct measurement [242] allows an indirect determination of m^-^ with a precision of better 
than 5% in the GigaZ case. By comparing this indirect prediction for m^^ with direct 
experimental information on the mass of this particle, the MSSM could be tested at its 
quantum level in a sensitive and highly non-trivial way. 

As a further example [63] for the potential of a precise measurement of the EWPO to 
explore the effects of new physics, we show in Fig. I3.25l the predictions for Mw and sin^ 6'eff in 
the SM and the MSSM in comparison with the prospective experimental accuracy obtainable 
at the LHC and the ILC without GigaZ option (labelled as LHC/LC) and with the accuracy 
obtainable at the ILC with GigaZ option (labelled as GigaZ). For the assumed experimental 
central values of Mw and sin^ 6^^ the current central values [18] are used. For the Higgs- 
boson mass a future measured value of rrih = 115 GeV has been assumed. The MSSM 
parameters have been chosen in this example according to the reference point SPSlb [34]. 
In Fig. 13.251 the inner (blue) areas correspond to 6mf^^ = 0.1 GeV (ILC), while the outer 
(green) areas arise from 6mf^'^ = 2 GeV (LHC). For the error of Aahad we have assumed 
a future determination of 7 x 10~^. In the SM, this is the only relevant uncertainty apart 
from 6mt (the remaining effects of future intrinsic uncertainties have been neglected in this 
figure). The future experimental uncertainty of rrih is insignificant for this kind of electroweak 
precision tests. For the experimental errors on the SUSY parameters we have assumed a 5% 
uncertainty for m^-^ , m^-^ , m^^ , m^^ around their values given by SPSlb. The mixing angles in 

the t and b sectors have been left unconstrained. The mass of the CP-odd Higgs boson Ma 
is assumed to be determined to about 10%, and it is assumed that tan/? 30 ± 4.5. 

The figure shows that the improvement in 6mt from 6mt = 2 GeV to 6mt = 0.1 GeV 
strongly reduces the parametric uncertainty in the prediction for the EWPO. In the SM 
case it leads to a reduction by about a factor of 10 in the allowed parameter space of the 
Mw - sin^ Oes plane. In the MSSM case, where many additional parametric uncertainties 
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Figure 3.24: Indirect constraints on the MSSM parameter space in the m^^-\cos6^\ plane 
from measurements of ruh, Mw, sin^^es, n^'t and m^^ in view of the prospective accuracies 
for these observables at the ILC with and without GigaZ option and at the LHC The direct 
information on the mixing angle from a measurement at the ILC is indicated together with 
the corresponding indirect determination of m^^ ■ 

enter, a reduction by a factor of more than 2 is obtained in this example. The comparison 
of the theoretical prediction in both models with the GigaZ accuracy on sin^ ^eff and Mw 
illustrates how sensitively the electroweak theory will be tested via EWPO (for a compari- 
son with the current experimental errors, which are not shown in Fig. 13.251 see Fig. 13. 111) . 
The simultaneous improvement of the precision on m^, sin^ (by an order of magnitude 
compared to the situation at the LHC) and Mw (by a factor of two compared to the LHC 
case) will greatly enhance the potential for establishing effects of new physics via EWPO. 

As mentioned above, the precision observable ruh will allow to set very stringent con- 
straints on the MSSM parameters, in particular in the scalar top sector (for large values 
of tan/3 also in the scalar bottom sector). This can be crucial for determining the mixing 
angle in the scalar top sector, and (as a related quantity) the trilinear Higgs-stop coupling. 
At. If the scalar top quarks are too heavy to be directly produced at the ILC, only rather 
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Figure 3.25: The predictions for Mw and sin^^es in the SM and the MSSM (SPSlb) [63]. 
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limited information on the mixing in the stop sector will be available from the LHC [10]. 
The prospects for an indirect determination of At within the MSSM from a precision mea- 
surement of rrih are illustrated in Fig. 13.261 A precise knowledge of the parameter At turned 
out to be crucial for global fits of the MSSM to the data [53,54], which will be necessary in 
order to determine the low-energy SUSY Lagrangian parameters, and for an extrapolation 
of the results obtainable at the next generation of colliders to physics at high scales [63]. 

Fig. 13.261 shows the prediction for m/j as a function of At, where the parametric uncer- 
tainties induced by all other MSSM input parameters are taken into account according to the 
prospective experimental information on the SUSY spectrum from the LHC and the ILC in 
the SPS lb scenario [34] (see Ref. [10]). The impact of the LHC and the ILC precision on the 
top-quark mass is indicated. The sensitivity for an indirect determination of At follows from 
intersecting the MSSM prediction for rrih with the experimental value. This comparison is 
affected, however, by the intrinsic theoretical uncertainties of the ruh prediction. The effect 
of the intrinsic theoretical uncertainties is shown by two horizontal bands illustrating the 
present intrinsic uncertainty of 3 GeV and a prospective uncertainty of 0.5 GeV. While the 
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Figure 3.26: The prediction for rrih within the SPS lb scenario, assuming experimental 
information from the LHC and the ILC on the SUSY spectrum with experimental errors 
according to Ref. [10], is shown as a function of A^. The light shaded (green) band indicates 
the uncertainty induced by the experimental errors of all MSSM input parameters (except At) 
and an assumed error on the top-quark mass of 6mf^^ = 2 GeV. The dark shaded (blue) band 
shows the parametric uncertainty induced by the experimental errors of all input parameters 
for the case of Smf^^ = 0.1 GeV. The experimental error of a prospective measurement of 
ruh is shown as a horizontal band. Two further bands are shown, demonstrating the effect 
of an intrinsic theoretical uncertainty on rrih of 3 GeV (today) and 0.5 GeV (future). 

present intrinsic uncertainty on rrih would not allow to obtain a reliable indirect determina- 
tion of At, a future theoretical uncertainty of 0.5 GeV together with a precision measurement 
of rrit at the ILC would allow an indirect determination of At to better than about 10%, up 
to a sign ambiguity. The sign ambiguity can be resolved using precision measurements of 
Higgs branching ratios at the ILC, see Ref. [243]. 

Likewise, it has been shown in Ref. [243] that an indirect determination of Ma can be 
performed (investigated in the case of the SPS la scenario [34]) from Higgs boson branching 
ratio measurements at the ILC combined with a precision measurement of and information 
on the SUSY spectrum from the LHC and ILC. 
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Chapter 4 

Implications in soft SUSY-breaking 
scenarios 

The fact that no SUSY partners of the SM particles have so far been observed means that 
low-energy SUSY cannot be realized as an unbroken symmetry in nature, and SUSY models 
thus have to incorporate additional Supersymmetry breaking contributions. This is achieved 
by adding to the Lagrangian (defined by the SU(3)^ x SU(2)^ x U(l)y gauge symmetry and 
the super potential W) further terms that respect the gauge symmetry but break SUSY 
(softly, i.e. no quadratic divergences appear), so called "soft SUSY-breaking" (SSB) terms. 
The assumption made in the MSSM that the i?-parity symmetry is conserved reduces the 
amount of new soft terms allowed in the Lagrangian. 

In the previous sections the EWPO have been discussed within the unconstrained MSSM. 
In the MSSM, no further assumptions are made on the structure of the soft SUSY-breaking 
parameters, and a parametrization of all possible SUSY-breaking terms is used [244,245]. 
This gives rise to the huge number of more than 100 new parameters in addition to the SM, 
which in principle can be chosen independently of each other. A phenomenological analysis 
of the EWPO in this model in full generality would clearly be very involved, and one usually 
restricts to certain benchmark scenarios, see e.g. Refs. [32-35]. On the other hand, models 
in which all the low-energy parameters are determined in terms of a few parameters at 
the Grand Unification scale (or another high-energy scale), employing a specific soft SUSY- 
breaking scenario, provide an attractive framework for investigating SUSY phenomenology. 
The most prominent scenarios in the literature are minimal Supergravity (mSUGRA) [26,27], 
minimal Gauge Mediated SUSY Breaking (mGMSB) [28] and minimal Anomaly Mediated 
SUSY Breaking (mAMSB) [29-31]. 

The Higgs boson sector has been analyzed in all three soft SUSY-breaking scenarios, 
see Refs. [16,17,187,246,247] and references therein. For a comprehensive analysis of EWPO 
within the mSUGRA scenario see Ref. [248]. 

4.1 The soft SUSY-breaking scenarios 

The three most commonly studied soft SUSY-breaking scenarios are 
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• mSUGRA (minimal Super Gravity scenario) [26,27]: 

Apart from the SM parameters (for the experimental values of the SM input parameters 
we use Ref. [3]), 4 parameters and a sign are required to define the mSUGRA scenario: 

{ mo , mi/2 , Aq , tanp , sign(/i) } . (4.1) 

The parameter mo is a common scalar mass, mi/2 a common fermion mass and Aq a 
common trilinear couplings, all defined at the GUT scale (~ 10^^ GeV). On the other 
hand, tan (3 (the ratio of the two vacuum expectation values) and sign(/x) are defined 
at the low-energy scale. ^ 

• mGMSB (minimal Gauge Mediated SUSY-Breaking) [28]: 

An interesting alternative to mSUGRA is based on the hypothesis that the soft SUSY- 
brcaking occurs at relatively low energy scales and is mediated mainly by gauge inter- 
actions through the so-called "messenger sector" [28,249,250]. Also in this scenario, 
the low-energy phenomenology is characterized in terms of 4 parameters and a sign, 

{ Mmess, -^mess, A, tan/3, sign(//) } , (4.2) 

where Mmess is the overall messenger mass scale; Amess is a number called the mes- 
senger index, parametrizing the structure of the messenger sector; A is the universal 
soft SUSY-breaking mass scale felt by the low-energy sector. The phenomenology of 
mGMSB is characterized by the presence of a very light gravitino G with mass given 

by m3/2 = uIq — ~ ^ loo^ev ) "^-^"^ ^ [251], where is the fundamental scale 

of SUSY breaking and M'p — 2.44 x 10^^ GeV is the reduced Planck mass. Since s/F 
is typically of order 100 TeV, the G is always the LSP in the GMSB scenario. 

• mAMSB (minimal Anomaly Mediated SUSY-Breaking) [29-31]: 

In this model, SUSY breaking happens on a separate brane and is communicated to 
the visible world via the super- Weyl anomaly. The particle spectrum is determined by 
3 parameters and a sign: 

{maux, mo, tan/3, sign(//)}. (4.3) 

The overall scale of SUSY particle masses is set by maux, which is the VEV of the 
auxiliary field in the supergravity multiplet. mo is introduced as a phenomenological 
parameter to avoid negative slepton mass squares, for other approaches to this problem 
see Refs. [29,252-255]. 

^More precisely, the scenario where universaUty of the soft SUSY-breaking parameters mo, rni/2 and 
Aq at the GUT scale is assumed should be called the constrained MSSM (CMSSM). An economical way 
to ensure this universality is by gravity-mediated SUSY breaking in a minimal supergravity (mSUGRA) 
scenario, but there are other ways to validate the CMSSM assumptions. The mSUGRA scenario predicts 
in particular a relation between the gravitino mass and mo, which is not necessarily filfillcd in the CMSSM. 
For simplicity, we do not make the distinction between the CMSSM and the mSUGRA scenario but use the 
phrase "mSUGRA" for both. 
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4.2 Ap in mSUGRA, mGMSB, mAMSB 



In order to compare the prediction for Ap in three soft SUSY-breaking scenarios, a scan has 
been performed over the parameters defined in eqs. ()4.1|) - ()4.3|) . For our numerical analysis, 
the scan has been done in the following ranges: 

• mSUGRA: 



50 GeV < 
50 GeV < 
-3 TeV < 
1.5 < 



mo 
mi/2 

Ao 
tan/3 
sign /i 



< 1 TeV , 

< 1 TeV , 

< 3 TeV , 

< 60 , 
= +1. 



(4.4) 



GMSB: 



10^ GeV < 



1 < 



A 


< 


2 X 10^ GeV , 


^mess 


< 


10^ A, 


Ajness 


< 


8, 


tan P 


< 


60 , 


sign /i 




+1. 



(4.5) 



AMSB: 



20 TeV < maux < 100 TeV, 

< mo < 2 TeV, 

1.5 < tan/5 < 60, 

sign/i = +1. 



(4.6) 



For each scan point the full low-energy spectrum of the MSSM has been evaluated. It 
has been checked that the low-energy result respects the existing experimental constraints 
(for a more detailed discussion, see Ref. [246]): 

• LEP Higgs bounds: 

The results from the Higgs search at LEP have excluded a considerable part of the 
MSSM parameter space [12,13]. The results of the search for the MSSM Higgs bosons 
are usually interpreted in three different benchmark scenarios [32]. The 95% C.L. 
exclusion limit for the SM Higgs boson of Mf^ > 114.4 GeV [13] applies also for 
the lightest CP-even Higgs boson of the MSSM except for the parameter region with 
small Ma and large tan/3. In the unconstrained MSSM this bound is reduced to 
nih > 91.0 GeV [12] for Ma ^ 150 GeV and tan/5 > 8 as a consequence of a reduced 
coupling of the Higgs to the Z boson. For the CP-odd Higgs boson a lower bound 
of Ma > 91.9 GeV has been obtained [12]. In order to correctly interpolate between 
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the parameter regions where the SM lower bound^ of Mf^ ^113 GeV and the bound 
rrih ^ 91 GeV apply, we use the result for the Higgs-mass exclusion given with respect 
to the reduced ZZh coupling squared (i.e. sin^(/3 — aeff)) [256]. We have compared the 
excluded region with the theoretical prediction obtained at the two-loop level for rrih 
and sin^(/? — ttefr) for each parameter set (using = 175 GeV). 

• Experimental bounds on SUSY particle masses 

In order to restrict the allowed parameter space in the three soft SUSY-breaking sce- 
narios the current experimental constraints on their low-energy mass spectrum [3] have 
been employed. The precise values of the bounds that we have applied can be found 
in Ref. [246]. 

• Other restrictions 

As mentioned above, the top-quark mass is fixed to = 175 GeV in our analysis. 
While Ap^^ depends quadratically on at the one-loop level, the impact on Ap^^^^ 
is relatively mild. 

We briefly list here the further restrictions that we have taken into account for the 
analysis in this section. For a detailed discussion see Ref. [246]. 

— The GUT or high-energy scale parameters are taken to be real, no SUSY CP- 
violating phases are assumed. 

— In all models under consideration the i?-parity symmetry is taken to be conserved. 

— Parameter sets that do not fulfil the condition of radiative electroweak symmetry 
breaking (REWSB) are discarded (already at the level of generating the model 
parameters) . 

— Parameter sets that do not fulfil the constraints that there should be no charge or 
color breaking minima are discarded (already at the level of generating the model 
parameters) . 

— We demand that the lightest SUSY particle (LSP) is uncolored and uncharged. 
In the mGMSB scenario the LSP is always the gravitino, so this condition is 
automatically fulfilled. Within the mSUGRA and mAMSB scenario, the LSP is 
required to be the lightest neutralino. Parameter sets that result in a different 
LSP are excluded. 

— We do not apply any further cosmological constraints, i.e. we do not demand a 
relic density in the region favored by dark matter constraints [257]. 

— The scan has been stopped at high squark masses, since the contributions of heavy 
particles to Ap^^^^ decouple [97,119]. No parameter points with > 1.5 TeV 
have been considered. 

^Instead of the actual experimental lower bound, Mf^^ > 114.4 GeV [13], we use the value of 113 GeV in 
order to take into account some effect of the uncertainty in the theoretical evaluation of m/i from unknown 
higher-order corrections, which is currently estimated to be ~ 3 GeV in the unconstrained MSSM (see 

eq. 12313). 
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If a point has passed all constraints, the results for the masses and mixing angles have been 
used to determine Ap, based on the one-loop result given in eq. (j2.64|) and the SUSY two- 
loop contributions described in Sect. 12.4.21 The result is shown in Fig. 14.11 where Ap^^^'*' is 
plotted as a function of the lightest scalar top quark mass, m^^ . In general, mSUGRA allows 
smaller scalar quark masses than mGMSB and mAMSB, and correspondingly larger values of 
^pSUSY pg^j^ realized. For m^^ < 300 GeV values of Ap™^^*-^^^ ^ 7 x 10""^ can be reached. 
For larger m^-^ values all three soft SUSY-breaking scenarios result in Ap^^SY < i ^ iq-4 
(a shift in Ap^^^^ of 1 x 10~^ corresponds to shifts in and sin^ ^eff of about AM^y = 
6 MeV and Asin2 6'eff = -3 x 10-^ respectively). No part of the mSUGRA, mGMSB, or 
mAMSB parameter space that fulfils all other phenomenological constraints (see above) can 
be excluded with the current precision on the EWPO. On the other hand, for m^-^ > 500 GeV 
all three scenarios result in roughly the same prediction, i.e. it would be very challenging in 
this case to obtain information on the soft SUSY-breaking scenario with the help of Ap. 
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Figure 4.1: Ap is shown in the three soft SUSY-breaking scenarios as a function of the 
lightest scalar top quark mass. 

Using eq. ()2.58p the SUSY contribution to Ap can be translated into a shift in the 
prediction of and sin^6'cfr. For m^^ < 300 GeV the shift induced within the mSUGRA 
scenario can amount up to 

^^mSUGRA < 35 ^gy^ 1^ gi^2 ^mSUGRA| < 3 x lO'^ , (4.7) 

which corresponds roughly to one standard deviation of the current experimental uncertain- 
ties. For larger m^^ , m^^ > 500 GeV, the shifts induced in Mw and sin^ 9cs for all three soft 
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SUSY-breaking scenarios fulfil 

< g Mgy^ 1^ gi^2 ^susY| < 3.5 X 10-5 . (4.8) 

While for the possible shift in this case is about one standard deviation of the GigaZ 
precision, for sin^ 9^^ deviations of 2-3 a of the GigaZ precision could be realized. 

4.3 Prediction for ruh in mSUGRA, mGMSB, mAMSB 

We now turn to the prediction of the lightest Higgs-boson mass for the case where the 
low-energy parameters are obtained from high-scale parameters within specific soft SUSY- 
breaking scenarios. Since the low-energy parameters are connected to each other via the 
renormalization group equations, they cannot be chosen independently. This results in a 
reduction of the upper bound on rrih compared to the unconstrained MSSM. As an example, 
we show in Fig. 14.21 the allowed values of tan/5 as function of Xtop = Xf/MsusY in the 
mGMSB scenario [258]. The high-energy scan parameters are chosen as in eq. ()4.5p (but 
with both signs of /i). It can be seen that large values of tan/3, which are necessary for large 
ruh values, can only be reahzed for Xt/MsvsY between —0.3 and —1. On the other hand, the 
largest values for rrih are obtained for Xt/Msusv ~ +2 [166,167], which cannot be realized 
in the mGMSB. Similarly, also the variation of the upper bound on rrih with rrit turns out to 
be somewhat different in the soft SUSY-breaking scenarios compared to the unconstrained 
MSSM (see below). 




''top 



Figure 4.2: Allowed tan/? values as a function of Xtop = Xt/MsvsY in the mGMSB scenario 
[258]. The high-energy scan parameters are chosen as in eq. ()4.5j) (but with both signs of fi). 

In the following we refer to the results of Ref. [187], which are in agreement with the 
previous results in Refs. [17,246], but use the most recent experimental value of the top- 
quark mass. In Tab. 14.11 the maximum values of rUh for = 178.0 GeV in mSUGRA, 
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mGMSB and mAMSB are compared. In order to have comparable numbers an upper limit 
on the scalar top masses in all scenarios has been chosen, ^m^^ m^^ < 2 TeV. No theoretical 
uncertainties are included. One can see that all three scenarios result in significantly lower 
maximum rrih values than the uncontrained MSSM, where masses up to ~ 138 GeV can 
be reahzed for Mgusy ~ 2 TeV and = 178.0 GeV (see Figs. 13.15} 13. 18^ . The variation 
of this maximum value with is also shown. In the unconstrained MSSM one has 
Smh/Snit ~ 1 [188]. In the mSUGRA, mGMSB and mAMSB scenarios this is reduced down 
to ~ 0.58-0.7. 





maximum rrih [GeV] 


Sruh/Smt 


mSUGRA 


129.0 


0.65 


mGMSB 


123.7 


0.70 


mAMSB 


124.6 


0.58 



Table 4.1: The maximum rrih values (for rrit = 178.0 GeV and ^m^^ m^^ < 2 TeV) and 
the variation of this maximum value with rrit are shown in the three soft SUSY-breaking 
scenarios. No theoretical uncertainties are included. See Refs. [17,187,246]. 

These results have an interesting consequence for the Higgs search at the Tevatron. The 
Tevatron has the potential to exclude a SM-like Higgs boson with a mass of Mf^^ ^ 130 GeV 
with an integrated luminosity of 4-8 fb~^ [14] per experiment (and it will furthermore reduce 
the experimental error on rrit). Since the coupling of the lightest CP-even Higgs boson to 
gauge bosons is close to the SM value for essentially all the parameter space of the three soft 
SUSY-breaking scenarios [246], the Tevatron should either observe an excess of Higgs-like 
events over the background expectation or rule out the mSUGRA, the mGMSB and the 
mAMSB scenarios. 

4.4 EWPO in mSUGRA 

In this section we review the prediction for Myy, sin^^efr, the lightest Higgs boson mass, 
the anomalous magnetic moment of the muon, = {g — 2)^/2, and BR(6 — *■ 57) within 
the mSUGRA scenario, taking into account constraints on the cold dark matter (CDM) 
relic density from WMAP and other cosmological data [257]. More details can be found 
in Ref. [248]. The results have been obtained by scanning the universal soft supersymmetry- 
breaking gaugino mass mi/2 and scalar mass mo for different representative values of tan/3 
and the trilinear soft supersymmetry-breaking parameter Aq. The sign of the supersymmetric 
Higgs parameter n has been chosen to be positive. 

We require the cosmological relic density f^^^^ to the neutralino LSP to fall into the 
range 

0.094 < Vt^K^ < 0.124 . (4.9) 

Lower values of f^^^^ would be allowed if not all the cosmological dark matter is composed 
of neutralinos. However, larger values of are excluded by cosmology. The CDM con- 
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straints have the effect within the mSUGRA scenario, assuming that the dark matter consists 
largely of neutralinos, of restricting itlq to very narrow allowed strips for any specific choice of 
Aq, tan/3 and the sign of fi [260,261]. Thus, the dimensionality of the mSUGRA scenario is 
effectively reduced, and one may explore SUSY phenomenology along these "WMAP strips" . 
We furthermore take into account the constraints on the parameter space from the direct 
search for supersymmetric particles [3] and Higgs bosons [12,13]. 

For tan/3 two values have been chosen, tan/3 = 10,50, representing values in the lower 
and the upper part of the (experimentally and theoretically) allowed parameter space. For 
the GUT-scale parameter Aq five different values have been investigated (below also a scan 
over Aq is performed), Aq = (—2, —1, 0, 1, 2) xmi/2, in order to cover the allowed parameter 
space. The top-quark mass has been fixed to rrit = 178 GeV. Since the results are analyzed 
along the WMAP strips, they are given as a function of mi/2- The corresponding itlq values 
(for fixed Aq and tan/3) follow from the CDM constraint. The non-excluded values for mi/2 
start at around mi/2 ~ 200 GeV for both values of tan/3. While for tan (3 = 10 mi/2 is 
restricted by the CDM constraint to be mi/2 ^ 900 GeV, for tan/3 = 50 the allowed values 
exceed mi/2 ^ 1500 GeV. 

We start with the prediction for Mw The evaluation is based on the corrections described 
in Sect. 12.51 We display in Fig. I4.3l the mSUGRA prediction for Myy and compare it with the 
present measurement (solid lines) and a possible future determination with GigaZ (dashed 
lines). Panel (a) shows the values of Mw obtained with tan/3 = 10 and |y4o| < 2, and 
panel (b) shows the same for tan/3 = 50. It is striking that the present central value of 
Mw (for both values of tan/3) favours low values of mi/2 ~ 200-300 GeV, though values as 
large as 800 GeV are allowed at the 1-a level, and essentially all values of mi/2 are allowed 
at the 90% confidence level. The GigaZ determination of Mw might be able to determine 
indirectly a low value of mi/2 with an accuracy of ±50 GeV, but even the GigaZ precision 
would still be insufficient to determine mi/2 accurately if mi/2 ^ 600 GeV (in accordance 
with the discussion in Sect. 14. 2() . 

The situtation is similar for the prediction of sin^ ^eff shown in Fig. 14.41 The results are 
based on the corrections described in Sect. 12.61 and are given for the same values of Aq and 
tan/3 as in Fig. 14.31 As in the case of Mw, low values of mi/2 are also favoured by sin^ 9cs- 
The present central value prefers mi/2 = 300-500 GeV, but the 1-a range extends beyond 
1500 GeV (depending on Aq), and all values of mi/2 are allowed at the 90% confidence level. 
The GigaZ precision on sin^ 9es would be able to determine mi/2 indirectly with even greater 
accuracy than Mw at low mi/2, but would also be insufficient if mi/2 ^ 700 GeV. 

Next the prediction of within mSUGRA is analyzed. The evaluation is based on 
the full one-loop result [198], the corresponding QED two-loop corrections [204] and the 
two-loop corrections from the closed SM fermion/sfermion loops [205]. The very recent 
two-loop corrections of Ref. [197] have been included via an approximation formula. For 
older evaluations of within mSUGRA (mostly based on the full one-loop result and the 
corresponding QED corrections), see Refs. [200-203]. 

As seen in Fig. 14.51 the mSUGRA prediction for is almost independent of Aq for 
tan/3 = 10, but substantial variations are possible for tan/3 = 50, except at very large mi/2. 
In the case tan/3 = 10, mi/2 ~ 200-400 GeV is again favoured at the il-cr level, but this 
preferred range shifts up to 400 to 800 GeV if tan/? = 50, depending on the value of Aq. 
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Figure 4.3: The mSUGRA prediction for M]^ as a function of mi/2 along the WMAP strips 
for (a) tan/3 = 10 and (b) tan (3 = 50 for various Aq values [248]. In each panel, the 
centre (solid) line is the present central experimental value, and the (solid) outer lines show 
the current ±1-0" range. The dashed lines correspond to the anticipated GigaZ accuracy, 
assuming the same central value. 
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Figure 4.4: The mSUGRA prediction for sin^ ^cfr as a function of mi/2 along the WMAP 
strips for (a) tan (3 = 10 and (b) tan /3 = 50 for various Aq values [248] . In each panel, the 
centre (solid) line is the present central experimental value, and the (solid) outer lines show 
the current il-cr range. The dashed lines correspond to the anticipated GigaZ accuracy, 
assuming the same central value. 
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Figure 4.5: The mSUGRA prediction for Aa^ as a function of mi/2 along the WMAP strips 
for tan/3 = 10,50 and different Aq values [248]. The central (solid) line is the central value 
of the present discrepancy between experiment and the SM value evaluated using e+e~ data, 
and the other solid (dotted) lines show the current ±1(2)-(T ranges, see eq. ()2.147|) . 

For the two tan/5 values the requirement of agreement of the mSUGRA prediction with the 
experimental data at the 95% CL. restricts mi/2 to 

tan/? = 10 : 200 GeV < mi/2 < 600 GeV , (4.10) 
tan/? = 50 : 350 GeV < mi/2 < 1100 GeV . (4.11) 

Now we turn to the decay b s'-f. Since this decay occurs at the loop level in the SM, 
the MSSM contribution might be of similar magnitude. The most up-to-date theoretical 
estimate of the SM contribution to the branching ratio is [262] 

BR(6 S7) = (3.70 ± 0.30) x 10-^ (4.12) 

where the calculations have been carried out completely to NLO in the MS renormalization 
scheme, and the error is dominated by higher-order QCD uncertainties. However, the error 
estimate for BR(6 — > 57) is still under debate, see e.g. Ref. [235]. The MSSM evaluation 
shown below is based on Refs. [262,263]. 

For comparison, the present experimental value estimated by the Heavy Flavour Aver- 
aging Group (HFAG) is [264] 

BR(6 ^ S7) = (SMt'ofs) X 10-^ (4.13) 

where the error includes an uncertainty due to the decay spectrum, as well as the statistical 
error. The very good agreement between eq. ()4.13p and the SM prediction eq. ()4.12|) imposes 
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important constraints on the MSSM. The uncertainty range shown in Fig. 14.61 combines 
hnearly the current experimental error and the present theoretical uncertainty in the SM 
prediction. Since the mSUGRA corrections are generally smaller for smaller tan/3, even 
values of mi/2 as low as ~ 200 GeV would be allowed at the 90% confidence level if tan/5 = 
10, whereas mi/2 ^ 400 GeV would be required if tan/5 = 50. These limits are very sensitive 
to Aq, and, assuming that in the future the experimental and theoretical uncertainty in 
BR(6 s'-f) can be reduced by a factor ~ 3, the combination of BR(6 — > 57) with the other 
precision observables might be able, in principle, to constrain Aq significantly. 
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Figure 4.6: The mSUGRA predictions for BR(6 — > s'-f) as a function of mi/2 along the WMAP 
strips for (a) tan j3 = 10 and (b) tan /5 = 50 and various choices of Aq. The uncertainty shown 
combines linearly the current experimental error and the present theoretical uncertainty in 
the SM prediction. The central (solid) line indicates the current experimental central value, 
and the other solid (dotted) lines show the current ±1(2)-(T ranges [248]. 



Finally we present results for the lightest Higgs boson mass in the CDM allowed strips 
of the mSUGRA parameter space. In Fig. 14.71 we show the results for rrih. A hypothetical 
measurement at mh = 120 GeV is shown. Since the experimental error at the ILC will 
be smaller than the prospective theory uncertainties (see Sect. 12.7."^ . we display the effect 
of the current and future intrinsic uncertainties. In addition, a more optimistic value of 
Aruh = 200 MeV is also shown. The figure clearly illustrates the high sensitivity of this 
electroweak precision observable to variations of the supersymmetric parameters (detailed 
results for Higgs boson phenomenology in mSUGRA can be found in Ref. [16,17,246,247]). 
The comparison between the measured value of rrih and a precise theory prediction will allow 
to set tight constraints on the allowed parameter space of mi/2 and Aq. 
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Figure 4.7: The mSUGRA predictions for mh as functions of mi/2 with (a) tan/5 = 10 and 
(b) tan/? = 50 for various Aq [248]. A hypothetical experimental value is shown, namely 
nih = 120 GeV. We display an optimistic anticipated theory uncertainty of ±0.2 GeV, as 
well as a more realistic theory uncertainty of ±0.5 GeV and the current theory uncertainty 
of ±3 GeV. 



4.5 Fits in mSUGRA 

The results for EWPO presented in the last section have been used to perform a fit for 
the mSUGRA parameter space with CDM constraints [248]. We first review the fit using 
the currently existing data on Myj/, sin^6'efr, and BR(6 57). Secondly, we show the 
precision that can be obtained in the future, using improved measurements of the EWPO 
and including also the measurement as well as the measurement of Higgs boson branching 
ratios. More details can be found in Ref. [248] 



4.5.1 Present situation 

We now investigate the combined sensitivity of the four low-energy observables for which 
experimental measurements exist at present, namely M^^, sin^^efr, (5'~2)^ and BR(6 — >■ 57). 
We begin with an analysis of the sensitivity to mi/2 moving along the WMAP strips with fixed 
values of Aq and tan/?. The experimental uncertainties, the intrinsic errors from unknown 
higher-order corrections and the parametric uncertainties have been added quadratically, 
except for BR(6 57), where they have been added linearly. Assuming that the four 
observables are uncorrelated, a fit has been performed with 



N / T-.ovr^ T-)theo 



^(flr^) . (4.14) 



n=l 



Here R^^^ denotes the experimental central value of the nth observable, so that = 4 for 
the set of observables included in this fit, E}^'^" is the corresponding mSUGRA prediction 
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and cr„ denotes the combined error, as specified above. 

The results are shown in Fig. 14.81 for tan/3 = 10 and tan/3 = 50. They indicate that, 
already at the present level of experimental accuracies, the electroweak precision observables 
combined with the WMAP constraint provide a sensitive probe of the mSUGRA scenario, 
yielding interesting information about its parameter space. For tan/3 = 10, mSUGRA 
provides a very good description of the data, resulting in a remarkably small minimum 
value. The fit shows a clear preference for relatively small values of mi/2, with a best-fit 
value of about mi/2 = 300 GeV. The best fit is obtained for Aq < 0, while positive values of 
^0 result in a somewhat lower fit quality. The fit yields an upper bound on mi/2 of about 
600 GeV at the 90% C.L. (corresponding to Ax^ < 4.61). The mass spectrum favored at 
the 90% C.L. contains many light states that should be accessible at the LHC and the ILC, 
offering good prospects of the direct detection of SUSY. 

For tan/3 = 50 the overall fit quality is worse than for tan/3 = 10, and the sensitivity to 
mi/2 from the precision observables is lower. This is related to the fact that, whereas 
and sin^ O^s prefer small values of mi/2 also for tan/3 = 50, as seen in Figs. 14.31 and 14.41 the 
CMSSM predictions for {g — 2)^ and BR(6 — 57) for high tan (3 are in better agreement 
with the data for larger mi/2 values, as seen in Figs. 14.51 and 14.61 Also in this case the best 
fit is obtained for negative values of Aq, but the preferred values for mi/2 are 200-300 GeV 
higher than for tan/3 = 10. The mass spectrum favored at the 90% C.L. is heavier than 
for tan/3 = 10. However, still several SUSY particles should be accessible at the ILC. Since 
colored SUSY particles should be within the kinematic reach of the LHC, also in this case 
there are good prospects of the direct detection of SUSY. 

We now turn to the results obtained from a scan over the mi/2^^0 parameter plane. 
Fig. 14.91 shows the CDM-allowed regions in the mi/2-Ao plane for tan (3 = 10 and tan jS = 50. 
The current best-fit values obtained via fits for tan /3 = 10 and tan /3 = 50 are indicated. 
The coloured regions around the best-fit values correspond to the 68% and 90% C.L. regions 
(corresponding to Ax^ < 2.30,4.61, respectively). 

For tan/3 = 10 (upper plot of Fig. 14. 9p . the precision data yield sensitive constraints on 
the available parameter space for mi/2 within the WMAP-allowed region. The precision data 
are less sensitive to Aq. The 90% C.L. region contains all the WMAP-allowed Aq values in 
this region of mi/2 values. As expected from the discussion above, the best fit is obtained 
for negative Aq and relatively small values of mi/2- At the 68% C.L., the fit yields an upper 
bound on mi/2 of about 450 GeV. This bound is weakened to about 600 GeV at the 90% 
C.L.3 

As discussed above, the overall fit quality is worse for tan (3 = 50, and the sensitivity 
to mi/2 is less pronounced. This is demonstrated in the lower plot of Fig. 14.91 which shows 
the result of the fit in the mi/2^^0 plane for tan/3 = 50. The best fit is obtained for 
mi/2 ~ 500 GeV and negative Aq. The upper bound on mi/2 increases to nearly 1 TeV at 
the 68% C.L. 

■^A preference for relatively small values of mi/2 within the mSUGRA has also been noticed in Ref. [239], 
where only (g — 2)^ and BR(6 —^ 57) had been analyzed. 
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Figure 4.8: The results of fits based on the current experimental results for the precision 
observables M^, sin^^efr, (g — '^)fi and BR(6 — > 57) are shown as functions of mi/2 in the 
mSUGRA parameter space with CDM constraints for different values of Aq [248]. The upper 
plot shows the results for tan j3 = 10, and the lower plot shows the case tan j3 = 50. 
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Figure 4.9: The results of fits for tan/? = 10 (upper plot) and tan/? = 50 (lower plot) 
based on the current experimental results for the precision observables M^r, sin^ 6'efr, — 
and BR(6 s'-f) are shown in the planes of the mSUGRA scenario with the WMAP 

constraint [248]. The best-fit points are indicated, and the coloured regions correspond to 
the 68% and 90% CL. regions, respectively. 
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4.5.2 Future expectations 

We now investigate the combined sensitivity of the precision observables Mw, sin^^^efi, 
{g - 2)^, BR(6 — > 57), rrih and the ratio BR(/i bb)/BR{h WW*) in the ■mi/2-Ao plane 
of the mSUGRA scenario using ILC (and GigaZ) accuracies. For {g — 2)^ we assume a 
reduction of the error by two, for BR(6 57) by a factor of three. At the ILC with 
^/s = 1 TeV a measurement of BR(/i bb)/BR{h — * WW*) with an accuracy of ~ 1.5% 
can be envisaged [265]. Fig. 14.101 shows the fit results for tan/5 = 10, while Fig. 14.111 shows 
the tan /? = 50 case. 

In each figure we show two plots, where the WMAP-allowed region and the best-fit point 
according to the current situation (see Fig. 14. 9p are indicated. In both plots two further 
hypothetical future 'best-fit' points have been chosen for illustration. For all the 'best-fit' 
points, the assumed central experimental values of the observables have been chosen such 
that they precisely coincide with the 'best-fit' points^. The coloured regions correspond to 
the 68% and 90% CL. regions around each of the 'best-fit' points according to the ILC 
accuracies. 

The comparison of Figs. 14.101 I^TTl with the result of the current fit. Fig. 14.91 shows that 
the ILC experimental precision will lead to a drastic improvement in the sensitivity to mi/2 
and Aq from comparing precision data with the mSUGRA predictions. For the best-fit values 
of the current fits for tan j3 = 10 and tan j3 = 50, the ILC precision would allow one to narrow 
down the allowed mSUGRA parameter space to very small regions in the plane. 
The comparison of these indirect predictions for mi/2 and Aq with the information from 
the direct detection of supersymmetric particles would provide a stringent test of the model 
at the loop level. A discrepancy could indicate that supersymmetry is realised in a more 
complicated way than assumed in mSUGRA. 

The additional hypothetical 'best-fit' points shown in Figs. 14.1'm 14.111 illustrate the in- 
direct sensitivity to the mSUGRA parameters in scenarios where the precision observables 
prefer larger values of mi/2- Because of the decoupling property of supersymmetric theories, 
the indirect constraints become weaker for increasing mi/2- 

For tan/9 = 10, we have investigated hypothetical 'best-fit' values for mi/2 of 500 GeV, 
700 GeV (for Ao > and Aq < 0) and 900 GeV. For mi/2 = 500 GeV, the 90% CL. 
region in the mi/2-Ao plane is significantly larger than for the current best-fit value of 
mi/2 ~ 300 GeV, but interesting limits can still be set on both mi/2 and Aq. For mi/2 = 
700 GeV and mi/2 = 900 GeV, the 90% CL. region extends up to the boundary of the 
WMAP-allowed parameter space for mi/2- Even for these large values of mi/2, however, the 
precision observables (in particular the observables in the Higgs sector) still allow one to 
constrain Aq. 



^It was checked explicitly that assuming future experimental values of the observables with values dis- 
tributed statistically around the present 'best-fit' points with the estimated future errors does not degrade 
significantly the qualities of the fits. 
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Figure 4.10: The results of a fit based on the prospective experimental accuracies for the 
precision observables Mw, sin^^efr, {g — 2)^, BR(6 — > 57), rrih and Higgs branching ratios 
at the ILC are shown in the plane of the mSUGRA with WMAP constraints for 

tan/3 = 10 [248]. In both plots the WMAP-allowed region and the best-fit point according 
to the current situation (see Fig. 14. 9|) are indicated. In both plots two further hypothetical 
future 'best-fit' values have been chosen for illustration. The coloured regions correspond to 
the 68% and 90% CL. regions according to the ILC accuracies. 
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Figure 4.11: The results of a fit based on the prospective experimental accuracies for the 
precision observables Mw, sin^ 9es, {g — 2)^, BR(6 57), rrih and Higgs branching ratios at 
the ILC are shown in the mi/2-Ao plane of the mSUGRA scenario with WMAP constraints 
for tan/3 = 50 [248]. In both plots the WMAP-allowed region and the best-fit point for 
tan /3 = 50 according to the current situation (see Fig. 14. 9p are indicated. In both plots two 
further hypothetical future 'best-fit' values have been chosen for illustration. The coloured 
regions correspond to the 68% and 90% CL. regions according to the ILC accuracies. 
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For tan/3 = 50, where the WMAP- allowed region extends up to much higher values 
of mi/2, we find that for a 'best-fit' value of mi/2 as large as 1 TeV the precision data 
still allow one to establish an upper bound on mi/2 within the CDM-allowed region. This 
indirect sensitivity to mi/2 could give important hints for supersymmetry searches at high- 
energy colliders. For 'best-fit' values of mi/2 in excess of 1.5 TeV, on the other hand, 
the indirect effects of heavy sparticles become so small that they are difficult to resolve even 
with ILC accuracies. To conclude, the indirect sensitivity from the measurement of precision 
obscrvablcs at the ILC have a potential even to exceed the direct search reach of both the 
LHC and ILC. 
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Chapter 5 
Conclusions 



An overview of the current status of precision tests of supersymmetry has been given, and 
future prospects have been discussed. We have mainly focused on the W boson mass, Mw, 
the effective leptonic weak mixing angle, sin^^eflf, the anomalous magnetic moment of the 
muon, {g — 2)^, and the lightest CT'-even MSSM Higgs boson mass, irih, but constraints 
from b physics, direct collider searches and cosmological data have also been included in the 
discussion. 

Precise experimental data are available for M^, sin^ 9^^ and {g—2)n, while nih is expected 
to become a precision observable if a supersymmetric Higgs sector is realised in nature. 
Confronting the high experimental precision with the theory predictions provides sensitivity 
to quantum corrections of the theory, where the whole structure of the model enters. This 
allows to set indirect constraints on the properties of particles even if they are too heavy 
to be produced directly. In order to exploit the experimental precision, the theoretical 
predictions for the electroweak precision observables in supersymmetry (or in other models 
that are confronted with the data) should be at least at the same level of accuracy. Ideally, 
the remaining theoretical uncertainties should be so small that they are negligible compared 
to the experimental errors. Sophisticated higher-order calculations are necessary in order 
to match this demand, and a considerable effort will be required for keeping up with the 
prospective improvements of the experimental accuracies in future experiments. 

We have briefly discussed the necessary ingredients of higher-order calculations in super- 
symmetry, focusing in particular on regularisation and renormalization, and have pointed 
out important differences compared to the case of the SM. The large number of parameters 
in the MSSM, most of which are not directly related to any particular physical observable, 
and the relations imposed by the underlying symmetry make it quite involved to formulate 
a coherent and easily applicable renormalization prescription for the whole MSSM. Differ- 
ent prescriptions exist in the literature for various sectors of the MSSM, but no common 
standard has emerged yet. 

The current status of the theoretical predictions for the most important precision ob- 
servables has been revieved, and estimates of the remaining theoretical uncertainties from 
unknown higher-order corrections and from the experimental errors of the SM input pa- 
rameters have been given. The theoretical predictions have then been compared with the 
current experimental results (in the case of nih the MSSM prediction has been confronted 
with the exclusion bounds from the Higgs search at LEP) . The resulting constraints on the 
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MSSM parameter space have been analyzed. We have investigated how well the MSSM 
describes the data and whether the data give some preference for the MSSM as compared 
to the SM. This has been analyzed both for the unconstrained MSSM and for specific soft 
SUSY-breaking scenarios. The mSUGRA scenario, characterised by four parameters and 
a sign, can still simultaneously satisfy the constraints from the electroweak precision data, 
direct collider searches and the stringent bounds on cold dark matter in the universe from 
WMAP and other cosmological data. It turns out that the mSUGRA scenario with cosmo- 
logical constraints in fact yields a very good fit to the data. The fit results indicate a clear 
preference for a relatively light mass scale of the SUSY particles, offering good prospects for 
direct SUSY searches at the LHC and at the ILC. 

We have investigated future prospects of electroweak precision tests of supersymmetric 
models. Anticipated improvements in the experimental precision have been discussed in 
view of the LHC and the ILC, and the prospects for a further reduction of the theoretical 
uncertainties have been analyzed. Based on these estimates of future experimental and 
theoretical precisions, we find that the sensitivity of the precision tests will improve very 
significantly, leading to stringent constraints on the MSSM parameter space (and on any 
other conceivable model of new physics). If supersymmetric particles are discovered at the 
next generation of colliders, the combination of information from the direct observation of 
SUSY particles and the indirect information from electroweak precision obscrvables will allow 
very powerful tests of the model. This can lead to a discrimination between the minimal and 
non-minimal models, a distinction between different SUSY-breaking scenarios, and indirect 
predictions for parameters or particle masses that are not directly experimentally accessible. 
These consistency tests at the quantum level using all available experimental information 
will be crucial in the quest to extrapolate the results of the next generation of colliders to 
physics at high scales. 
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Appendix A 
Loop integrals 



In this appendix we present the loop integrals needed for the two-loop evaluation of the 
SUSY contributions to the EWPO. D denotes the space-time dimension and 5 = ^{A — D). 
In the following formulas we neglected the terms proportional to 7^ — ln(47r //^) , which are 
connected to the divergent parts. They always cancel for physical observables. 

The analytical formulas for Aq and Bq are taken from Ref. [141], T134 and T234/ are taken 
from Ref. [142], the other integrals can be found in Ref. [126]. The notation for the integrals 
is as in Ref. [214]. 
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Here contains a small imaginari part, ze, e > 0. 
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(A.24) 



where ri and r2 are given by eq. ()A.7I 



A. 5 T123/4 

The following formula, looking at the series expansion in 1/6, are correct up to 0{6^). 



Ti23'4{ml,ml,0,ml) = Ti23'4("^?, ^4, m^, 0) = (A.25) 
(1 - 2(ln(p2) - tn)5 + ^ (21n(p2) - ^nY S'^ l±l(B{p',ml,ml)Y 



1 

+2 



' Xi4ln^(Xi4) X4iln^(x4i) X , ^/ N 

- (1-X,4)^ - (1-0:4,)^ " ^^""^ + 
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with the following functions 
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(A.27) 
(A.28) 

(A.29) 
(A.30) 

(A.31) 



A. 6 rii23/4 

The following formula, looking at the series expansion in 1/5, correct up to 



Tii23'4(m?, ml 0, m?) = | (l - 2(ln(p') - iT^)5 + ^ (2 ln(p2) - mf 



(-1 - 6)B{p^, ml, ml)B'{p^, mj, ml) 



y-p -p J 



ln(a;i4) 1 + ln(a;i4) 0:41 ln(a;4i) 1 + 2:41 ln(a;4i) 



1 — xu 1 — Xu 



IV 



2 mi 



777^ TJlP' 



1 — X41 1 — X41 
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(A.32) 



A. 7 T123/45 

The following formula, looking at the series expansion in 1/5, correct up to 0{S^). 



^ (F(l) + F(m4i) - F(xi4) - F(x4i)) (A.33) 



This function is finite in the limit 5^0. 
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Appendix B 

Input parameters and benchmark 
scenarios 



For our numerical results, the following values 
other quark and lepton masses are negligible): 

Gf = 1.16639 X 10"^ 
Mw = 80.450 GeV, 
Mz = 91.1875 GeV, 
Tz = 2.4952 GeV, 



of the SM parameters have been used (all 
rur = 1.777 GeV, 

mt = 174.3 GeV, , , 

B.l 

nib = 4.25 GeV, ^ ' 

= 1.5 GeV. 



The predictions for the observables in this report are in some cases expressed in terms of 
running bottom- and top-quark masses in order to absorb QCD corrections. The numerical 
values of these running masses differ from the pole masses given in eq. ()B.1|1 . 

For our numerical evaluation we often refer to four benchmark scenarios that have been 
defined in Ref. [33] for Higgs boson searches at hadron colliders and beyond. The four 
benchmark scenarios are (more details can be found in Ref. [33]) 

• the "m™'^^" scenario, which yields a maximum value of Mh for given Ma and tan f3, 

mt = 174.3 GeV, Mgusy = 1 TeV, /i = 200 GeV, M2 = 200 GeV, 
Xt = 2MsusY, Ar = Ab = At, m^ = 0.8MsusY, 

• the "no-mixing" scenario, with no mixing in the t sector, 

mt = 174.3 GeV, Mgusv = 2 TeV, /i = 200 GeV, M2 = 200 GeV, 

Xt = 0, Ar = Ab = At, mg = 0.8 MsvsY, (B.3) 

• the "gluophobic-Higgs" scenario, with a suppressed ggh coupling, 

mt = 174.3 GeV, Mgusv = 350 GeV, /i = 300 GeV, M2 = 300 GeV, 
Xt = -750 GeV, = Ab = At, m-g = 500 GeV , (B.4) 
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• the "small- CKeff" scenario, with possibly reduced decay rates for h ^ bb and h t'^t , 

mt = 174.3 GeV, MgusY = 800 GeV, fj, = 2.5 Msusy, M2 = 500 GeV, 
Xt = -1100 GeV, Ar^Ab^ At, m~g = 500 GeV . (B.5) 

As explained above, for the sake of simplicity, Msusy is chosen as a common soft SUSY- 
breaking parameter for all three generations. 
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